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PREFACE 


Le  but  de  ce  document  est  essentiellement  de  rassembler  sous  forme  didactique  tes  connaissances 
actuelles  sur  les  mithodes  experimentales  de  determination  des  coefficients  airodynamiques  necessaires 
aux  etudes  de  mecanique  du  vol. 

A  ma  demande,  I'auteur  a  accepte  d'eiargir  quelque  peu  son  sujet,  de  maniere  : 

-  d'une  part  a  exprimer  les  donnees  theoriques  fondamentales  en  langage  matriciel,  ce  qui, 

malgre  les  avantages  evidents  de  clarte  et  de  concision,  6tait  rarement  le  cas  dans  les  documents 
deja  publies  sur  ce  theme  ; 

-  d'autre  part,  a  souligner  les  lacunes  des  methodes  actuelles  lorsqu'il  s’agit  d’analyser  des  ecoule- 

ments  fortement  non  lineaires. 

J'espere  que,  sous  cette  forme,  cet  ouvrage  pourra  servir  a  I'accomplissement  des  proves  indispen- 
sables  et  a  I'etablissement  d'un  langage  commun  entre  specialistes  de  la  Mecanique  du  vol  et  adrodyna- 
miciens. 


P.  CARRIERE 

Haut  Conseiller  Scientifique  a  I'ONERA 


Juin  1979 

Ingenieur  General  P.  CARRIERE,  Membra  du  Fluid  dynamic  Panel,  Editaur. 


DETERMINATION  DANS  LES  INSTALLATIONS  AU  SOL 
DES  PARAMETRES  AERODYNAMIQUES  DE  STABILITE 
DES  AERONEFS 


RESUME 

Un  apergu  da  l'dtat  actual  das  moyens  d ' experimentation  au  sol,  servant  A  la  determination  das  paramAtres 
a&rodynamiques  de  stability  das  a^ronefs,  a  i6gag4  das  comptes— rendus  das  reunions  FDP  at  FMP  AGARD, 
tenues  pendant  cen  quatre  demiferes  anndes.  Un  expose  critique  das  methodes  axperimantalas  d' analyse 
das  phenomena,  a4rodynauiques  instationnaires  met  en  evidence  las  lacunas  das  methodes  at  moyens  actuals 
notamment  en  ce  rui  concerns  les  informations  necessairas  A  la  mise  en  oeuvre  correcte  das  simulateurs 
da  vol  dans  les  cas  do  compcrtements  non  lineairas  (decollemants,  ondes  de  choc).  Das  suggestions 
sur  das  orientations  possibles  A  dormer  aux  recherches  pour  combler  cea  lacunas,  sont  presentees. 


DETERMINATION  IN  GROUND  FACILITIES 
OF  AERODYNAMIC  STABILITY  PARAMETERS  OF  AIRCRAFT 


ABSTRACT 

A  survey  of  the  present  state  of  experimental  ground  facilities  fcr  determining  the  aerodynamic  stability 
parameters  of  aircraft  has  been  deduced  from  the  proceedings  of  the  AOARD/fDP  and  FMP  meetings  of  the 
last  four  years.  A  critical  study  of  the  experimental  methods  for  analyzing  unsteady  aerodynamic  pheno¬ 
mena  brings  to  light  the  insufficiencies  of  present-day  methods  and  means,  especially  as  regards  infor¬ 
mation  necessary  for  the  correct  implementation  of  flight  simulators  in  case  of  non-linear  behaviour 
(separation,  shock  wave).  Suggestions  are  presented  on  the  possible  orientations  of  research  for 
filling  these  gaps. 
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NOTATIONS 


1  -  NOTATIONS  NORMALISEES 

!Hm  sent  conforr.es  aux  normes  ISO  R  1151,  R  1152,  R  1153,  et  au  projet  de  norms  ISO/DIS  2765.  [1*6] 
TRmans  . 

X  y  z  trlAdre  lid  A  l'avlon,  dont  l'orlglne  0  est  usuellement  le  centre  de  gravitA  de 

l'avlon.  II  est  constltuA  par  I'axa  longitudinal  Ox,  l'axe  transversal  Oy  et  l'axe 
normal  Ox  (fig.  1),  oe  demler  dtant  orlentd  vers  le  ventre  de  l'avlon. 

Xaya*a  trlAdre  adrodynamlque  ayant  usuellement  la  mSme  orlglne  que  le  prdoddent.  L'axe  x& 

est  confondu  en  direction  et  sens  aveo  la  Vitesse  air,  l'axe  z  est  dans  le  plan 
de  symAtrie  (fig.  1). 

W  trlAdre  terrestre,  appelA  trlAdre  normal  terrestre  lore que  l'axe  zQ  eat  vertioal. 

0  l’axe  zQ  est  orientA  vers  le  bas,  L'orlglne  est  soit  un  point  fixe  par  rapport  A  la 

terre,  choisl  sulvant  les  besolns,  solt  le  centre  de  gravltd  G  ;  dans  ce  demler  cas 
et  si  de  plus  l'axe  zQ  est  vertioal,  le  trlAdre  est  appeld  trlAdre  normal  terrestre 
portd  par  l'avlon  (fig.  2  et  3). 

Voir  aux  notations  vectorlelles  les  symboles  des  veoteurs  unltalres  de  oes  trlAdres. 


Remarques  :  d'autres  trlAdres,  non  normalises ,  sont  couramment  utilises  : 

1)  les  axes  de  "stabilitd”  lies  A  1'aAronef,  oonfondus  avec  le  triAdre  adrodynamique  lorsquo  1’aAronef 
est  dans  sn  position  de  rAgime, 

2)  le  triedre  "aArobalistique’’,dAsigne  souvent  par  xyz,  comprend  l’axe  longitudinal  X  et  est  anime 
de  la  vitesse  angulaire  -p  pur  rupport  au  triedre  lie  xyz  (non  rolling  axes). 


ANGLES 

a)  DAfinisaant  la  position  angulaire  du  vecteur  vitesse-air  par  rapport  au  triAdre-avion  (fig.  l)  : 

P  angle  de  dArapage  ;  par  convention  ^  ^  "f  * 

a  angle  d'incidence  (ou  d’attaque)  par  convention  -  s  <  a  <  *, 

b)  ddfinissant  la  position  angulaire  du  trlAdre  avion  par  rapport  au  triAdre  normal  terrestre  (fig.  2)  i 

41  azimut , 

(J  assiette  longitudinale , 

<I>  angle  de  glte  , 


c) 


4<2  (Ce  symbole  n’est  pas  normalisd)  angle  de  la  rotation  autour  de  l'axe  z  fixe,  amenant 
l'axe  x  d'une  position  x-^  A  une  position  Xg  dlffArente  de  x^. 

4^=  4*  si  et  seulement  si  z  est  vertioal. 

definiasant  la  position  dn  In  trujectoire  par  rapport  au  triedre  terreatre  *oyQz0  (fig.  3)  ! 


7a  pente  de  la  trajectoire,  par  convention 
la  azimut  aArodynaoiique , 


-5.  <  5 

2  *  2 


> 


pa  angle  de  glte  adrodynamlque . 


GRANDEURS  GEQMETRIQUES  DE  REEEHENCE 


S  surface  de  reference  ;  gAnAralement  pour  un  avion  la  surface  de  l'aile  prolongee  A  l'interieur 
du  fuselage  et  des  fuseaux, 

t  longueur  de  reference, 

b  envergure  de  l'aile. 


Note  1  :  de  nombreux  auteurs,  le  plus  souvent  de  langue  anglaise  adoptent  deux  longueurs  de  reference  : 
C  longueur  de  reference  longitudinale,  gdncralement  la  corde  aArodynamique  moyenne  de  l'aile, 

b  longueur  de  reference  transversale  dAja  symbolisAe  ci-dessus. 


De  plus,  chacune  de  ce s  longueurs  est  divisAe  par  2  dans  les  expressions  sans  dimension  des  vitesses 
angulaires  (voir  note  3). 


Cette  pratique  n'est  oas  con forme  aux  recommendations  de  In  norme  ISO. 
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GRANDEURS  CINEMATIQUES 
Vit«sses  lineal res  et  angulaires  i 

V  vi tease-air  (  vecteur  &  3  eoapoaantes  dont  la  grandeur  eat  le  scalaire  V  ) 

u  v  h  eoapoaantes  de  la  vitesae-air  sur  lea  axes  du  trifedre  avion, 
a  cdldritd  du  son, 

^  noabre  de  Hach 


Note  2.-  lea  coaposantea  de  la  vitesae-air  dans  le  trifedre  terreatre  et  le  tri&dre  a^rodynamique 
ont  lea  aeaes  systoles  affect^s  respectivement  des  indices  infdrieurs  "o”  et  "a"  . 


p  q  r 


P 

p* 

q* 

r* 


composantes  du  vecteur  M  aur  les  axes  du  tri&dre-avion,  respectivement  :  vitesse  de  roulis, 
vitesse  de  tangage,  vitesse  de  lacet 

derivee  reduite  par  rapport  au  temps  de  1 

ddrivde  reduite  par  rapport  au  temps  de  1 

acceleration  tangentielle  reduite  ( 
vitesse  reduite  de  roulis (  p»  _  Et  )  , 
vitesse  reduite  de  tangage  (  q  B  f 

vitesse  reduite  de  lacet  (  r,  a  rt.  j 


'angle  d' incidence  ( 
'angle  de  ddrapage  (  p 

V{  vv  dV  > 

VT>  ’  =  dt  }  ' 


a  = 


at 

*T, 

V  ’ 


da  i 
“  dt  1 
41  ) 
dt 


Note  3  :  ces  expressions  sans  dimension  sont  souvent  rapportees  a  d'autres  longueurs  de  reference,  deja 
signalees  dans  le  prec^dente  note  1 . 


&*  =  as  ,  9£ 

2V  *  2V 

=  is ,  P*=  & 

2V  *  p  2V 


vitesses  angulaires  longitudinales 


rb 

2V 


vitesses  angulaires  transversales. 


MASSE  ET  INERTIE  DE  L'AVION 


■  masse  de  l'avion 
Ix(ou  *)  moment  d'inertie  par 
X  (oil  B)  moment  d'inertie  par 

w  . 

I  (OU'.C)  moment  d'inertie  par 

i 


rapport  N 
rapport  a 
rapport  k 


l’axe  longitudinal  de  l'avion 
X'axe  transversal  de  l'avion 
l'axe  normal  de  l'avion. 


I  (ouD)  Droduits  d'inertie 

I„(ouE) 

Ixy(ou*‘) 


de  l'avion 


p  masse  rdduite  (  p  =  2m  /pS?) 


( 

( 

(  I 


I  » 

,yz 


*y 


Jyz  d*  )  , 
Jzx  dm  )  , 
fxy  dm  )  , 


FORCES,  MOMENTS  ET  COEFFICIENT  SANS  DIMENSION 

R  rdsultante  du  systeme  de  forces  s'exerqant  sur  l'avion,  incluant  les  forces  aerodynaraiques  et 

de  propulsion,  mais  excluant  les  forces  de  pesanteur,  d'inertie  et  les  forces  de  contact  avec 
le  sol, 

X  composantes  de  R  sur  les  axes  du  trifcdres  avion 

Y 

Z 

Y  y  z 

C  coefficients  sans  dimension  :  n  =  A  C,r  =  C_  *  .  _ 

c£  x  ^pV2S  Tf  ±pV2S  Z  *pv23 

ou  p  est  la  masse  voluraique  de  l’air  non  influence  par  le  champ  aerodynamique  de  l’avion. 


Corapcsanten  du  moment  resultant  du  systeme  R  sur  les  axes  du  triedre  avion,  usuclleme.it  par  rapport 
au  centre  de  gravity. 


L  moment  de  roulis 

M  moment  de  tangage 

N  moment  de  lacet 


coefficients  sans  dimension  : 


o* 


_* _ 

±pvzst 


Note  4  :  ces  coefficients  sont  souvent  exprimes  ovec  d'outres 

la  pr4c4dente  note  1 . 


|pv*se 


c 


n 


_N _  . 

ipV*St 


longueurs  de  reference,  dejN  signalres  dans 


C 


( 


L 

ip7*Sb 


Cm 


M 

ipV*Sc 


H 

ipV^Sb 


DERIVEES  PARTIELLES  DES  COEFFICIENTS  SANS  DIMENSION  DE3  COMPOSANTES  DE  LA  FORCE.  RESULTANTS  ET  DU 
MOMENT  RESULTANT. 
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Appelees  d^rivees  de  stability  ou  deriv4es  a^rodynamiques . 

0  Derivees  premieres, 

C^E  d^rivee  du  coefficient  par  rapport  a  la  variable  E 

R  designe  l*un  des  symbolcs  X,  Y,  Z,  {  ,  m,  n 

E  designe  l'un  des  symboles  tt  *,  0  »  p,  q  •  r  etc.  repr^sentant  une  des  variables  d'etat 

de  I'avion,  par  exemple 

r  iCu  dM/da  4Cm  4M/d^  _  ^  4M/6^ 

°ao  =  a  a  =  i pV2  Sf  *  ®q  4q  *pv2S?  ’  <*“  =  *$*  =  *pV2Sf 

2)  Derivees  secondes. 

C|»  eie2  d^rivee  seconde  du  coefficient  Cj,  par  rapport  aux  variables  E.  et  E?  i  E.  peut  §tre  egale  N  E?. 
par  exemple 


* C_ 

.  =  _ ■  _ 

42M/da2  d2C, 

d2L/4a4  f 

'ma2  4a2 

ipV2Sf  *  °laP  da  dp 

4*0*  42M/da  dq^ 

~TevTs T  ’ 

]  —  C  — 

mqa  maq 

dq*da  “  ■J-pV2St 

Relations  entre  lec  coefficients  des  derivees  de  stability  definis  :  uivant  les  norraes  ISO  et  ceux  utilises 
courammer.t  par  les  auteurs  americains. 


exemple 

de 

coefficient  Cjy 
oonforme  aux  notations 

rapport 
coeff.  a 

derivee 

a  :  IoO 

b  :  U.d 

cceff.  b 

derivt'e  do  position: 

longitudinale 

trar.Hversale 

4  m/ 4a 
di/dp 

ipV2secmao 
±pV2StCtp  p 

ipv2Sc  C^a 
ipv2  sb  aep  P 

e/t 

b/e 

d-'rivee  par  re.porrt 
a  la  vi tease 
angulaire  : 
longitudirale 
transverse le 

4M/dq* 
4L/ 4p* 

ipV2StCmq  qr/V 
*pV2SfCtp  pC/V 

ipv2ScCBq  qo/2V 
ipV2SbCfp  pb/2V 

2(c  ID2 
2(b/t)z 

apoelo  "  C  dynaniquc-"  dor.t  1' expression  eat  : 


pD 


Cnpc08a  ~  j  c^psina 


VuRUS  A3RC  0Y!«  AMIt&E  (  DU  FLANEUR) 

Coefficients  sar.r.  dimension,  facteur  de  charge. 

Rn  frreo  nr  rodyartmique  du  plar.eur  d/finie  arbitrairement  comru  la  difference  vectorielle  entre 

la  f  <rce  R  et  la  pour  sc' e  F. 

X*\  Y“  oor.pcsantes  de  KA  sur  les  axes  du  triedre  avion, 

X^f  Y‘\  .fl  con non antes  do  RA  sur  les  axes  du  triedre  aorodynamique,  auxquelles  sent  donnees  los 
d<Vo  ninntionr.  suivnr.ter.  : 

trafi.i'e  (eu  resistance*  n»'rouynamiquo)  compnnnnte  de  RA  sur  l’axe  x^,  change e  de  signe. 
force  Intervale,  composantc  deRa  sur  y^  non  changer  de  signe, 
nrjrtance,  eumpor.nnte  de  RA  sur  z^  chang^e  de  signe . 

> 

4 

V 
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Coefficients  sans  dimension  ^ 

A 

C  coefficient  de  trainee  ou  de  resistance  aErodynamique  Cv  »-  Cy. - ,  dEsignE  par  C„ 

1  1  II  y*  fl  D 

dans  les  publications  de  langue  an?1 lise.  *  ’  ® 

Cy  coefficient  de  face  latErale  Cy»Cy#»  Yjf S 

C  coefficient  de  portanceCe  V*  par  C,  dans  lea  publications  de  langue 

anglaise.  **  */  *  L 

Facteur  de  charge. 

~S  vecteur  facteur  de  charge,  c'est  le  rapport  de  la  force  rEsultante  au  module  du  poids  de 

1' avion  : 

n  dans  le  texte  est  la  oomposante  du  facteur  de  charge  relative  it  la  portanoe,  aoit  i 

ni-ni  ■  -  Z«/ mg 

g  acceleration  de  la  pesanteur.  *  ' 

2  -  NOTATIONS  NON  NOHMALISEES  j 

Coefficients  a4rodynamiques  des  Equations  des  petits  mouvements: 

Mouvements  longitudinaux. 

Xu  Xw  Zu  Zw  voir  les  definitions  chapitre  I  §4-5 

Mouvements • transversaux  (ou  lateraux)  : 

y  *  -  n«  >  f*  >  no  •  fp  >  nr  voir  les  definitions  chapitre  I§  4-6 

P  P  P  "  r 

aquation  diffErentielle  du  second  ordre  a  coefficients  constants  et  racines  complexes  : 

X  ♦  2  S ca,x*  u>*x  «0 

x(t)  Elongation  du  mouvement 

x  derivee  premiere  par  rapport  au  temps 

x  dErivEe  seconde  par  rapport  au  temps, 

fo»  respectivement  pulsation,  frequence  et  pEriode  propre 

»  f.  T  respectivement  pulsation,  frequence  et  ps eudo-pEriode 
£  coefficient  d'amortissement  rEduit, 


t. 


temps  de  passage  de  1* amplitude  de  sa  valeur  initiale  a  la  valeur  kx^  ! 

.  -t/T,  , 

constante  de  temps  (  X=Xe8  ). 


NOTATIONS  VECTORIELLES  -  MATRICIELLES  -  OPERA TEURS. 
Yecteurs  : 


I 

[J,] 

UJ 

Veotours  unitaires  : 


vecteur  de  composantes  dans  un  certain  rEfErentiel  (ou  une  certaine  base 

ou  un  certain  repere) . 

matrioe  oolonne  oorrespondante. 


matrice  transposEe ,  la  notation  ^J^jT  est  aussi  utilisEe. 


x  y  z 

relatifs  au 

repEre  Gxyz 

VA 

relatifs  au 

repere  Gx.y.z, 

Vo‘o 

relatifs  au 

repere  Gx,y,Z, 

h 

l'indioe  x 

dEslgne  le  repEre  utilise 

Gx  y  z 

la 

1' Indies  a 

dEsigne  le  repEre  utillsE 

Gxf.z. 

Jo 

l'indioe  o 

dEslgne  le  repEre  utilisE 

G  JqjfeZo 

II 


«  toraeur  des  forces  adrodynaaiques  et  de  propulsion  dans  Is  re pi re  C  x  y  z 
d' expression  transposde  : 

Cg  torseur  des  coefficients  sans  dimension  correspondents 
8  torseur  des  forces  d'inertie 

«x  vecteur  pesanteur  dans  r,  x  y  z 

6k  vecteur  d'dtat  (  ou  torseur  distributeur  des  vitesses  )  dans  G  x  y  z 
Ex  ddrivde  par  rapport  au  temps  de  Ex 
At  )  accrolssement  petit  du  premier  ordre 
()*  vecteur  dont  les  diluents  sont  sans  dimension. 

Matrices 

1  matrice  unite 

•  l_  matrices  unitd  ddgdndrdes 

•  ,  matrice  diagonals  d’dldments  unite  affectds  de  signes  let  - 

0  vecteur  ou  matrice  nulle 

■n  tenseur  de  masse 

In  tenseur  d'inertie 

Mll  Mr  IV|L 

CRL,EL  CBTET 

CrL.ET  C«TEL 

Opdrateurs 

D  operateur  de  derivation  (  matrice  6x6  ) 

d  operateur  de  derivation  (  matrice  3x3  ) 

A  prodult  vectoriel 


matrices  d'dldments  sans  dimension  reprdsentant  des  termes  d'inertie 
matrices  ddfinies  au 

matrices  dont  les  elements  sont  les  coefficients  sans  dimension  des 
derivdes  adrodynamiques  de  staMHtd 


D'autres  notations,  non  couramment  utilisdes,  sont  prdcisdes  dans  le  texte. 
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INTRODUCTION 


La  determination  des  dArivAes  aAro dynamiques  de  stability,  en  soufflerie  et  au  tunnel  de  tir,  sur  dee 
maquettes  animAes  de  petite  mouvements  oscillante,  ou  de  mouvements  continus  de  rotation  uniforme  de 
faible  vitesse  angulaire,  remonte  A  une  trentaine  d’annAes  environ. 

Theoriquement ,  les  rAsultats  de  ces  experiences  donnent  l'ensemble  des  dArivAes  des  oomposantes  du 
torseur  des  forces  aArodynamiques  par  rapport  &  celles  du  vecteur  d'etat  qui  dAf in it  le  regime  de  vol 
de  1 ' aeronef .  Cependant,  il  faut  rappeler  que  les  dArivAes  de  ce  torseur  prises  par  rapport  aux  varia¬ 
bles  autres  que  la  vitesse  angulaire  peuvent  @tre  obtenues  A  partir  des  mesures  sur  maquette  fixe. 

Le  professeur  VALENSI  est  l'auteur  d'un  des  premiers  documents,  publiA  en  1954,  faisant  le  point  des 
methodes  et  moyens  utilises  pour  oes  essais  dans  les  souffleries  europAennes. 

Puis,  en  1955,  paratt  un  premier  AGARDograph  d' ARNOLD  Lee  sur  les  mesures  dynamiques  en  soufflerie. 

En  1959,  ORLIK-H0CKEMANN  publie  une  revue  des  methodes  de  mesure  des  dArivAes  de  stability  dynamique  des 
aAronefs. 

Quelques  annAes  aprAs,  en  1967,  le  m6me  sujet  est  traite  dans  un  second  AGARDograph  de  SCHUELER,  WARD  et 
HODAPP. 

Les  references  de  ces  publications  sont  classees  par  ordre  chronologique  [  7  A  10  ]. 

Au  cours  de  ces  dix  deraieres  annAes ,  un  accroissement  notable  des  besoins  d' informations  sur  la  stabi¬ 
lite  s' est  manifeste  par  suite  du  dAveloppement  des  techniques  aArospatiales  dans  les  di verses  directions 
caracterisees  par  : 

-  1' extension  du  domaine  de  vol  des  avions  de  transport  en  transsonique  et  en  supersonique  ; 

-  la  mise  en  service  d' avions  A  decollage  et  atterrissage  court  ou  vertical  (ADAC-ADAV),  en  anglais 
STOL-VTOL  ; 

-  l'extension  du  domaine  de  vol  des  avions  de  combat  aux  evolutions  A  facteur  de  charge  Aleve  en 
transsonique  j 

-  les  recherches  sur  la  rentree  balistique  puis  sur  la  navette  spatiale  ; 

-  les  recherches  de  configurations  d'avion  A  commande  automatique  gAnAralisAe  (CAG),  en  anglais  CCV 

[  16  ]. 

La  complexity  des  questions  posees  par  les  qualites  de  vol  de  ces  aAronefs ,  la  necessity  de  rAduire  les 
coflts  et  dAlais  de  leur  mise  au  point  imposaient  d' entreprendre  1 ' etude  approfondie  de  leur  stability  dAs 
le  dAbut  du  projet,  d ' ou  la  necessity  de  recourir  systematiquement  aux  essais  sur  maquette  dans  les 
installations  au  sol. 

Ces  essais  ont  AtA  effectuys,  pendant  cette  pAriode,  dans  les  installations  dycrites  dans  le  prAcAdent 
AGARDograph  [ 1 0  ].  La  plupart  d'entre  elles  sont  encore  en  service  A  l'heure  actuelle. 

Des  descriptions  plus  dAtaillAes  de  ces  installations  et  d' autres  mises  en  service  depuis  ont  yty  prAsen- 
tdes  dans  diverses  publications  d'ORLDC-RtJCKEMANN  [11A15],  A  la  suite  d'une  enquSte  sur  les  besoins  d' in¬ 
formations  sur  la  stability  aux  grands  angles  d'incidence  et  sur  les  moyens  d'essai  permettant  de  les 
satisfaire. 

Plus  rAcemment,  les  comptes  rendus  des  communications  aux  deux  derniAres  rAunions  de  spAcialistes  de 
l'AGARD,  sur  le  DECROCHAGE  ET  LA  VRILLE  DES  AVIONS  MILITAIRES  (PMP  1977)  et  sur  les  PARAMETEES  DE  STABI¬ 
LITE  (FDP  1978)  foumissent  les  AlAments  necessaires  A  une  vue  d' ensemble  sur  ces  questions.  II  apparatt 
A  la  lecture  de  ces  documents  que  des  progres  considyrables  ont  AtA  accomplis  en  mAtrologie  grfice  A  la 
mise  en  service,  au  cours  de  la  demiAre  dAcennie  dans  les  centres  d'expyrimentation  en  vol  et  au  sol. 
de  moyens  modernes  d' acquisition  et  de  traitement  numArique  des  informations  d'essai.  On  observe  que  ces 
moyens  ont  At A  mis  intensAment  A  contribution  dans  1'expArimentation  en  vol  sur  la  stability  par  les 
mAthodes  "d* identification  de  paramAtres".  Mais  leur  utilisation  dans  les  installations  au  sol  n'avait 
AtA  prAvue  que  pour  obtenir  en  temps  rAel  les  coefficients  aArodynamiques  sans  dimension  stationnaires 
classiques.  Toutefois  leur  disponibilitA  a  permis,  dans  un  premier  temps,  d'allAger  l'acquisition  et  le 
traitement  des  mesures  des  coefficients- de  stability  quasi-stationnaires  par  les  methodes  lineaires 
d' oscillations  de  faible  amplitude,  de  sorte  que  les  mesures  entrent  maintenant  dans  le  cadre  des  essais 
oourants  sur  maquette  en  soufflerie. 

Ces  mAthodes  ont  AtA  rapidement  Atendues  A  de  plus  larges  domaines  de  variation  de  paramAtres  :  par 
exenple  les  intervalles  d'incidence  et  de  dArapage  des  balances  stationnaires  ont  AtA  accrus,  d'autre 
part,  des  souffleries  ont  AtA  AquipAes  de  balances  animAes  d'lin  mouvement  de  rotation  uniforme  n 
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II  est  ainsi  devenu  possible  de  determiner  les  composantes  3tationnaires  du  torseur  d’ effort  aerodynami- 
que  en  fonction  de  quatre  param&tres  d’etat  «  »  p  ,  et  Cl  (dans  le  cas  ou  ce  vecteur  est  paral- 
Iblo  A  l’ecoulement).  Ces  composantes  ne  varient  lineairement  que  dans  un  petit  domaine  de  valeurs  de  ces 
paramAtres.  Dans  le  reste  du  domaine  explore  (grandes  incidences,  ecoulement  transsonique)  les  composantes 
ne  sont  plus  lindaires,  elles  peuvent  aussi,  ainsi  que  leurs  ddrivees,  presenter  des  di3continuites. 

II  arrive  m@me  d* observer  des  phenomAnes  d' hysteresis. 

Or,  les  deriveea  determinees  par  les  methodes  lindaires  dans  les  zones  de  discontinuites  sont  sujettes  & 
caution  ;  et  l’on  constate  qu’il  n'existe  pas  &  l'heure  actuelle  de  moyens  au  sol  traitant  de  faqon 
satisfaisante  1* ensemble  des  cas  de  vol  instationnaires  dans  lesquels  les  forces  aerodynamiques,  A  un 
instant  donne,  dependent  non  seulement  des  valeurs  instantandes  des  paramAtres  d’etat,  mais  aussi  des 
valeurs  prises  par  ces  paramAtres  pendant  une  oertaine  periode  prdcddtuit  1' instant  consider.  C’est  le 
cas  du  vol  A  haute  incidence  lorsque  1' ecoulement  sur  l’adronef  est  le  siAge  de  decollements  et  de 
recollements,  ou  encore  des  ecoulements  transsoniques  avec  ondes  de  choc  instationnaires. 

Les  forces  exercees  en  vol  sur  l'aeronef  au  cours  d’une  evolution  determinee  de  ce  genre  ne  peuvent  dtre 
simu^ees  avec  certitude  sur  maquette  qu’en  imposant  a  celle-ci  une  evolution  angulaire  semblable  et  non 
un_jnouvement  elementaire  a  _priori.  Les  techniques  de  vol  libre  semblent  rdpondre  a  cette  necessite,  toute- 
fois  raalgre  les  perfect  ionnanents  apportes  ces  dernieres  annees  aux  installations  de  vol  libre,  coniine  par 
exemple, ^celles  de  la  NAbA  a  LANGLEY  FIELD  [i7]  et  celles  de  l'Institut  de  Mecanique  des  Fluides  A  Lille 
L19J,  1' ecoulement  reste  encore  du  domaine  de  1 ’incompressible  et  l’influence  du  Reynolds  y  est  negligee. 

L' experience  en  soufflerie  sur  une  maquette  disposant  d’un  nombre  de  degr4s  de  liberte  suffisants  et 
couplee  avec  un  ordinateur  effectuant  le  calcul  de  la  trajectoire  peut  permettre  une  simulation  plus  fiddle, 
sous  reserve  bien  entendu  de  respecter  le  nombre  de  Mach  et  un  nombre  de  Reynolds  suffisamment  (Sieve.  Dans 
un.cas  comme  dans  1* autre,  1* analyse  des  resultats  exige  la  connaissance  et  l'emploi  des  techniques  de 
mecanique  du  vol. 


II  a  done  paru  intAressant,  suivant  les  conseils  de  M.  i'Ingenieur  General  CARRIERE,  qui  a  bien  voulu 
susciter  et  orienter  ce  travail,  de  presenter  dans  cet  AGARDograph  destine  a  une  synthese  des  methodes 
et  moyens  d’essais  au  sol  un  expose  prAliminaire  des  techniques  modernes  de  calcul  de  la  dynamique  du 
vol. 

* 


* 


* 


Ie  sujet  est  traite  dans  deux  chapitres  precedes  d'une  liste  de  notations  et  symboles.  Des  complements 
theoriques  et  pratiques  sont  indiques  en  annexe. 

Deux  ensembles  de  notations  sont  utilises  : 

-  le  premier  coraprend  les  notations  conf«>rmes  aux  recommandations  ISO  sur  les  "Terraes  et  Symboles  de  la 
Mecanique  du  Vol",  officiellement  adoptees  par  1* ensemble  des  nations  europeennes.  0  A  4J 

La  generalisation  de  leur  emploi  devrait  faciliter  la  lecture  des  documents  et  des  publications  sur  la 
dynamique  du  vol.  Elle  permettrait  aussi  d'eviter  des  erreurs  d' interpretation  dues  a  la  diversity  des 
notations  et  des  grandeurs  de  reference,  choisies  de  faqon  plus  ou  raoins  arbitraire,  principalement 
lorsqu’il  s'agit  d'engins. 

Quant  aux  notations  relatives  aux  avions,  elles  sont  habituellement  conformes  a  une  symbolisation  dejA 
ancienne.  La  symbolisation  ISO,  plus  complete,  parait  mieux  adaptee  aux  besoins  actuels. 

La  correspondance  entre  ces  deux  symbolisations  est  precisee  dans  la  liste. 

-  Le  second  contient  les  notations  matricielles ,  non  encore  normalisees ,  utilisees  dans  l'espace  de  confi¬ 
guration  des  torseurs  A  six  dimensions  :  torseurs  d' effort,  torseur  distributeur  de  vitesse  (ou  vecteur 
d’etat)  et  torseur  d ’acceleration. 

Le  premier  chap*  tre  est  un  rappel  de  notions  elementaires  de  dynamique  du  vol,  rAdige  plus  partieuliere- 
ment  a  1* attention  des  experimentateurs  de  soufflerie,  generaleraent  peu  informes  sur  ce  sujet. 

Les  composantes  du  vecteur  d'etat  qui  caracterise  A  un  instant  donne  les  conditions  de  vol  sont  prAcisAes 
et  les  definitions  des  "deriveea  de  stability"  ou  "derivees  aerodynamiques"  sont  rappelees  au  debut  de 
ce  chapitro. 

Le3  considerations  theoriques  qui  suivent  sont  limitAes  A  1 ’ Atablissement  des  equations  des  petits  mouve- 
inents  dans  les  domaine3  lineaire  et  non  lineaire. 

Dans  le  domaine  lineaire  cla33ique,  les  deux  systfemes  d’ equations  lineaires,  a  coefficients  constants, 
des  petits  mouvemenfcs  longitudinaux  et  transversnux ,  decouples  entre  eux,  sont  etablis  en  partant  des 
equations  gAnerales. 

Puis,  aprAs  quelques  transformations,  ils  sont  prAsentAs  sous  la  forme  cunonique  utilisee  habituellement 
en  dynamique  du  vol.  Cette  analyse  montre  comment  les  coefficients  sans  dimension,  qui  expriment  les 
derivees  aerodynamiques,  interviennent  pratiquement  dans  les  calculs  de  dynamique  du  vol. 

Dans  le  domaine  non  lineaire,  deux  exemples  sont  prAsente3,  comportant  chacun  une  non- linearite  parti- 
culiAre,  observee  dans  des  etudes  de  missiles. 


14 


Dans  chacun  des  cas  traitds  ci— dessus,  la  recherche  des  solutions  n'est  pas  abordde. 

Les  expressions  des  Equations  gdndrales  de  la  dynamique  du  corps  solide  sont  dgalement  rappeldes  en 
annexe,  k  titre  d' aide-memoire.  Biles  sont  ddcrites  sous  3  formes  :  vectorielle,  matricielle  et  ddve- 
loppde. 

Les  notations  matricielles  utilisdes  pour  presenter  ces  notions  fondamentales  permettent  d'ecrire  sous 
forme  condensee  et  sans  confusion  possible  des  expressions  mathematiques  dont  les  developpements  sont 
lourds  et  penibles  k  lire.  Elies  facilitent  ainsi  l'analyse  des  phenomknes  k  etudier  et  servent  avanta- 
geusement  k  etablir  les  programmes  de  traitement  sur  ordinateur,  des  informations  acquises  en  vol  et  en 
soufflerie. 

Un  langage  commun  aux  disciplines  qui  regissent  la  dynamique  du  vol  devrait  necessairement  renforcer  la 
collaboration  entre  les  specialistes  mecaniciens  et  adrodynamiciena, rendue  indispensable  par  la  complexity 
des  questions  qui  se  posent  k  l'heure  actuelle. 


Les  techniques  d' experimentation  qui  sont  exposdes  succinctement  dans  le  second  chapitre  repo3ent 
sur  la  remarque  suivante  :  dans  le  domaine  de  valeurs  du  vecteur  d'dtat  dans  lequel  le  torseur  des  forces 
adrodynamiques  est  continflment  ddrivable,  l'dquation  des  petits  mouvements  de  stability  peut  6tre  prdsen- 
tde  comme  un  systkme  d'yquations  lindaires  dont  les  ddrivdes  adrodynamiques  sont  les  inconnues. 

II  en  rdsulte  que  les  ddrivdes  adrodynamiques  peuvent  6tre  ddtermindes  en  traitant  par  des  mdthodes 
lineaires  des  informations  obtenues  en  soufflerie,  soit  sur  maquette  fixe,  soit  sur  maquette  en  rotation 
uniforme  parallels  a  l'dcoulement,  soit  encore  sur  maquette  animde  de  mouvements  k  nombre  restreint  de 
degrds  de  libertd. 

Des  descriptions  ddtailldes  des  montages  d'essais  dtant  ddjk  rassemblds  dans  les  documents  rdfdrencds  ci- 
dessus,  il  a  paru  suffisant  de  ne  rappeler  ici  que  les  principes  des  montages,  surtout  pour  signaler  les 
p r ogres  rdalisds  dans  les  sys terr.es  de  mise  en  mouvement  des  maquettes.  Quelques  schdmas  d 1  installations 
typiques  figurent  cependant  en  annexe.  Apres  avoir  prdcise  la  nature  des  informations  d'essais,  les 
principes  des  systkmes  d'acquisition  et  de  traitement  analogiques  et  numdriques  sont  dnoncds  en  montrant 
les  progres  apportes  par  les  moyens  numeriques.  Des  indications  sont  donndes,  k  titre  d'exemple,  sur  une 
mdthode  de  traitement  non  lindaire  appliqude  k  des  essais  en  soufflerie  pour  lesquels  den  cycles  limites 
ont  ete  observes  [2  0  ]. 

II  est  fait  dgalement  mention  des  perfectionnements  apportes  aux  moyens  d'essais  ausol  sur  maquette  en 
vol  libre. 


Un  travail  de  synthkse  sur  un  sujet  aussi  vaste  et  complexe  comporte  necessairement  des  omissions  dont 
l'auteur  prie  le  lecteur  de  bien  vouloir  l'excuser. 


IS 


CHAPITR*  I 

elements  de  dynamique  DU  VOL 


I.I.-  GENERA  LLTES  . 

Le  present  chapitre  contient  las  notions  thAoriques  essentielles  qui  sarvent  de  base  A  l’Atude  de  la 
stability  des  aAronefs. 

Cette  Atude  est  menAe  couraament  de  faqon  approchAe  en  admettant  les  hypotheses  suivantes: 

La  variation  de  la  masse  volumlque  de  l'air  dans  le  domains  de  vol  considArA  est  nAgligeable, 
la  nasse  de  1'aAronef,  les  AlAments  du  torseur  des  forces  de  propulsion  dans  un  repAre  liA  A 
l'aeronef  sont  constants  et  1'aAronef  est  indAformable. 

Ces  hypotheses  rApondent,  A  de  rares  exceptions  prfes,  aux  conditions  expArimentales  usuelles  dans  les 
installations  au  sol. 

L'action  des  gouvernes  n'est  pas  AtudiAe. 

Les  definitions  des  dArivAes  de  stabllltA  sont  rappelAes  ainsi  que  celles  des  coefficients  sans  diaension 
qui  les  reprAsentent  habituellementj  ensuite  11 'eat  montrA  comment  elles  sont  utilisAee  pratiquement  dans 
les  cas  les  plus  courants. 

Des  publications  rAcentee  de  dynamique  du  vol,  de  H.H.2.M.  THOMAS  [21][22)et  de  J.  C.  VAHMER  ( 23 }  pouraont 
etre  consul  tees  pour  des  informations  plus  approfondies  et  dAtalUAea  sur  oe  sujet. 

Ia  dynamique  du  vol  d'un  aAronef  est  un  cas  particulier  de  la  dynamique  du  corps  Bolide  dont  les  Aquations 
gAnArales  sont  rappelAes  en  annexe  A. 1.1. 

Les  composantes  du  vecteur  d'Atat  E-=  [u  v  M  P  1  r3  et  celles  du  torseur  des  forces 
sont  rapport Aes  A  un  triAdre  G  x  y  z  lie  A  1'aAronef  et  ayant  son  origins  au  centrd  de  gravitA  G. 

Dans  le  cas  de  1'aAronef  le  torseur  des  forces  est  la  somme  de  trois  torseurs  ! 

-  des  forces  d'inertie, 

-  de  pesanteur,  Gg 

-  des  forces  aArodynamiques  et  de  propulsion  ft. 

L'Atude  de  la  stabllltA  d'un  rAglme  de  vol  de  1'aAronef  (  croiaifere,  montAe,  descents,  arrondi,  virage . . . ) 
est  conduits  de  facon  classic ue  en  analysant  son  comportement  lorsque  qu'll  est  lAgArement  AcartA  de  sa 
position  de  rAgime. 

Un  Acart  donnA  de  cette  position  est  caractArisA  par  un  accroissement  du  vecteur E*  auquel  correspondent  des 
accroissements  des  trois  torseurs  v  ,  GgSt  It  .  Pour  Aviter  des  confusions  ces  accrolssements  sont  dABimAs 
par  la  notation  A( )  ,qui  reprAsente  une  auantltA  petite  du  premier  ordre. 

Cruelaues  indications  sont  donnAes  sur  les  rAglmes  de  vol  avant  d*  aborder  l'Atude  thAorique  de  la  stabllltA. 

I.  2 .  REGIMES  DE  VOL.-  Ils  peuvent  etre  dAfinis  et  classAs  par  ordre  de  dlfficultA  crolssante  en  ce  qui  con¬ 
cerns  l'Atude  de  leur  stabllltA.  Cette  Atude  est  classique,  depuis  de  nombreuses  annAes,  dans  le  cas  des 
rAglmes  HnAalres.  mais  elle  n’a  pu  etre  menAe  A  bien  A  l'heure  actuelle,  que  dans  quelques  cas  partlculiers 
de  regimes  non-llnAalres . 

1.2.1-  RAglmes  HnAalres. »  Le  torseur  des  forces  aArodynamiques  et  de  propulsion  de  oes  rAglmes  es;  une 
fonctlon  llnAalre  du  veoteur  d'Atat.  soit  : 

"  =  . 

dans  cette  relation  ■<x  est  une  matrice  dont  tous  les  termes  sont  des  constantes  par  hypothAse. 


Ces  £gns$antes  80nt.  par  dAfinition,des  jfT\vies  ,Ae  stabllltA  _lii-Jirealar_arArfl. . 


Parmi  les  rlgimes  rApondant  A  ces  conditions  se  trouvent  i 

-  les  rAglmes  de  vol  normal  des  avions  ayant  un  plan  de  symAtrie  (croisiAre,  montAe,  descents,  arrondi). 
Le  vecteur  d’Atat  et  le  torseur  aArodynamique  de  ces  rAglmes  ont  leuxs  composantes  transversales  nulles. 

La  trajectoire  du  centre  de  gravitA  eat  rectlligne,  ou  de  faible  courburs,  et  de  faible  pente.  La  vitesse  ast 
constants  .  Ces  rAglmes  sont  appelAs  dans  ce  qui  Buit,  rAglmes  symAtriques  HnAalres. 

-  les  rAglmes  de  manoeuvre  normals  des  avions  rendus  lAgArement  dissymAtriques  par  les  braquages  des  gu- 
vernes  ou  le  rAglage  des  forces  de  propulsion  (virage,  glissade  ). 

Les  vecteur* Exetfl  de  ces  rAglmes  ont  des  composantes  longitudinales  et  transversales ^non-nulls*!  mals, 
gAnAralement  les  composantes  longitudinales  de  R  ne  dApendent  que  de*  composantes  de  meme  nature  de  Cx ,  il 
en  eat  de  meme  de  leure  composantes  transversales. 


Lorsque  cette  condition  set  vdrifide,  le  regime  de  vol  peut  etre  oonsiddrd  comma  un  rdglne  dlaeymdtrlcue 
sana  couplaga  adrodynaalaue. 

1.2.2-  Regimes  non-llndalres. 

Ce  sont  des  regimes  pour  lesquels  il  n'existe  pax  de  relation  lindaire  entre  R  et  Ex  •  II  eat  tou jours 
possible  d'dcrire  formellement  R-  t^Ex  >  ®ai8  dans  ce  cas  les  dldments  de  sont  des  fcnctions  de 

Ex  et/ou  du  temps. 

Les  rdgimes  suivants  rdpondent  gdndralement  &  ces  conditions  i 

-les  rdgimes  de  vol  des  engine  de  forme  dlancde  dont  les  termea  de  rappel  et  d'amortissement  adrodyna- 
miques  n'dvoluent  pas  lindairenent  avec  l'angle  d'incidence. 

-  les  manoeuvres  A  facte ur  de  charge  dlevd  des  avlons  de  combat  en  transsonique . 

-  les  figures  acrobatiques, 

-  les  cas  de  vol  marginaux  dangereux  pour  la  sdcuritd  de  1* avion  (  ddcrochage,  vrille.etc.  ). 

Dans  ces  divers  rdgimes,  les  accelerations  lindaires  et  angulaires,  les  angles  d'incidence  et  de  ddrapage 
atteignent  des  valeurs  importantes,  pour  lesquelles  les  dldments  de  Kx>  c'est-k-dire  les  ddrivdes  de  sta- 
bilitd,  dvoluent  de  faqon  importante  en  fonction  de  Ex . 

II  arrive  meme  que  les  variations  de  R  en  fonction  de  Ex» 

-  no  sont  pas  des  fonctions  ddfinies,  mais  des  fonctionnelles  qui  ddpendent  des  valeurs  de  ff ■„  fonction 

du  temps,  * 

-  prdsentent  pour  certalnes  valeurs  de  Eg  des  dlscontinuitds  pour  lesquelles  des  dldments  de  K,  ne 
sont  pas  ddfinis, 

-  sont  multiformes,  solt  it  la  suite  du  passage  par  une  dlscontlnultd,  solt  par  effet  d'hystdrdsls. 


1.2.3-  Expressions  du  vecteur  d’dtat,  du  torseur  adrodynamlaue  et  de  leurs  accrolssements  du  premier  ordre. 
Variables  adrodvnanlques. 

-  Le  vecteur  d  'dtatE*  est  compoedde  deux  vecteurq,  la  vitesse  du  centre  de  gravitdV  et  la  rotation  £2  . 
Son  expression  dans  les  axes  avion  peut  etre  dcrite  sous  diverses  formes  transposes  : 

EX=VUQ  3  [  u  v  w  p  q  r]  =  [Vcosdcosp  Vsinp  Vslnacos0  p  q  r)  . 

Les  variables  couramment  utilisdes  dans  les  dtudes  d'adrodynamique  V,  d,0  apparaissent  dans 
1' expression  dcrite  k  droite. 

Le  torseur  des  forces  adrodynamlques  et  de  propulsion  a  pour  expression  transposde  dans  lee  mSmes 

axes:  *-lx  I  Z  L  M  V  ]  . 

-  L'accroissement  du  vecteur  d'dtat  peut  etre  prdsentd  sous  deux  formes  transposes  : 

AEX“  [du  Av  A»  Ap  Aq  Ar]  et  AE*-  [AV/V  dp  Aa  Ap*  Aq*  Ar*]  , 

qui  sont  lides  par  la  relation: 

AE*=I*  AE„/v, 

L' expression  ae*  'cans  dimension  est  utilisde  oouranment  dans  les  mesuree  des  ddrivdes 

de  stabilitd.  L'exprer  sion  ddveloppde  de  f?  est  dcrite  en  annexe  A. 1.4. 


■L'accroissement  du  torseur  adrodynamique  a  pour  expression  transposde  : 

AR  =  [dX  AT  AZ  AL  AH  AH] . 

S'il  existe  une  relation  lindaire  entre  R  et  E  , leurs  accrolssements  sont  relids  par  : 

AR  =  K  AE  avec  K  =  KxlY’V  . 

lee  expressions  prdnddswtmm  ne  sont  valables  que  pour  des  valeurs  deA(  )  petites  du  premier  ordre. 
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1.3-  STARILITS  13' UN  REGIME  OB  VOL. 

L'Atude  thAorique  de  la  stability  d'un  rAgime  de  vol  constate  it  rechercher  led  solutions  des  Aquations 
diffArentielles  des  petits  mouvements  de  1'aAronef  qui  donnent  Involution,  en  fonction  du  temps,  de 
l'Acart  du  vecteur  d'Atat  A  sa  valeur  de  rAgime. 

a)  -  Si  le  regime  est  llnAalre.  les  Aquations  diffArentielles  sont  linAaires  &  coefficients  constants. 

Dans  la  premiere  parti e  de  ct  paragraphs,  l'Atablissement  de  ces  Aquations. est  rappelA  dans  les  caa  les 
plus  classiques.  Elies  sont  prAsentAes  sous  la  forae  canonique  couramment  utilisAe  dans  la  pratique. 

Les  relations  entre  les  expressions  des  dArivAes  de  stabilitA  figurant  dans  les  Aquations  et  les  coefficients 
aArolynamiques  sans  dinension  qui  les  reprAsentent  sont  prAcisAes. 

b)  -  Les  moyens  d'Atude  de  la  stabilitA  des  rAgimes  non-llnAaires  paraissent  avoir  AtA  trfes  peu  dAveloppAs 
A.  l'heure  actuelle  pour  diverses  raisons.  En  ce  qui  conceme  les  avlons.  la  complexitA  des  phAnombnes 
aArodynamiques  se  prete  difficllement  au  choix  d'un  modAle  mathAmatique  les  reprAsentant  correctenent  et 
les  moyens  de  calcul  nAcessaires  aux  applications  numAriques  d'approximations  thAoriques  viennent  A  peine 
d'entrer  ans  la  pratique  courante  .  De  plus,  lesnoyens  d'essais  susceptibles  de  foumir  les  valeurs  des 
coefficients  aArodynamiques  dans  ces  rAgimes  sont  encore  k  crAer. 

II  a  done  fallu  se  contenter  d'Atudier  la  stabilitA  avec  les  moyens  expArimentaux  existants  ,  par  exempje 
avec  des  mesures  en  souff lerie  A  grande  incidence  sur  maquette  fixe,  en  oscillations  de  faible  amplitude  ct 
en  soufflerie  verticale  de  vrille,  sans  porter  d'effort  thAorique  particulier. 

Par  contre,  enee  qui  conceme  les  missiles  1'expArience  a  montrA,  dAs  les  annAes  1950,  que  les  mAthodes 
linAaires  Ataient  trop  insuffisantes  et  qu'il  Atait  indispensable  de  faire  appel,  dans  la  mesure  du  possible 
aux  mAthodes  non-linAaires  existantes. 

Ces  mAthodes,  basAes  sur  les  travaux  de  P0INCARE(  i88l)  et  de  LJAPOUNOV  ( 1893)  ont  AtA  dAveloppAes  en  vue  de 
rApondre  aux  besoins  crAAs  par  les  techniques  d1 automatisms;  elles  ne  peuvent  Avidemment  etre  Atendues  A 
l'ensemble  des  systAmes  non-linAaires,  mais  elles  permettent  de  traiter  de  nombreux  cas  particuliers  ren- 
contrAs  dans  la  pratique. 

Des  traitAs  d’initiation  A  ces  mAthodes  ont  AtA  publiAsi  parmi  lesquels  celui  de  W. CUNNINGHAM  ,  en  langue 
anglaise,  paru  en  1958 1 24  ) .  En  langue  franqaise,  la  nouvelle  Adition  du  traitA  de  CH.  GILUC  ,  P.  DECAUIME 
et  M.  FEIECRIN,  publiAe  en  1975.  donne  un  aperqu  complet  de  l'At&t  actuel  de  la  question  et  contient  une  im- 
portante  aocumentation  1251. 

Quelques  remarques  sur  des  applications  de  mAthodes  non-linAaires  aux  Atudes  de  dynamique  du  vol 
et  un  exemple  simple  d'Atablissement  des  Aquations  de  stabilitA  d'un  rAgime  de  vol  non-linAaire  d'un 
engin  de  rAvolution  AlancA  sont  prAsentAs  au  paragraphe 

1.4  -  STABIUTE  D'UN  REGIME  DE  VOL  L1NEAIRE  SYMETRIQUE. 

Dans  ce  rAgime  le  plan  de  symAtrie  de  1 'avion  reste  vertical  et  fixe  par  rapport  A  la  terre.  Le  cas  le  plus 
courant  est  celui  du  vol  A  vl tease  constants  oar  une  trajectoire  rectiligne  de  faible  pente,  qui  eet  le  Beul 
analysA  ici. 

Ce  rAgime  est  caractArisA  par  les  valeurs  suivantes; 

Vconstant  ,  n= B  ,  u  constant  ,  p  =  0  ,  constant, 

le  vecteur  d'Atat  E^et  le  tomeur  R  ont  pour  expressions  transpoaAes  : 

ix=[u  0  w  0  0  0]  - )  1=[V  0  o  0  0  o] , 

et 

ft  -  [X  0  Z  0  0  0]  . 

lea  expressions  des  accroissementsAEet  AReont  Acrites  au  paragraphs  I.2Jci  deaeus. 

U  valeur  de  la  pente  de  la  trajectoire  ou  ceil*  de  l'assiette  longitudinals  •  llAes  par, 

e  =  T^+  a  , 

complete  les  caxactArlstiquee  du  rAgime.  Ia  valeur  de® intervient  dans  les  projections  de  la  peeanteur  sur 
les  axes  du  tribdxe  avion}  elle  a  done  une  influence  sur  les  conditions  de  stabilitA  en  vol  et  en  soufflerie. 

En  position  de  rAgime,  1'accAlAration  est  nulls,  1' Aquation  du  aouvement  exprime  que  le  torseur  dee  forces 
extArieures  eat  nul. 


R  +  G*  =  0  , 

GreprAaente  la  peeanteur,  son  moment  jar  rapport  au  centre  de  gravitA  est  nul. 

-  dans  le  trlAdre  normal  terreetre, 

[0  0  mg  0  0  0]. 

-  dans  le  triAdre  avion  et  dans  une  attitude  quelconque  0/  0,  0 

=  mg  [-slid)  slrifcoa®  cosBeoe®  0  0  0]  < 
en  rAgime  de  vol  symAtrlqua ,  ♦  -  0  , 

Gx  ®  mg  [  -ein®  0  co#®  0  0  0  3  . 
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1.4. I  -  Equation  des  torsgure  d*  effort  des  petite  mouvemente  da  stability. 


L*s  torseurs  d' efforts  resultant  d'un  petit  dcart  de  l'adronef  de  sa  position  de  revise  sont  lies  par: 

(i)  A  <7=  AR  ♦  A6X  . 

I*  torseur  des  forces  d'inertie  peut  s'dcrire  sous  la  forme: 

^  =  D.m.Ex  , 

dont  l'expression  ddveloppde  est  donnde  en  annexe AI.  1  t  son  accroissement  a  pour  expression: 

A 9  *  A(D.m.  Ex)  =  Ao.m.  Ex  +  D.m.AEX| 

En  regime  de  vol  symdtrique  k  tra jectoire  rectiligne ,  D  est  au  second  ordre  pr&s  une  matrice  diagonal* 
i'opdrateurs  d/dt  d'orire  6,  qui  peut  etre  reprdsentde  par  : 


D  *  () . 


Le  symbole  ( 'f  utilise  dans  la  suite  ddsigne  l'opdrateui  sans  dimension  :  d/dt  r()  . 


le  terms  : 


A££i sAD.ro. Ex  , 


iont  le  ddveloppement  est  lndiqud  en  annexe  A1.5,  est  le  torseur  des  forces  d'inertie  d'entrafnement  dues 
i  la  rotation. 

im  et  {)  etant  nymdtriques  sont  commutables,  de  sorte  que: 

D.m  .  AEX-  ro.A  Ex  , 


1* equation  (l)  se  presente  alors  sous  la  forme: 

a 

(2)  m.AEx  s  An  *  AGK  -  A9n 

avecf 

AR  =  2  VlSkACR  et  AGX  =  mg  j[-cose  0  -sine  0  0  O]A0+[O  cose  0  0  0  0]A4>{. 

ACr  est  1* accroissement  du  torseur  des  coefficients  adrodynamiques  sans  dimension  et  It  une 
eatrice  diagonale  d' ordre  6  dont  les  £ laments  sont  dgaux  A  1  et  -1  pour  les  forces,  etf  pour  les  moments, 
conformdment  aux  ’-ecommandations  des  normes,  rsoit: 


k  = 


-1  . - 

N.’X  0 

I  %  »\ 

!  0  \f 


Dans  le  cadre  des  hypotheses  dcrites  au  ddbut  de  ce  chapitre,  le  vecteur  d'etat  est  ndcessairement  une 
fonction  continue  du  temps,  et  dans  la  pratique  on  admet  que  ACpest  une  fonction  lineal re  de  AE*et  de  sa 
ddrivde  premiire  par  rapport  au  temps,  Afesoit: 

ACr5[  CREt.  C^(f]AE*. 

Ties  ddrivdes  adrodyraunlques  de  stability  sont  les  dldments  des  matrices  carries  d' ordre  6 

qui  reprdsentent  les  ddrivdes  partielles  du  torseur  sans  dimensionCfjpar  rapport  aux  dldments  deEet  E  . 
Elies  sont  constantes  en  rdgime  lindaire, 

XI  en  rdsulte  l'expression  suivante  de  AR  : 

AR  =  yV2Sk  CREA£*+  jV'SkCj^iAE* 

1.4.  I  -  Tenseur3  de  masse  entrafnde  et  de  masse  apparente. 

[a  presence  d'un  terme  enAEdans  l'expression  de  ar  conduit  k  introduire  la  notion  de  masse  entrainde ■ 
la  ddrivde  AE„qui  intervient  en  dynamique  du  vol  s' exprime  en  effet  par  : 

AE„  =  ^V’AE*, 

de  sort?  que  1'dquation  du  mouvement  nrend  la  forme  : 


(3)  (»nem.)AE|  r|vlSlcCR(EAE*+ AGj-Ae^j , 

lens  laquelle  on  a  po*4i 


me  s-ffikCRgl5  y  avec  p  =2^  (  maB3e  rdduite  )■ 
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Compte  term  des  dimensions  de  k  et  de  R  (annexe  A. 1. 4)  les  elements  de  ITI^ relatifs  aux  viteases 
de  translation  ont  les  dimensions  d'une  masse  d'ou  l'appellation  traditionnelle  de  masse  en trainee  j 
les  elements  relatifs  aux  vitesses  angulaires  ont  les  dimensions  d'une  inertie.  La  matrice  *n0  peut  done 
Stre  appeiee  tenseur  de  ma3se  entratnee  et  le  tenseur  m+fn*  represente  la  masee  apparente  de  I'a^ronef. 

1.4.3^  ToiBeur  d 'accelerations  relatives. 

Par  ailleurs  (voir  annexe  A. I. 5)  est  une  fonction  lindaire  de  An  ,  done  de  Aft*  «  c©  term© 

peut  done  se  combiner  avec  af  pour  aboutir  It  la  presentation  de  (3)  sous  une  forme  compacte  : 

(A)  AEX*  MAE*  4-  Aflx 

Les  expressions  de  la  matrice  M  et  du  vecteur  Aflx  3 ont  l£velopp4©8  aux  paragraph© s  suivants  dans 
lesquels  les  Equations  des  petits  mouvements  de  stability  longitudinal©  et  de  stability  transversal©  sont 
d tab lies  sdpardment. 

Remarques  :  Les  matrices  Cfl.E  et  Cr  £  comportent,  en  principe,  chacune  36  elements,  mais  les  valeurs 
d'un  grand  nombre  d'entre  eux  sont  nuAes  par  raison  de  symdtrie  ;  dans  le  cas  de3  regimes  de  vol  oourants 
des  avions  et  des  missiles,  la  masse  fluide  entralnde  est  trbs  faible  au  regard  de  la  masse  de  l’a6ronef, 
le.  plupart  des  termes  de  CM  peuvent  Stre  negligSs.  Par  contre.  si  la  masse  volumique  de  1'aSronef 
est  voisine  de  celle  du  flui<Je  dans  lequel  ils  Svoluent  [26  ]  et  [  27  ],  il  faudra  tenir  compte  d'un 
plus  grand  nombre  de  termes  de  Cr(£  ,  e'est  le  oas  particulier  des  dirigeables. 


-  Le  to: 
de  Mach 
sible. 


torseur  sans  dimension  ACr  et  les  matrices  de  ddriv4es  Cop  et  Cr.  j  sont  fonction3  du  nombre 
ch,  sauf  aux  regimes  de  faible  vitesse  subsonique  ou  le  fluiae  peut  Stre  considers  comme  incompres- 


L' expression  de  ACr  en  regime  lineaire,  Scrite  ci-dessus,  doit  Stre  complSt4e  par  l'addition  du  torseur 
ACr(^T)  ,  d'  express  ion  : 

AC„(^)=AC„  ^  sJ?  =(Cr,e~£  +^E*+  Cr,^  +^)A^r 

avec 

i*=[v/v„  00000],  l*=[vf/v0l  00000] 

ou  V,  est  une  vitesse  de  reference  dont  la  valeur  est  fixSe  arbitrairement . 

Les  termes  entre  parentheses  restent  constants  par  hypothSse  au  cours  des  petits  mouvements  de  stability 
autour  de  E*  et  E*  donnes  ;  en  consequence  la  stability  d'un  regime  de  vol  fonction  du  nombre  de  Mach 
dans  ces  conditions  est  encore  rSgie  par  une  equation  linSaire  de  la  forme  (3)  ou  (4). 


1 . 4  A-  Bxnresalona  des  comnosan 
T  l'aSronef. 


tudinales  et  transversales  des  toraeurs  dans  le  renfere  114 


T'K  0  I 


avec  e.= 


En  decomposant  la  matrice  "unite"  en  une  somme  de  deux  matrices  degenerees  d’ordre  6  et  de  rang  3  s 

'1  v-— . ;  [o'-: . "»* 

f'.o\  0  I  0  [ 

1 =  ®L+®T  avec  ®L=  I  '  \b\  i  et  •t!=  I  \'K  i 

j  0  \i\‘  !  0  \cN 

_i . -voj  - ~m. 

Le  produit  A  gauche  d'un  torseur  (ou  d'unvecteur)  quelconque  de  rang  6  par  et  par  •-!-  permet  de 

disjoindre  les  composantes  longi tudinales  indicSes  L,  des  composantes  transversales  indiodes  T  de  ce 
torseur. 

A  titre  d'exemple,  le  calcul  est  applique  au  torseur  des  forces  aerodynamiques  et  de  propulsion  R  dont 
les  composantes  longi  tudinales  Rj_  et  les  composantes  transversales  R.y  ont  comme  transposees  : 


•|F s  FL  =[X  0  Z  0  MO] 
•^R=  Rj  *  [  0  Y  0  LON] 


et  l'on  retrouve  R  par  la  somme  : 

R  =  RL+  Rt  =  [  X  Y  Z  L  M  N] 

De  m£me,  les  transposes  longi  tudinales  et  transvarsales  de  AE*s  'ecrivent  : 

1  «lAE*=  AlE*  =  [Av*  0  A«  0  Aq*  0  ] 
(  •tAE*=AtE*=[  0  Afi  0  Ap*  0  A  r*J 
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at  l'on  retrouve  ae  an  faisant  la  aomme  : 

AE#=  AlE#+AtE*=  [AV‘  Ap  Ac  Ap*  Aq#  Ar*]  . 

Lorsque  le  torseur  est  exprime  par  la  produit  d'une  matrioe  par  un  vecteur,  par  exemple  A«MAE* 
da  l'dquation  (4),  la  disjonction  peut  Atre  effectude  an  ondrant  uniauement  aur  la  matrice.  ai  l'on  tient 
compte  de  la  relation  : 


A=(«l  +  *t)M(*i+*t)AE#, 

qui  peut  s' exprimer  par  : 

A  =  (•lM  «l+,  «lM  9j  *  »T  IWI  eT  +  exM«i.)AE* , 

ce  qui  revient  &  dcrire  M  30US  la  forme  : 

M  *  +  MTL  t  +  IW»T  • 


Compte  tenu  de  la  presence  des  colonnee  de  zdro  dans  lee  quatre  matrices  composantes  de  M  ,  lee 
expressions  des  composantes  longitudinales  et  transversales  de  A  s'ecrivent  : 


Al  =  M  L  A  lE*+  m[AlE* 

Af  ■  My  AjE  +  My  AxE* 


Des  exemples  de  M  [  et  Mj  sont  donnes  aux  paragrapr.es  suivants. 

Des  exemples  de  matrices  croisdes  sont  donnas  en  annexes. 

Dans  le  cas  du  vol  symdtrique  Mj  =  =  0  ,  et  liquation  (4)  peut  6tre  ecrite  sous  la  fonme 

de  deux  Equations  k  variables  inddpendantes  : 


(4.1) 

(4.2) 


A|EX  =  m£alE*+  A,.*,  , 

Axtx  -  MxAxE*+  Ax9x  . 


Leo  equations  4.1  et  4.2  sont  les  Equations  des  petits  mou-'ewents  de  stabilite  longitudinale  et  de  stabi- 
lite  transversale.  Les  torseurs  d1 accelerations  relatives  dans  le  repkre  lid  k  l’adronef  figurent  aux 
premiers  membres  des  dquations. 

1.4.5-  Stability  longitudinale. 

Les  equations  des  petits  mouvements  de  stability  longitudinale  sont  gdndralement  rapportdes  aux  axes  de 
"stabilite"  lids  a  l'adronef  et  confondues  en  position  de  rdgime  avec  les  axes  adrodynamiques  Gx.y^z,  . 

L'intdrSt  de  cette  disposition  est  de  simplifier  les  expressions  de  la  matrice  ML  et  du  vecteur 
A^g  ,  qui  peuvent  ktre  ensuite  rapportdes  a  des  axes  X  Z  diffdrents  par  un  changement  de  repkre. 

Dans  le  repkre  de  "stabilite",  1' angle  d'incidence  est  nul  et  l'angle  d'assiette  longitudinale  est  dgal  a 
la  pente  de  la  trajectoire  : 

«=0  ,  0  =  Y, 

Les  expressions  du  vecteur  d'dtat  Ex  et  E  sont  identiques,  ainsi  que  celles  de  leurs  accroissements  : 

Ex=  E  =  [V  0  0  0  0  0]  ,  Aix=  AE  =  [aV  0  VAC  0  Aq  0]  . 

L'dquation  (4.1)  peut  alors  s'dcrire  sous  la  forme  : 

(  ;  )  ALE=  M^aALE*+ ALgw 

avec:  Nlfca=  f  V*8  [(m+nwfrlM  CBlEL+  et  ALg,a=  [( m  ♦  ALGI8  . 


L' indice  a  qui  designe  les  axes  de  "stability",  n'est  pas  r£p4t£  de  fagon  sys tematique  aux  paragraphes 
a)  et  b)  ci-dessous  pour  soulager  l'^criture. 
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"fc.  se  prdsente  sous  la  forme  d'un  produit  de  matrices  par  un  scalaire,  a  savoir  : 

a)  ®RL  EL  ddrivdes  des  coefficients  a^rodynaraiques  longitudinaux  par  rapport  aux  elements 
*  longitudinaux  de 

*fc  coefficient  de  lq*  de  la  composante  longitudinale  de  1’ acceleration  d’entrainement 

de  G  ,  dirig^e  suivant  l’axe  z  .  dont  le  calcul  est  indique  en  annexe  A. I. 5. 

Le  produit  de  la  somme  de  ces  deux  matrices  par  k  ,  ddfini  ci-dessus  A  la  suite  de  liquation  (2)  a 
pour  expression  ddveloppee  : 

r  -2C«  0  -Cyn  0  -Cv«  0  1 


k(CRL,EL+*^t^aF 


0  0  0 

•2Cx  0  -Cxoe, 
0  0  0 


0  .Czq  +1*  0 

0  0  0 


0  0  ?Cmtt  o  fCmq 


[(m+mJ^T  matrice  inverse  de  la  masse  apparente  longitudinale, 


m  tenseur  de  masse  defini  en  annexe  A  1.1.2. 

m«L  tenseur  de  masse  entrainee  longitudinale  d’ express ion  : 

«n«L  =  “fs/kC^iL^a-SkCRL^i.^ 

et  sous  forme  developpee  : 


*tCRL,EL*«  = 


0 

0 

0 

a 

0 

0 

0 

0 

0 

•c  ti 

0 

-Czi 

0 

-Kii, 

0 

0 

0 

0 

0 

0 

C 

tcm* 

0 

fC»na 

0 

^rnq 

0 

0 

0 

0 

0 

0 

0 

Pratiquement  1‘ordre  de  grandeur  de  la  masse  reduite  H  ,  dont  la  valeur  se  situe  dans  l'intervalle 
50  <  |i  <  500  dans  les  exemples  de  la  reference  [23  ]  permet  de  n^gliger  -£xq__  et  J?zq__  dans  CRL,EL 
et  la  plupart  des  coefficients  de  Cpipj  a  l1 exception  de  C^q 

II  est  done  possible  d'obtenir  une  bonne  approximation  de  la  masse  apparente  longitudinale  avec  1* expres¬ 
sion  suivante  : 


m  - . - 

r-?’\  • 

s 

I  •  \9 


r-i  v . ir  i 

!*0\  0  |  *  S 

_m  i  \-1  >.  !  . . t, 

m  }  ;  o  o  o  : 

i  0  \  j  •  *  ^ 
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et  son  inverse  : 


[(m+me.)^7,=  l 


1 

, 

!  V 


0 


:  I  mf  r  K  m 

i  ift! ' 

L . . a  0 


Dans  l'expression  de  ecrite  oi-dessous,  le  terme  A _ C^g  (ligne  5,  colorme  2)  est  multiplie 


(My 


Dar  u  seul  &L4ment  non  nul  de  • 

II  faut  toutefois  signaler  que  les  valeurs  de  Cio  .  «, 
(voir  annexe  A. 1.2.4)  ;  il  est  done  utile  de  les  connaitre. 

Expressions  de  M^_  et  m£ 

Compte  tenu  de  la  relation  : 

£V2g  -El’ 

2  V  5  ~  m 

L'expression  approchde  de  M^J  a  pour  forme  developpee  : 


et  Cjq  entrent  dans  le  caloul  de  Cma  Cmq 


2C* 

0 

CX£L 

0 

0 

0 

2CZ 

0 

C2et 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 1» 

0 


f^-ma  0  ^5(Cmq+Cma)  0 
0  0  0 


avec  i  1  _ 

ly 

On  remarque  que  les  dldments  de  la  matrice  ddveloppee  au  second  membre  sont  sans  dimensions. 

Les  coefficients  C*  et  C j  sont  ceux  de  la  polaire  Eiffel  du  planeur  symbolisds  par  -Cj,  et  -Cj, 
au  paragraphe  1.7.1  des  normes  ISO  et  franqaises  ;  ces  coefficients  sont  representds  par  les  symboles 
C0  et  q  dans  les  publications  en  langue  anglaise. 

En  regime  de  vol  symdtrique  (  p  =  0),  les  deux  reperes  Gx.y.Z,  et  Gxyz  ont  en  commun  l'axeljy 
Une  rotation  a  dans  le  sens  positif  amAne  ,XaZa  sur  XZ  .  II  en  resulte  que  les  dldments  de 
rapportds  A  Gxyz  sont  lids  A  ceux  de  IVlL*par  la  relation  suivante  (voir  annexes  A. 1. 2. 2.  et 
A.I.2.3)  : 


ML=P«MLaPa 


qui  peut  §tre  repr^sentee  sous  forme  d£velopp£e  par  l'expression  : 


Xw  0  0  0 

0  0  0  0 
Zw  0  0  0 

0  0  0  0 
mw  0  mq  0 
0  0  0  0 


Les  auteurs  des  r£f£rences  6i]  6?  et  C23]  utilisent  les  m@mes  symboles,  mais  avec  des  definitions  diff4- 
rentes  h  un  facteur  scalaire  prfes. 
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Dana  cette  expression  lea  elements  de  la  matriee  d6veloppee  au  second  membre  sont  sans  dimension.  Ils 


s'expriment  en  fonction  des  Elements  de  C 
par  la  relation  : 


RL,EL  et  CftL,EL 


xu 

2CX 

^z'Cxa 

Cz* 

0 

Xw 

Cx# 

2Ck-ClB 

-2&z 

0 

Zu 

2Cx-CXa 

-Cxa 

0 

Zw 

^za 

2Cj+  Cxa 

20* 

0 

mw 

0 

0 

0 

Cma 

mq 

0 

0 

0 

Cmq  + 

,  retenus  apres  les  approximations 


cos2a 

sinacosa 

*m3<r 


Le  vecteur  AL9>a  a  pour  expression  transpose  : 

ALg,a=  gf-cos y,  0  -siny,  0  0  O]i0 

en  rappelant  qu*en  position  de  regime,  dans  le  repkre  , 

0  =  Ya 

Par  ailleurs  (annexe  A.I.3-1)  10  est  relie  a  Aq  par  la  relation  cinema tique  : 

AO  -  lq  =  0 

Et  dans  le  repere  Gx  y  z  : 

*  *a  ALBx# 


d' expression  transposee  : 


1.4.6-  Stability  transversale . 


ALgx  =  g  [-cos(  y,«-a)  0  sin(ya*a)  0  0  0]  A©  . 


Les  equations  des  petits  mouvements  de  stabilite  transversale  ne  sont  pas  n^cessairement  rapportees  aux 
axes  de  "stabilite"  mais  les  aeronefs  ayant  generalement  un  plan  de  symetrie,  les  axes  a  et  z  sont 
choisis  dans  ce  plan. 

Les  parametres  du  regime  de  vol  se  presentent  alors  de  la  fapon  suivante  : 

a  et  V  constants,  ©sy.+  a  constants,  fi-fi- 0  . 
d'ou  les  expressions  transposes  du  vecteur  d'etat  : 

Exs  [u  0  w  0  0  0]=  [Vcosa  0  Vsino  0  0  Oj 

et  de  son  accroissement  : 

A^Ex  =  [0  vAp  0  Ap  0  Ar]  A7i;s^i*=[0  Ap  0  AP*0  Ar*] 

L' equation  des  petits  mouvements  transversaux  s'£crit  sous  la  forme  : 

(6)  Aji=  A-,€*+ATgx 


=  f  V’s  [(m+me^rWCRT  E^T, 

Les  matrices  qui  forment  la  matriee  mt  ont  pour  definitions  : 

a)  CflTfET  matriee  des  derivees  des  coefficients  transversaux  par  rapport  aux  elements  de  Af®* 

&  coefficients  de  Ap*  et  Ar*de  la  composante  transversale  d'acceleration  d'entralne- 

ment  de  G  ,  dirigde  ouivant  l'axe  y  dont  le  ca.lcul  est  indique  en  annexe  A  1.5. 


24 


Le  produit  de  la  somme  de  ces  deux  matrices  par  k  a  pour  expression 


b)  [(m+me)T] 


TV 


0 

0 

0 

0 

0 

0 

0 

cyP 

0 

Cyp+  psinu 

0 

V 

ET+"^T>  = 

0 

0 

0 

0 

0 

0 

0 

Kq 

0 

«*p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

fCnp 

0 

ftnr 

matrice  inverse 

de  la 

masse 

anparente  transversale. 

m  tenseur  de  masse, 

tt)a  tenseur  de  masse  entralnee  transversale  d' expression 


et  sous  forme  d^veloppde  : 


»ne=-?ltCRXixPI, 


_,T  _  m 
m»T  =  j 


0 

0 

0 

0 

0 

0 

0 

H 

0 

cyp 

0 

Cyf 

0 

0 

0 

0 

0 

0 

0 

lcri 

0 

*c*p 

0 

% 

0 

0 

0 

0 

0 

0 

0 

fCnp 

0 

^np 

0 

*C  „f 

cosa  0 
0  1 
-sina  0 


stna 

0 

cosa 


N 


Pratiquement.  l'ordre  de  p:andeur  de  p  permet  de  negliger  Cyp  et  Cyr  dans  CRX  ff  et  la  plupart 
des  terraes  de  CRX>£X  a  1' exception  de  Cfp  et,  Cnj  •  II  est  done  possible  d'ottenir  une  bonne 
approximation  de  la  masse  apparente  transversale  avec  1' expression  suivante  : 


(m  +  me^= 


0  0  0  I 

0  m  0  '  0 

0 

* 

0  0  0  [ 

m 

. 

. 

j  lx  0  -l« 

0  0 

0  1 0  0  o 

0  0  0 

0 

l-Ixz  0  Iz 

0  *Cn$  0  i 

. 

'  0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

i 

o 

» 

£ . 

0 

:  ix 

I  m 
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Et  la  matrice  inverse  a  pour  expression  d£velopp4e 


[(m  +  *  1 
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Dans  l'expression  de  Mj  Norite  ci-dessous,  les  dl&nents  en  1/H  de  la  matrice  pr^c^dente  sont  multi¬ 
plies  par  les  elements  de  proportionnels  a  |l 


Expression  de  Arly 
Compte  tenu  de  la  relation 


VS 

rv  b  u  7 


l'expression  developpee  de  Mj  peut  s  * ecrire  : 


0  0 
0  y» 


M  5 


Les  elements  de  la  matrice  developpee  au  second  membre  de  cette  expression  sont  sans  dimension.  Ils 
o'expriment  en  fonction  des  Elements  de  Chi’  g-j  ^RTET  retenus  apres  les  approximations  par  : 

yp=cyp,  yp  =  »*"«  ,  yr  =  -cosa, 

V  -  i^Clp  +  rfV*  fP=  ^ein«+j55(c„pdcn.sina)],  2r=  ^^cosa+^C^-C^cosa)] , 

nP  ~  7|j(Cnp  +  T~Ct^  *  "p  =  7fJCnp+  Cnp8in  “  Ctp4C|p8lna)]  •  =  r^nr"^08  0 + ]”(Cfr  -C^c°«u)]  . 

r*2  =  et  r£J=lz/m/} 

Le  vecteur  At«x  a  pour  expression  transpos^e  : 

AT*x=  [o  cose  0  0  0  0]  gA<t>  , 

\®  e3t  reli6  a  S?  et  Sr  par  la  relation  cin£matique  (Annexe  A  1.3)  : 

\4>  -  Ap-  Artge  =  0  . 

Si  le  calcul  avait  ete  conduit  dans  les  axes  de  "stability"  (  a=  0  ),  les  termes  en  a  auraient 

dicparu  et  les  coefficients  de  ma3se  entrainee  ne  figureraient  que  dans  les  elements  fr  et  Of  Or, 

cette  representation  ne  paralt  pas  int^ressante  car  elle  neces3ite  de  reprendre  pour  chaque  regime  4tudi£ 
le  calcul  du  tens* ur  d* inertie  en  fonction  de  a  .11  semble  plus  simple  de  passer  de  Mk  A  Hflr 

et  de  con3erver  Mt  sous  la  forme  pr^sentdc  dans  ce  paragraphe. 

L' etude  qui  vient  d'etre  expo3($e  raontre  comment  sont  ^tablies  les  Equations  lineaire3  utilia^es  couramment 
en  dynamique  du  vol  en  precisant  les  simplifications  adoptees  dans  la  pratique.  Ce3  Equations  expriment 
le  torseur  des  acceleration  relatives  en  fonction  des  coefficients  a^rodynamiques  de  stability  3ans 
dimension.  En  particulier  la  presence  de3  termer  Cpnq+Cpna  et  Cfp  +  Cri  % intt  ,  Cnr-CnpcOta 
est  raise  en  Evidence.  **  etc.. 
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La  presentation  des  equations  de  stabilite  longitudinale  et  de  stabilitd  transversale  avec  des  matrices 
d'ordre  six  peut  sembler  fastidieuse  dans  le  cas  d'un  regime  de  vol  symdtrique  k  trajectoire  rectiligne, 
m tuu  son  intdrSt  apparatt  immediatement  dks.  qu'il  s'agit  d'une  manoeuvre^au  cours  de  laquelle  inter— 
viennent  simultandment  des  elements  de  M£  ,  Ml  .  MT  et  My 

Quelques  indications  sur  1' etude  de  la  stabilite  de  ces  regimes  particuliers  sont  presentees  aux  para- 
graphes  suivants. 

I  .5  -  CAS  DE  VOL  NON-LINEAIRES  A  STABILITE  LINEAIRE 
Les  torseurs  d'effort  resultant  d'une  manoeuvre  quelconque  de  l'aeronef  verifient  liquation  : 

(7)  «7  =  |v2SkCR+  Gx  , 

Dans  la  plupart  des  cas  cette  equation  est  non  lineaire,  mais  si  les  angles  d* incidence  et  de  derapage 

restent  moderes,  la  non  linearity  n'est  pas  de  nature  a£rodynamique,  elle  est  due  au  couplage  par  inert ie. 

Les  coefficients  aerodynamiques  sont  pratiqueznent  les  m@mes  qu'en  regime  lineaire  syiodtrique  £tudi£  ci- 
dessus. 

L' equation  de  stabilite  ecrite  sous  la  forme  (4)  est  encore  lineaire  puisque  les  6ldments  de  la  matrice 
M  sont  constants,  mais  la  manoeuvre  est  dissymdtrique 

*  p  0  at  /  ou  My  *  0  . 

Des  elements  de  nature  aerodynamique  dus  a  la  dissymetrie  de  1* ecoulement  peuvent  aussi  figurer  dans 
Mp  et  .  L' equation  du  mouvement  (4)  est  lineaire  si  ces  Elements  sont  constants. 

Par  ailleurs  il  est  possible  d'etendre  k  certains  cas  de  vol  1' etude  approchee  de  la  stabilitd  en  deux 

groupes  independants  soit  en  analysant  separement  la  stabilitd  longitudinale  et  la  stabilitd  transversale 
(vol  rectiligne  avec  un  derapage  de  quelques  degres),  soit  en  decouplant  1' etude  des  mouvements  de  rota¬ 
tion  de  celle  des  mouvements  de  translation  (vrille).  Dans  ce  dernier  cas  on  peut  admettre  que  la 
direction  de  V  reste  constante  au  cours  d'un  petit  mouvement  de  rotation  (les  relations  entre  p  ,  q  , 

r  et  a  ,  £  repondant  k  cette  condition  sont  indiquees  en  annexe  A  1.3). 

Les  deux  groupes  d' equations  lindaires  siraplifides  ainsi  determinees  sont  gdneralement  discutds  sous 
forme  analytique,  alors  que  le  cas  gdneral  ne  peut  pratiquement  ?tre  aborde  qu'a  l'occasion  d'applica- 
tions  numeriques  sur  ordinateur. 

Le  calcul  approche  des  elements  de  M  expose  ci-dessus,  prdsente  l'interSt  de  mettre  en  dvidence  les 
termes  des  matrices  m  ,  ma  ,  CR  £  et  dont.  l'influence  sur  la  stabilitd  est  prdponddrante. 

Mais,  dans  les  applications  il  est  desormais  beaucoup  plus  aise  d'effectuer  directement  le  calcul  complet 
sur  ordinateur. 

Si  la  matrice  M  est  connue,  les  solutions  de  l'dquation  lineaire  (4)  sans  second  membre  donnent  les 
amortissements  et  les  pdriodes,  constantes  caracteristioues  du  mouvement. 

1.6-  REMARQUES  SUR  LES  DC”  A  INKS  DE  VOL  A  STABILITE  HON  LINEAIRE. 

Tout  domaine  de  vol  pour  lequel  les  valeurs  de  M  sont  variables  dans  le  temps  et/ ou  ne  sont  pas 

ddfinies  est  non  lindaire.  Quelques  remarques  sur  les  cas  de  vol  rdpondant  k  ces  conditions  sont 
exposdes  cl-dessous. 

l.ft.VCas  de  vol  oil  le  torseur  d33  forces  adrodvnamlques  e3t  contlnument  ddrlvable  par  rapport 
au  vecteur  d'dtat. 


Deux  exemples  sont  prdsentds  : 

1.6. 1.1-  o  et  V  variables,  Cp  f  constants. 

L'dtude  de  la  dynamique  de  ce  cas  de  vol  non-lindaire  n'est  pas  du  domaine  de  l’adrodynanicien, 
sous  rdserve  que  lee  coefficients  puissent  Stre  considdrds  conune  inddpendantades  nombres  de  Mach 
et  de  Reynolds. 

L'exemple  du  vol  ascendant  d'une  fusde  sonde  dtudid  par  M.  Bismut[28),  qui  vdrifie  cette  condition, 
n'est  citd  ici  que  pour  mdmoire.  la  figure  ci-dessous  montre  les  variations  non-lindaires  des 
ddrivdes  dimensionnelles  de  stabilitd  dues  k  la  variation  de  la  masse  volumiaue  de  l'air  p  avec 
l'altitude,  dont  11  faut  tenir  compte  dtant  donnd  l'onlre  de  grandeur  de  dp/dt. 


rappel  adrodjmamlque 
effet  Magnus 


All 
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|.  6,1.2  -  p  et  V  variables .  CR.E  et  CR.E  fonctlons  continues  et  ddrlvables  de  E  et  E  . 

Les  S16mnts  de  M  sont  alors  complbtement  ddfinls  dans  le  doma^ne  de  vnleurs  de  E*  et  E*  >  si,  de  plus, 
11  est  suppose  que  1'adronef  ne  rencontre  pas  de  rafale  atmosphdrique  se  traduisant  par  des  discontinuitds 
sur  a  et  (l  ,  E*  et  £*  sont  des  fonctlons  continues  du  temps  . 

II  existe  alors,  k  tout  instant  t  ,  une  solution  lindaire  tangente  k  la  solution  exacte. 

Les  pdriodes  et  amortissements  de  ces  solutions  lineaires  varient  de  faqon  continue  en  fonction  du 
temps,  Dans  certains  cas  particuliers  les  solutions  de  1* Equation  du  mouvement  rendent  compte  de  eye les— 
llmltes  dontla  frequence  ,  l'amplitude  et  la  stability  dependent  essentiellement  des  valeurs  des  coeffi¬ 
cients  adrodynamiques . 

L' etude,  faite  par  Murphy  et  Bradley  [29J  du  mouvement  "conique”  1'un  missile  dlancd  observe  sur 
maquette  au  tunnel  de  tir,  en  donne  un  example. 

La  mise  en  equation  de  ce  mouvement,  exposde  ci  dessous,  est  limitee  au  cas  du  missile  de  revolution 
pour  res ter  dans  le  cadre  dldmentaire  du  present  expose. 

Les  forces  adrodynamiques,  la  force  de  propulsion,  la  pesanteur  et  les  forces  d'inertie  en  acceleration 
lindaire  sont  supposdes  constamment  dquilibrdes  suivant  la  direction  de  la  trajectoire.  Le  nodule  de  la 
vitesse  peut  done  9tre  quelconque  et  les  equations  du  mouvement  exprimentseulement  l'equilibre  dec 
moments  par  rapport  au  centre  de  gravi 1 4 .  De  plus  si  1'  axe  x  est  pris  comme  axe  de  revolution,  li 
tenaeur  d'inertie  de  l'engin  a  pour  expression  : 


In - 


!x  0  0 
0  1  0 
0  0  I 


et  etant  donne  la  forme  dlancde  *x<  1  • 


De  m@me,  le  coefficient  du  moment  adrodynamique,  autour  de  G  x  ,  C.  est  suppose  nul,  de  sorte  que  le 
coefficient  du  torseur  des  forces  adrodynamiques  a  pour  transpoade  :  “ 

[0  Cm  Cn  ]  . — Cp  =  [  Cm  Cn  ]  , 

Ainsi  les  equations  du  mouvement,  au  nombre  de  deux,  expriment  l'dquilibre  des  moments  autour  de  deux 
axes  rectangulaires  Gy  ,  Gz  choisis  arbitrairement  dans  le  plan  normal  &  l'axe  de  revolution  Gx  ,  soit 

(8)  dm. 12  =— VzS£Cr  , 

(  0)  peut  £tre  remplacde,  compte  tenu  de  la  forme  diagonale  de  In  ,  par  1' expression  sans  dimension 
suivante  : 

(  9  )  d* .  12*  =  H  C  p  avec  H  -  S  ^  / 1  , 

Les  variables  d'dtat  inddpendantes  dtant  au  nombre  de  deux,  on  peut  choisir  par  exemple: 
l'angle  ostJfl/X)  de  la  trajectoire  avec  l'axe  de  revolution  ,  et  la  composante  Pg  de  la  vitesse 
angulaire  de  rotation  sur  le  vecteur  vitesse  de  sorte  que  le  vecteur  12  a  pour  transposde  dans  les 
axes  lids  : 


pjeosa  a*  p’ji  sinaj 

L' expression  ddveloppde  de  ( S )  s'dcrit,  en  faisant  figurer  les  ddrivdes  des  coefficients  de  stabilite 
adrodynamique  : 


ra 

dT* 

-  pj  cosa 

r  •  * 

a* 

=  H 

Cma 

^mq  + 

0 
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a* 

p*  coma 

1* 

p*  sina 

0 

0 

cnr 

p*sina 

(10) 

La  symdtrie  de  revolution  entrains  p  =  0  et  il  n'y  a  aucun  terme  en  Cnp  et  Cnp  da"8  les  Equations 

MURPHY  et  BRADLEY  supposent  que  l'engin  est  dynamiquement  instable  dans  un  secteur  de  o  symdtrique  par 
rapport  &  l'origine.  Cette  instabilitd  s'explique  si  l'on  admet  les  relations  quadratiques  : 

(») 


‘-mq  +  ^mot  =  ^nr+  Cnrc^  1  +  *  )  a 


r 


par  raison  de  sym4trie 
et  l1 on  a  pos£ 


>  Cnr+  ^nra2  a1 

Cmq+Cma  *Cnr 
Cmqa*+  Cmaa1  =0**"*)  Cnra* 


L' instabilitd  A  l'origine  s’exprioe  par  l'indgalitd  C„r  >0,  qui  correspond  A  des  coefficients  d'amor- 
tissement  de  tangage  et  de  lacet  positifs  a  l'origine. 
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D' autre  part,  en  supposant  C„ra*<0  ,  aes  coefficients  s'annulent  pour  des  valeurs  de  a  c&lcuiees 
d 'apres  (11).  Ils  devlennent  successivement  ndgatifs  pour  des  valeurs  de  a  crolssantes.  D'autre  part 
Cm(t  est  suppose  n^gatlf  et  constant  en  fonctlon  de  o  . 


Liquation  du  mouvement  s 'ecrlt  alors  : 

T  m2 


(12) 


■^1^2  *■  cnral  (l  +  •)]  ^  ^  ma  cos  a  sina 

2P2eo#a  •"»«(«►*  °«r  +  <W«J  ) 


=0 


Cette  Equation  n’est  pas  lineal re  ;  elle  comprend  des  termes  en  uJ  ,  p*  ,  sina  et  COSO  .  Elle 
adraet  deux  solutions  particuliferes  :  * 

1°  p*=  0  qui  v^rlfie  la  seconde  Equation  quelle  que  solt  la  valeur  de  u  ;  dans  ce  cas  le  mouvement 

de  l'engln  est  une  oscillation  plane  r^gie  par  une  Equation  de  Van  der  Pol  : 

“*-H[cSr  +CnraK1  +  a)  0  ]«*-HCma  a  =  0  , 

dont  la  solution  classique  conduit  k  un  cycle  Unite  pour  une  valeur  convenable  de  1’ amplitude. 

2°  p*  =  constante  et  u  a  la  valeur  qui  annule  C^r  »  c'est-i-dlre  : 

Cnr  =  0  =>  a  =-  <%r/ 

Cette  solution  indique  que  le  mouvement  "conique"  est  possible.  On  vd rifle  qu'il  est  stable  si  |>1  j 
cette  condition  sigralee  par  les  auteurs  se  retrouve  aisdment  k  partir  de  1' Equation  des  petlts 
mouvements  de  stability  qui  dans  ce  cas  particulier  est  lin^aire. 

I.  6. 2 -Cas  de  vol  ou  le  torseur  des  forces  aerodynamlques  est  derivable  par  morceaux,  statlonnalre  ou  non. 

Certains  cas  de  vol  k  incidence  elevee ,  comportent  des  discontinultes  aerodynamlques  provenant  de  d£col- 
lements,  le  plus  souvent  instationnaires,  qui  provoquent  des  perturbations  des  quality's  de  vol  de 
l'aeronef . 

II  n'exlste  pas  encore  de  methode  theorique  d'analyse  de  ces  discontinuity  et  il  est  n^cessaire  de  falre 
appel  a  l'experience  pour  tenter  d'y  remedier,  voir  une  synthese  recente  sur  ce  sujet  de  W.J.Mc  CROSKEY[30]. 

Quelques  exemples  de  discontinultes  de  cette  nature  d^terminees  par  des  mesures  en  soufflerie  3ont 
exposes  ci-dessous . 

1.6. 2.1  -  Dans  le  domalne  de  valeurs  de  E  ou  des  decollements  et  recollements  se  manifested,  les 
relations  entre  R  et  E  sont  generalement  rnultiformes,  plusieurs  valeurs  de  R  correspondent  alors 
ci  une  valeur  donn£e  de  E 

S'il  existe  un  cycle  llmite,  l'un  au  moins  des  elements  de  E  varle  de  faqon  cycllque  et  les  variations 
de  R  en  fonction  de  E  presented  souvent  une  (ou  des)  boucle  (s)  d 'hysteresis,  meme  si  la  variation 
de  E  dans  le  temps  est  tree  lente  ;  c'est  le  cas  par  exemple  des  mesures  statlonnalres  effectuees  par 
valeurs  discretes  en  rotation  unlforme. 

Lorsque  E  ^  O  ,  mais  varie  de  faqon  quelconque,  la  forme  de  la  boucle  depend  generalement  de  l' evolu¬ 
tion  de  E  dans  le  temps. 

Ia  figure  1.2,  tlree  de  la  reference  [31]  en  donne  un  exemple  dans  le  cas  simple  du  mouvement  harmonique 
d'un  profil  en  courant  plan,  pour  lequel  la  frequence  reduite  suffit  a  definir  cette  Evolution. 


Fig.  1.2 

Les  boucles  represented  le  coefficient  de  force  normale  Cm  en  fonction  de  a  ,  mesure  &  dif^drentes 
frequences  (  k«  nff/V  ,  2h  )  sur  un  profil  anime  d'oscillations  forcees  de  tangage  de  grande 

amplitude  (  A  a  *15')  autour  de  1' incidence  moyenne  de  decrochage  statlonnalre. 
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-  Ia  figure  I.  3  ,  tir^e  de  la  reference  [32]  repr^sente  les  cycles  d'hyst^rdsls  des  coefficients 

rodynamlques  transversaux  Cn  .  Cy  ,  Cf  en  fonction  de  l'angle  de  dkrapage  $  ,  raesur^s  en 

statlonnaire  ( n=0>  en  soufflerie  a  une  incidence  sup^rieure  k  celle  du  ddcrochage,  choisie  dans 
une  zone  de  decroissance  ou  d'dclatement  du  tourbillon  d'apex. 

-  Les  deux  exemples  suivants  representent  les  cycles  d'hystkrdsis  obtenus  en  rotation  unifonne  par 
valeurs  disci4tes  et  constantes  de  Q  ,  soit  ! 

SI* a  [p*  0  0 ]  et  0:0  avec  p*m  p*  , 

les  axes  x.  et  x,  ^tant  confondus. 


ng.  i.4 

La  figure  1.4  ,  tir^e  de  la  r^fkrence  [33]  raontre  Involution  du  coefficient  de  moment  Cf,  autour  de 

l'axe  de  rotation  6  X,  d'une  maquette  ADAV  d'allongement  7,  disposke  _ur  la  balance  reprksent^e  fig. 
k  une  incidence  correspondent  k  un  kcoulement  dkcollk  (  n«  1 6°30’,  0°)  • 

L' Evolution  de  ce  coefficient  Cj,  ,  quasi -syroktrique  par  rapport  k  l'origine  peut  etre  analyske  dans 
le  domalne  de  valeurs  positives  de  la  Vitesse  angulaire  rkduite  o‘=%b/2V  .  Les  flkches  lndiquent  le 

sens  d'exkcution  des  mesures  effectukes  d'abord  par  valeurs  crolssantes  de  fL  ,  k  partir  de  la  valeur 
nulls  Jusqu'k  la  valeur  maximale  puis  en  revenant  k  la  valeur  initiale  par  valours  dk crolssantes. 

La  boucle  d'hystkrksis  OABCDE  prksente  des  segments  de  courbe  continus  pouvant  §tre  assimilks,  avec  une 
approximation  acceptable  en  premikre  analyse,  k  des  segments  de  drolte  de  rofene  pente. 

Le  premier  segment  passe  par  l'origine  ;  sa  pente  fypo  est  constante  et  negative,  soit  : 

^jjp^  -  cste  d'ou  done  ^lg^~ ®  113118  l'lntervalle  0  ^  £  ^0fl7S  . 

Lorsque  p£  attelnt  en  A  une  valeur  volslne  de  0,075,  ^/o  augmente  brusquement,  sa  valeur  passe  de 

-  0,05  k  0,01.  Cet  accrolssement  tradult  vraisemblablement  un  dkcollement  de  l'kcoulement  sur  la  partie 
droite  de  l'alle,  provoquk  par  l'lncldence  induite  par  la  rotation  dont  la  variation,  lineaire  sulvant 
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l'envergure  est  crolssante  sur  la  partie  droite  de  l'aile  si  est  posltif.  A  une  certaine  incidence 

locale  un  d^collement  se  produit  entralnant  une  discontinuity  de  la  portance  locale  qui  impose  un  ychelon 
de  moment  positif  par  rapport  k  l'axe  de  rotation  (  Cf0">  0) 

Les  coefficients  du  moment  dfl  au  dycollement  et  le  coefficient  d'aroortiS3ement  C/po  restent  senslble- 
ment  constants  dans  l'intervalle  :  8  -* C 


dans  lequel  :  CfcS  0,1  +  crpo(0)Po*  . 

revyiant  ainsi  l'etablissement  d'un  equilibre  approxlmatif  du  regime  d'ycouleroent  dScoliy .  Lorsque  le 
sens  de  variation  de  la  rotation  s' inverse  en  C,  on  observe  un  nouvel  accroissement  de  plus 

important  que  le  pry cy dent.  Cette  nouvelle  discontinuity  CD  rysulte  de  la  brusque  variation  du  rygime  de 
decollement,  liye  au  changement  de  signe  de  l'acoyiyration  du  roouvement. 


Le  nouveau  rygime  yvolue  continQment  en  fonction  de  pj  et  le  coefficient  d'amortlssement  de  moment 
reste  sensiblement  constant  Jusqu'en  E,  a.  une  vitesse  de  rotation  assez  faible  pour  qu' intervienne  le 
reco1 lement .  Celui-ci  s ' effectue,  par  palters,  jusqu'A.  annulation  de  la  vitesse  de  rotation. 

II  faut  noter  que  ces  mesures  stationnaires  effectuyes  a  et  a  constants  ryvelent  incldem- 

ment  1' influence  de  1' acceleration  de  rotation  sur  le  Cfo 
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-  la  figure  I.  5  ,  tiree  de  la  reference  [34]  non t re  1  '^volution  du  coefficient  de  force  transversale  Cy 

en  fonction  de  sur  une  maquette  de  missile  en  rotation,  calee  k  1' incidence  de  90°  sur  le  montage 

represente  fig.  All. 12  L'ogive  de  cette  maquette  est  constituee  par  un  hemisphere  raccorde  a  un  paral- 
lelepipfede  droit  de  section  carree  a  aretes  arrondies. 


Les  mesures  effectuyes  k  des  nombres  de  Reynolds  R,  variables  dans  le  rapport  de  1  a  10  confirment 
l'importante  influence  de  ce  parametre  sur  la  nature  des  decollements. 

Elies  montrent  qu'un  essal  effectue  a  une  valeur  de  R«  trop  faible  risque  de  donner  une  yvolutlon  du 
torseur  R  completement  dlffyrente  de  celle  qui  se  presenteralt  sur  l'apparell  rdel. 

En  effet,  si  l'on  examine  les  evolutions  de  Cy  mesurees  par  valeurs  crolssantes,  puis  decroissantes 

de  p*  on  constate  &  partir  de  l'origine  : 

ro 

-  dans  la  zone  de  R.  la  plus  yievee,  une  yvolutlon  lincaire  de  Cy  dans  le  domaine  de  ()*  explore. 

-  au  plus  faible  R«  ,  Cy  d'abord  croissant  avec  p*  ;  vers  £2  0,04  ,  un  deuxieme  rygime  apparait, 

entratnant,  une  discontinuity  positive  AB  de  Cy  suivie  d'une  decroissance  BC  sensiblement  lincaire  ;  les 
valeurs  de  Cy  obtenues  par  valeurs  dycroissantes  de  p*  sont  les  memes  que  celles  obtenues  par 
valeurs  crolssantes, 

-  a  la  valeur  de  R»  intermydiaire,  Cy  dycrolt  lineairement  de  0  i  D  suivant  une  droite  de  mfeme  pente 
que  celle  observee  au  nombre  R«  plus  dlevy  ;  puis  en  D  une  discontinuity  de  Cy  sensiblement  dgale  k 
la  precedente  se  manifeste  ;  Cy  dicrolt  ensulte  suivant  la  droite  EC  pratiquement  confondue  avec  celle 
suivie  au  faible  R«  ;  de  C  en  E,  Cy  reste  sur  cette  me  me  droite  Jusqu' a  ce  qu'apparaisse  en  F  une 
discontinuity  FG  ramenant  au  rygime  Initial. 


Un  cycle  d' hysteresis  DEFG  exiote  done  a  cette  valeur  de  R«  . 

-  L'auteur  de  la  ryfyrence  C34]  associe  les  diccontinuitys  observyes  aux  faibles  nombres  de  Reynolds  A 
1’ extension  brutale  de  la  zone  decoliye  en  aval  de  la  maquette  vers  la  face  opposye  A  la  rotation  ;  le 
retour  k  la  configuration  initiate  se  faisant  avec  un  certain  retard  aux  R«  intermydiaires. 

Cette  extension  du  sillage  et  le  cycle  d'hystyrysis  qui  l'accompagne  disparalt  progressivement  lorsque 
R«  crott  (voir  les  schymas  reprysentys  sur  la  partie  droite  de  ]a  figure  1.5). 


1.6. 2. 2-  Phdnomenes  aerodynamiques  observes  sur  les  avions  au  eours  d'une  manoeuvre  &  grande  incidence 
en  ecouiement  de colie  stationnalre  ou  non. 

La  figure  1.6  .  extraite  de  la  reference  [353,  represente  le  coefficient  de  portance  Cj  *  CL 

equilibre  en  fonotion  de  l'incidenoe,  mesure  a  differents  nombres  de  Mach  sur  une  maquette  au  1/10 
de  l'ADAV/ADAC  HARRIER,  &  derapage  nul. 
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Une  discontinuite  de  la  pente  C*,  («  CLa)  est  visible  sur  chaque  courbe  a  une  Incidence  critique  qui 
decrolt  de  11°  a  8°  lorsque  le  nombre  de  Mach  crolt  de  0,6  k  0,9.  Dans  ce  domalne  d'incidence  des  decol- 
lements  commencent  a  apparaltre  puis  s'etendent  sur  une  partie  de  plus  en  plus  importante  de  l'aile  au 
fur  et  a  mesure  que  1' incidence  et  le  nombre  Mach  augmentent.  Ces  decollements  &  caractere  instationnaire 
sont  la  cause  des  perturbations  des  qualites  de  vol  apparaissant  a  incidences  croissantes  avant  que  ne  soit 
attetnt  le  Ci  maximal,  dans  l'ordre  suivant  : 

1)  Le  tremblement  (buffeting)  au  cours  duquel  le  torseur  R  reste  generalement  symetrlque,  mais  pr^sente 
une  premiere  phase  de  variations  d'intensite  relativement  falble  dont  le  spectre  en  forme  de  bruit  blanc 
excite  les  modes  de  structure,  et  une  seconde  phase,  a  caractere  instationnaire  nettement  plus  marque  et 
d'intensite  plus  importante  qui  s ' apparenterait  a  un  sillage  de  tourbillons  altem^s  avec  une  loi  de 
frequence  aleatoire .  [36] 

2)  Le  balancement  d'alle  (wing  rock),  oscillation  transversale  de  roulis,  lacet  et  ddrapage  d'amplitude 
importante,  partieulierement  en  roulis  sur  les  avions  de  faible  allongement,  accompagnee  d'une  liSgere 
composante  longitudinale .  Ce  mouvement  est  dfi  vraisemblablement  a  un  processus  de  d^collement-recollement 
presentant  un  cycle  d' hysteresis  noteur,  analogue  a  oelui  analyse  au  paragraphe  11.5.  2.  2  ,  du  second 
chapitre  mais  sous  une  forme  antisymetrique  dans  ce  cas  particulier. 

3)  L' inversion  de  la  commande  de  gauclilssement  (roll  reversal  )  qui  peut  s' interpreter  comme  la  conse- 

quence  d'une  variation  de  la  difference  des  efficacites  des  gouvemes  droite  et  gauche,  operant 

chacune  dans  un  domaine  non  lineaire  different  : 

Pour  6l>  o  on  obtient  suivant  les  cas  ACt>  0  OU  ACf  <  0  :  ceci  re  produit  par  exemple 

si  une  zone  dccollee  s'etend  syst^matlquement  sur  les  gouvemes. 

4)  L'lnstablllte  de  direction  et  le  depart  en  vrllle  lncontrolable ■ 


Des  indications  sur  la  determination  de  ces  perturbations  dans  les  installations  au  sol  sont  pr/ sente es 
ci-dessous  en  II .  5  .  2 
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CHAPITRI  II 


TECHNIQUES  D' EXPERIMENTATION  AU  SOL 


Q.l.ffiKERALTIES. 

Dans  son  Important*  oonnunl cation  au  symposium  FDP/AQARD  sur  les  PARAMETHES  EE  STABUJTE,  tenu  en  Mai  78 
A  A  thine  8,  ORUK-RUCKEWANN  a  prAsentA  les  Installations  de  mesure  des  dArivAes  aArodynamlques  en  souf¬ 
fle  rle,  utlllsAseA  l'heure  actuelle,  dans  la  plupart  des  centres  d'essals  aArodynamlques  d'AmArique  du 
Nord  et  d' Europe.  Des  Informations  supplemental  res  sur  des  Installations  mlses  en  service  ces  demlers 
mols,  ou  encore  en  cours  de  mlse  au  point,  ont  AtA  apporties  par  d'autres  auteurs  A  ce  mSroe  symposium. 


Ces  comaunl cations  montrent  que,  pratlquement,  ces  mesures  sont  effectuies  dans  des  soufflerles  classlques 
sur  des  maquettes  anlmAea  de  petlts  mouvements  dont  le  nombre  de  liberty  est  restrelnt. 

NAanmoins,  dans  certains  caa  partlcullers,  les  valeurs  de  ces  coefficients  sont  obtenues  par  des  mesures 
effectives  sur  le  mouvement  d'une  maquette  en  vol  llbre . 

On  peut  aussl  clter,  pour  mi moire,  une  tentative  de  mesure  des  dArivAes  de  R  par  rapport  A  12  sur 
maquette  fixe  dans  une  soufflerie  dont  l'Acoulement  est  rendu  artificiellement  rotatlonnel  [371. 

II  a  sentoJji  superflu  de  reprodulre  dans  leur  ensemble  les  descriptions  des  montages  et  1'exposA  des  mAtho- 
des  expirlmentales  dAjA  rassemblAs  dans  les  publications  du  symposium  PDP/AGAHD  d'Athines,  plus  partlcu- 
llirement  dans  celle  d '  ORUK-RUCKEMANN . 

C'est  pourquol  quelques  a chi mas  seulement  d1 installations  typiques,  4e s  indications  sur  les  ordres  de 
grandeur  et  les  risultats  d'essals  sont  prAsentAs  lei  A  titre  d'exemple.  Le  dAveloppement  de  ce  chapitre 
est  essentlellement  un  rappel  des  princlpes  et  des  mAthodes  d' experimentation,  prAsentA  dans  l'ordre 
sulvant  : 

f  11.2  Toraeurs  d 'efforts  sur  une  maquette  en  soufflerie. 
f  11.3  Princlpes  des  installations  d'essai. 

f  11.4  datura  des  Informations  (toraeur  d'efforts  et  vecteur  d'etat)  et  chalnes  d' acquisition. 

|  11.3  Procedures  expArimantalea  et  traitement  des  Informations. 

11. 2.-  TORSEURS  D' EFFORTS  SUR  UNE  MAQUETTE  EN  SOUFETERIE. 

Ia  maquette  est  dlapoaAe  dans  la  soufflerie  sur  un  montage  (ou  suspension)  qui  sert  soit  &  1 ' Immobiliser, 
aolt  A  l'anlmer  d'un  mouvement  A  un  ou  plusleurs  degrAs  de  libertA. 

L'Aquation  du  toraeur  d'efforts  sur  la  maquette  dans  le  repine  Gxyz  llA  A  celle-ci  a  pour  expression 
A  tout  Instant  t  : 

(  1  )  ^-R+Gj  +  Sj  , 

dans  laquelle  Sg  eet‘  le  toraeur  d'effort  ImposA  A  la  maquette  par  la  suspension. 

L'Aquation  prAcAdente  est  A  rapprocher  de  1 'Aquation  (1)  du  §  1.4.1  Atablie  dans  le  cas  du  vol  llbre. 

Solt  un  petit  lntervalle  de  temps  it  dAflnl  par  tf  — 1|,  la  diffArence  entre  les  valeurs  de  (1)  aux 
Instants  tf  et-  t|  s 'exprime  par  : 


(  2  )  m(igf_Eg,)  =  Rf-Ri  +  AGg-A^j+Sff-Sfj  . 


IL2  .2~Sl  le  toraeur  Rest  contlnu.  on  peut  Acrlre , A(  )  reprAsentant  une  variation  arbitralrement  petite  : 

*«f“*«l*  ,  R,-R,  ,  SM-SxjsiSx  • 

31  R  est  dArlvable  par  rapport  A  Es  et  fe»  (2)  s'Acrit  : 

(  3  )  ^ At  +*,J?AE  J+’AGg+iSg  . 

La  relation  entre  et  est  donnie  en  annexe  AI.5,  les  expressions  de  It  et  de  A Gg  sont  lndi- 

quAes  au  §  1.4.1. 
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« 

En  faisant  passer  les  tenses  en  CJpg  et  CM  au  premier  Mmbre,  il  vient  : 

(  4  )  CreAE*+CR(*\E*  *^k-’[n»AEJI-AG1,-ASl]-«^AE*, 

qui  s'ecrit  sous  forme  transpose  : 

(  5. )  AI*C^e + H*C^I  =  [I^m  -  AG,-  AsJk-'-Ai*  *•, , 

2 

■1  ,  le'et  •,  etant  symetriques  ont  les  mgmes  expressions  que  leurs  transpos^es. 

Liquation  (5)  peut  etre  consld^ree  oomme  un  svstfeme  d' Equations,  dont  les  inconnues  sont  les  vecteur3 
colonnes  des  matrices  transposes  ttCpJ  . 

Cft  dtant  suppose  derivable  par  rapport  a  E*  et  t*  ,  les  valeurs  de  AE  et  At  peuvent  etre  malnte- 
nues  dans  des  intervalles  suffisamment  petite  pour  que  Cr.e  et  Crje  reprdsentent  avec  une  bonne  approxi¬ 
mation  les  constantes  du  systfeme  lln^alre  tangent  du  oas  de  vol  deflni  par  Eg  et  Eg  • 

Les  valeurs  de  oes  constantes  ne  sont  done  pas  ndeessairement  les  me  me  s  dans  tout  le  domaine  de  valeurs 
de  Eg  et  Eg  . 

Dans  oes  conditions,  un  systeme  (5)  donnd  comprend  2  fois  36  inconnues  et  theorlquement  pour  le  resoudre 
il  faudralt  pouvoir  disposer  de  2  fois  six  ensembles  de  valeurs  deAE*et  Afllndairement  indd pendants . 
Cependant,  pour  les  raisons  indiquees  aux  paragraphes  14.5  etl4.6,un  grand  nombre  d ' elements  de  CRiE  et 
CR  i  peuvent  etre  negliges  dans  1' etude  de  la  stabilite  des  avions  et  des  missiles.  lies  prinoipes  des 
metliodes  experimentales  de  determination  des  coefficients  qui  interviennent  effectivement  sur  la  stabilite 
sont  exposes  au  paragraphe  3  de  ce  chapitre . 

Cette  determination  ne  p  re  sente  pas  de  diffioulte  theorique  parti culifere,  car  en  principe,  on  peut  conce- 
voir  une  experimentation  qui  donnerait  les  ensembles  de  vecteurs  AE*et  At  sous  la  forme  de  matrices 
diagonales,  ce  qui  permettralt  de  oalculer  un  par  un  les  elements  de  CRj:  et  CM.  C'est  pr^cisement  le 
cas  des  valeurs  des  derivees  aerodynamiques  deduites  des  essais  stationnaires  classiques,  par  exemple : 

Cgo  ,  Cma  ,  Cyp  ,  Cfp  ,  Cnp 


11.2.2-  Dans  le  cas  ofi  le  torreur  R  presente _ soit  des  discontinuites ,  soit  des  variations  trfes  rapides 
(  d^crochage )  l'equation  (5)  n'a  pas  de  sens,  raais  l'equation  (l)  reste  valide  quelles  que  soient  lss  valeurs 
de  Eget  Eg  . 

Les  valeurs  des  sauts  de  Or  ,  011  de  CRE  et  CR£  symbolisees  par  : 

<cR),CcR.E)  .(cm)  (page  36  §  II. 1.5  de  la  reference  2  indiquee  en  annexe  A.I.i  ) 
peuvent  etracalculees  d'apres  cette  equation  si  *7  ,  Gg  et  Sg  sont  connue  dans  un  aomaine  de  valeurs 
de  Eg  et  Eg  incluant  celles  des  discontinuites. 

Ce  mode  d' interpretation  peut  etre  applique  aux  mesures  comportant  des  variations  a  travers  des 

zones  de  discontinuite . 

Dans  tous  les  cas  ou  CreCrI  ne  sont  pas  constants,  la  methode  precedente  correspond  A  une  approxi¬ 
mation  "multilineaire",  s^il  y’a  discontinuite  elle  s'accompagne  generalement  d'un  cycle  d'hyseresis 
[  ref. 25  page 7  T.  1  §l,2].et  §1.6.2. 

II  .3  PRUl'IPES  DES  INSTALLATIONS  D'ESSAI 
ll.3.l-Hesures  stationnaires 

Elies  sont  effectuees  sur  les  balances  de  mesures  stationnaires  rappelees  icl  pour  memoire  [38]  et  sur  les 
balances  en  rotation  uniforme  parallele  a  1 'ecoulement  dont  quelques  schemas  sont  presentes  en  annexe  A 
titre  d'exemples. 

On  remarque  que  certaines  de  ces  balances  sont  mises  en  mouveroent  par  des  moteurs  hydrauliques  ou  A  air 
comprimd .  L'interet  de  cette  disposition  est  precise  ,1-dessous  en  II. 3.3.3 

(1. 3.2-Osclllatlon3  llbres 

Deux  categories  de  montages  sont  utilises  pour  anlmer  une  maquette  de  mouvements  d'oscillations  "llbres". 

-  les  degres  de  llberte  11.  'r4s  avec  les  montages  de  la  premiere  categorle  sont  effectivement  llbres  sans 

contralnte.  ce  qui  peut  s' exprimer  par:  _  — 

De  tels  montages  ont  ete  mis  en  service  il  y  a  plus  de  dlx  ans  [  3  9  ]  dans  les  grandes  soufflerles  super- 
sonlques  et  hypersoniques  de  l'A.E.D.C.  k  TUUAHOMA .  Les  axes  de  rotation  sont  materialises  par  des  paliers 
A  gaz  dont  la  valeur  du  coefficient  de  frottement  se  sltue  au-dessous  de  tout  seuil  mesurable.  Des  montages 
4  un,  deux  et  trols  degres  de  llberte  sont  en  service. 

-  Ceux  de  la  seconde  categorle,  dont  deux  modeles  figurent  comme  exemples  fig.  AIM  ,  sont  constitues  par 
des  articulations  eiastiquco  k  lanes  crones.  Ce  genre  d'articulatlons  ne  rc  defotme  pratlquement,  qu' 
autour  de  l'axe  represente  par  1 ' Intersection  du  plan  median  des  lames,  par  l'action  d'un  moment  proportion- 
nel  k  1' angle  de  deformation.  La  resistance  mecanique  de  1' articulation  Impose  toutefois  de  ne  pas  depasser 
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suivant  les  dimensions,  des  angles  de  valeur  maximale  comprise  entre  1  &  6  degree  ,  One 

maquette  dispos^e  sur  un  tel  montage  prend,  k  la  suite  d’un  4  cart  de  sa  position  d'kqulllbre,  un  mouve- 
roent  d' oscillations  ik  un  degry  de  liberty,  dits  libres,  mais  en  r^alit^  soumises  k  la  contralnte  due  k 
la  deformation  de  sorte  que  : 

A*,*  0  • 

Ia  valeur  du  coefficient  d'amortissement  de  ces  articulations,  quoique  nettement  supkrieure  k  celle  des 
pallers  k  gaz,  est  encore  trfes  petite  de  l'ordre  de  : 

0,0005  <  £,<  0,0015 

d'aprks  les  resultats  obtenus  par  X.VAUCHERET,  repr^sent^s  fig. AH 

Elies  sont  done  parfaitement  adaptees  k  la  mesure  des  falbles  amortissements. 

Des  dispositifs  de  conmande  k  distance  permettent  d'effectuer  des  lSchers  k  une  cadence  sufflsante  pour 
obtenir,  pendant  la  duree  de  1  'ecoulement  sain  d'une  soufflerie  interim  ttente  plusieurs  series ji' oscil¬ 
lations  k.  une  moyenne  donnee,  ou  me  me  k  des  incidences  dlf  fkrentes .  Un  schema  de  1 '  installation  utlllske 
dans  les  souffleries  de  l'ONERA  a  Modane  est  indiquk  fig .  All.  3  • 

Xa  figure AKZirepresente  une  variante  d'un  montage  de  cette  eategorie  sp^cialement  dtudiy  par  le  P.F.A. 
pour  la  mesure  de  l'amortissement  d'une  aile  ogivale,  dans  sa  soufflerie  transsonique  supersonique  S4 
de  section  de  veine  1  rr£. 

Ilj.3  -Oscillations  entretenues  et  forc^es 

Les  principes  de  ces  deux  modes  de  mise  en  mouvement  de  la  maquette  se  distinguent  entre  eux  de  la 
fagon  suivante  : 

IL3.3.1  -La  maquette  animee  d'oscillations  entretenues  est  disposee  sur  une  suspension  elastlque  dont  la 
rigidite  est  calculee  de  fagon  k  s 1  assurer  d'au  mo ins  une  frequence  de  resonance .  Uh  moteur  relid  k  la 
suspension  par  une  liaison  de  faible  rigidite,  imprevue  a  la  maquette  un  mouvement  d'oscillations  de 
frequence  voisine  de  cette  frequence  de  resonance. 

Dans  ces  conditions,  la  puissance  necessaire  a  l'entretlen  du  mouvement  est  gendralement  minimale.  II 
arrive  meme  que,  lorsque  les  conditions  correspondent  k  la  divergence  ou  au  flottement  sont  rdunies, 
cette  energie  soit  empruntee  a  l'ecoulement  de  la  soufflerie.  Le  moteur  doit  alors  etre  ddbraye,  ou 
meme  le  cas  echeant  pouvoir  fonctionner  en  recepteur  pour  s'opposer  k  la  divergence. 

■  Des  schemas  de  balances  du  R.A.E.  et  du  N.A.E.  sont  presentees  en  annexe  k  titre  d'exemples  fig. 
et  fig.  . 

Les  moteurs  de  ces  balances  sont  des  osclllateurs  electromagnetiques.  Le  dessln  de  la  balance  du  R.A.E. 
a  ete  etudie  de  fagon  k  s'assurer  de  l'existence  de  trois  modes  d'oscillations  particuliepes  pratique- 
ment  decouples  entre  eux  k  leur  frequence  de  resonance.  Cette  disposition  permet  de  determiner  les 
derivdes  aerodynamiques  transversales  sur  une  maquette  animee  de  deux  modes  d'oscillations  en  ddrapage 
et  en  lacet  autour  de  deux  axes  paralleles  entre  eux  et  perpendiculaires  a  l'axe  du  dard,  puis  d'un 
mode  de  rculis  autour  de  celui-ci. 

Ie  passage  d'un  mode  k  1 'autre  est  obtenu,  sans  necessiter  aucun  demontage,  en  faisant  unlquement  varler 
la  frequence  du  courant  d' alimentation  de  l'oscillateur  electromagnetlque . 

lie  N.A.E.  a  realise  et  perfectionnd  deux  modeles  de  balances-dard  en  oscillations  k  un  seul  degry  de 
liberty  en  rotation.  Sur  l'une  d'elles,  la  maquette  est  animee  d'oscillations  autour  d'un  axe  perpen- 
diculaire  k  l'axe  du  dard  correspondent  a  des  raouvements  de  tangage  ou  de  lacet  suivant  la  disposition 
de  la  maquette  sur  la  balance.  Sur  l'autre,  la  maquette  est  animee  d'oscillations  de  roulis  autour  de 
l'axe  du  dard. 

Dans  les  realisations  les  plus  recentes  1'attache  de  la  maquette  sur  la  partie  mobile  de  la  balance  est 
faite  par  1' intermediaire  d'un  dynamometre  a  faible  course  k  cinq  composantes  t  14  ]  ;  la  composante  X 
suivant  l'axe  du  dani  n'est  pas  mesurde. 

Les  deux  balances  du  NAE  sont  spdcialement  destinees  aux  mesures  k  grande  incidence.  Elies  donnent  les 
elements  des  matrices  et  CR,E  qui  ont  une  influence  sensible  sur  la  dynamlque  du  vol.  II  dolt 

cependant  etre  note  que  les  rdsultats  fournis  ne  sont  valables  que  dans  les  domaines  ou  le  torseur  CR 
est  continument  derivable . 

11.3.3.2-  Le  prlnclpe  des  balances  en  oscillations  forcees  consiste  a  imprimer  direotement  ,  k  une 
maquette  montee  sur  une  suspension  aussi  rigide  que  possible,  un  mouvement  d'oscillations  k  un  ou 
plusieurs  degre'e  de  liberty. 

L'attache  de  la  maquette  sur  la  suspension  est  faite  par  l'intermddialre  d'un  dynamomktre  k  faible 
course.  Les  schemas  de  la  plupart  de  ces  balances,  dont  certaines  sont  en  service  depuls  pres  de  vingt 
ans,  sont  presents  dans  la  rticente  publication  d ' ORUK-RUCKEMANN  deja  citee  en  II.  1  .  II  a  paru 
n^anmoins  utile  de  rappeler  k  titre  d'exemple  la  balance-dard  transsonique  ou  supersonique  du  CALSPAN, 
mise  en  service  en  1962.  qui  pry sente  1' originality  de  pouvoir  mettre  la  maquette  en  oscillations  de 
tangage  d'amplitude  importante  (  1  6°)  autour  d'un  axe  virtuel  confondu  avec  l'axe  transversal  0  de 
l'akrcnef  fig.AH.5 .Cette  disposition  offre  des  perspectives  intyressantes  k  l'heure  actuelle  pour  lVtude 
du  dycrochage  aux  grandes  incidences .  [S9J , 
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Dans  les  comptes  rendus  du  m§me  symposium,  trois  balances  de  stability  (rotation  unlforroea,  oscillations 
entretenues  et  oscillations  forcdes)  basse-vitesse  (  V 2 60  ms-1)  sont  present des  dans  les  communications 
du  Pr.  HA  PER  [  41  ]  et  de  Jan  derHSCKEN  et  ali.[  42  ) .  Ces  balances  sont  mises  en  mouvement  respectivement 
par  un  moteur  hydraulique,  un  oscillateur  4lectromagn4tique  et  des  moteurs  41ectriques  A  courant  continu. 

11.3.3.3-  Remarques  sur  la  motorlsatlon  des  balances. 

Les  progrks  accomplis  ces  dlx  demikres  ann4es  dans  le  d^veloppement  technique  des  g4n4rateurs  et  moteurs 
k  fluides  sous  pression  (pneumatique  et  hydraulique)  ont  conduit  k  utiliser  ces  raat^riels,  fabriqu4s 
industrlellement,  pour  la  motorisation  des  balances  en  rotation  et  en  oscillations. 

Ces  moteurs  p re sen tent  les  avantages  suivanta  : 

-  foumir  les  me  roes  puissances  que  les  moteurs  4lectriques  avec  des  masses  et  des  encombrements  beaucoup 
moins  importants. 

-  ne  pas  introduire  de  parasites  dlectriques  au  voisinage  des  lignes  conductrices  des  signaux  de  faible 
niveau  (  ljiV  k  quelques  millivolts)  4mis  par  les  divers  ddtecteurs. 

On  observe  4galement  que  le  pilotage  4lectronique  assure  a  l'heure  actuelle  avec  une  grande  precision, 
soit  la  regulation  de  la  vitesse  de  rotation,  soit  la  stability  de  la  frequence  et  ceci  quelle  que  soit 
la  source  d'4nergle  utilises. 

L'utilisation  de  l'alr  comprime  paralt  convenir  aux  moteurs  en  rotation  uniforme  de  faible  puissance, 
de  1 ' ordre  de  2  kw  ,  et  l'hydraulique  aux  puissances  sup4rieures  fig.AD«12- 

En  outre,  dans  le  cas  particulier  des  oscillateurs  hydrauliques  asservis,  1 'amplitude  et  la  position 
moyenne  des  oscillations  sont  reglees  continQment  a  distance. 

A  notre  connaissance  la  premifere  balance-dard  mue  par  un  vdrin  hydraulique  asservi  a  eter4alis4e  par  l’AEDC. 
Les  figures  de  la  page  171  de  la  reference  [10],  montrent  le  v4rin  et  sa  servovalve  disposes  a  l'extr4- 
amont  du  dard  imm4diatenient  en  aval  d'une  articulation  dont  l'axe  est  perpendiculaire  a  celui  du  dard. 

Les  faibles  dimensions  de  ce  dispositif  pemettaient  de  le  loger  k  l'interieur  du  fuselage  de  la 
maquette.  Cette  balance  a  4te  perfectionnee  depuis  comme  le  montre  la  fig.  9  de  la  reference  [  14  ]. 

Une  balance,  p revue  pour  etre  instance  dans  les  grandes  souffleries  subsonlque  de  l'ONERA  k  Modane 
(8m  de  diamitre)  et  au  FAUGA  (3,5  x  4,5  m2),  utlllsant  ce  mode  de  motorisation  est  en  cours  de  montage. 

Elle  comporte  deux  tetes  4quipees  de  dynamometres,  disposees  a  l'interieur  du  fuselage,  qui  peuvent 
etre  fixees  l'une  ou  l'autre  a  I'extr4mit4  d'un  mkt  vertical.  L'une  d'elles  est  destinee  aux  mesures 
stationnaires  et  l'autre  aux  mesures  instationnaires .  Ia  difference  essentlelle  entre  ces  deux  dispositifs 
reside  essentiellement  dans  les  valeurs  maxlmales  des  vitesses  angulaires,  de  l°/sec.  en  mesures  station¬ 
naires  et  15°/sec.  en  essais  instationnaires. 

Les  conditions  de  fonctionnement  de  cette  balance  sont  approximativement  : 

maquette  :  masse  m  S  300  kg  ;  inerties  ;  |j  ^  50  kg  m2  ,  lyet  lz  S  100  kg  ire, 

envergure  b  <  }  *  , 

plages  d' angles  d' Incidence;  de  derapage  et  de  gite  : 

-10’  $  a  5  )o‘  i  -  lo’  5  p  ^  180°  ;  -150  a?  O  <  15°  , 

L' amplitude  des  oscillation?  en  tangage,  lacet  et  roulis  est  de  +  2°  pour  un  lntervalle  de  frequence  de 
0,5  £  f  ^2  Hz. 

Enfln.le  mSt  vertical  peut  etre  deplace  suivant  son  axe  a  une  vitesse  de  A  ms  1  sur  une  course  de  3  metres. 
Vol  libre  danr  les  installations  au  r.ol. 

Les  m4 thodes  de  determination  des  coefficients  a4rodynamiques  d'apres  1 'analyse  du  mouvement  d’une 
petite  maquette  de  projectile  en  vol  llbre,  au  tunnel  de  tir  ou  dans  une  soufflerie  supersonique,  ont  4te 
decrltes  dans  deux  AGARDograph  publics  il  y  a  une  dizaine  d'ann4es  [441  et  [«]  .  Une  etude  sur  le  me me 

sujet  poursuivle  au  Laboratoire  de  Recherches  Balistiques  dr  Vernon  est  parue  en  1973  poj 

Les  objectlfs  recherches  dans  ces  moyens  d'essais  4tant  pratlquement  atteints,  les  efforts  de  modernisation 
de  ces  demleres  ann4es  ont  4t4  orl entes  princlpalement  sur  les  4quipements  des  installations  d'essais  au 
sol  sp4clalls4es  dans  l'4tude  du  vol  libre  d'une  maquette  d'avion,  en  vue  de  r4pondre  aux  questions  posees 
par  les  qualit4s  de  vol  brikvement  exposees  au  paragraphe  II. 5. 2. 

Paml  les  Installations  concem4es  on  peut  citer  notanment  : 

-  Les  souffleries  verticales  de  vrllle  de  la  NASA  k  Langley-Field  et  de  l'lnstitut  de  M4canique  des  Fluides 
de  Lille.  Les  dimensions  et  la  vitesse  de  ces  deux  souffleries  sont  voisines,  des  maquettes  de  1,5  m  d'en- 
vergure  peuvent  y  8tre  essay4es. 

-  L' installation  d'essai  en  vol  libre  de  la  grande  soufflerie  (9m  x  l8m)  de  Langley-Field  dans  laquelle  sont 
4tudl4es  des  maquettes  motoris4es  de  1,5  m  d'envergure  et  de  3  m  de  longueur. 

-  Le  nouveau  laboratoire  de  vol  libre  de  l'lnstitut  de  M4canique  des  Fluides  de  Lille  disposant  d'un 
domaine  de  vol  de  50  m  de  longueur  dans  lequel  les  maquettes  sont  catapul ties . 

Le  nombre  de  Mach  de  l'ecoulement  de  ces  installations  se  situe  aux  environs  de  0,1. 
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II. 4*  NATURE  DES  INFORMATIONS  ET  SYSTEMES  D' ACQUISITION 
Les  informations  recueillies  en  essai  sont  represent6es  par  : 

-  les  elements  &  vecteur  AC*  ,  de  s a  ddrivde  par  rapport  au  temps  At*  et  ceux  du  torseur 

dans  les  cas  ou  et  cm  existent  et  sont  calculdes  d'aprks  (5).  c'est-k-dire  dans  le  domaine 

ou  R  est  lindaire.  ' 

-  Les  elements  de  E«  ,  E»  et  Eg  en  fonction  du  temps  dans  les  autres  oas,  en  partioulier  dans  les 
domaines  ou  R  ost  instationnaire  ou  prdsente  des  discontinuitds. 

Dans  la  pratique,  1' acquisition  des  informations  est  faite  en  deux  dtapes  : 

a)  la  premiere  correspond  aux  mesures  stationnaires  classioues.  dans  lesquelles  il  faut  induce  les 
mesures  sur  balance  en  rotation  uniforme  d'axe  parallele  k  1' dcoulement.  Ces  mesures  donnent,  d'apres 

1' Equation  (l)  considdrde  dans  le  cas  &  S0  ,  les  vecteurs  R  et  par  consequent  CR  dans  un  domaine 
de  valeurs  de  E*  d1 expression  transposee  : 

I*=[/M  p'q*  r*]  E  =  fi  A V  =  0  H*=[p*q*r*] 

Les  informations  CR  en  fonction  de  E*  ainsi  obtenues  permettent  : 

-  de  calculer  d'apres  (5)  les  Elements  de  ^R  E  gradient  de  CR  suivant  E*  dans  les  parties  de  ce 
domaine  oh  II  est  continu  et  derivable  ; 

-  d' observer  et  d'evaluer  des  discontinuites  par  exemple  des  cycles  d* hysteresis  visibles  sur  les  figures 

1.3  et  I*  A 

b)  La  seconde  etape  comporte  1* acquisition  des  m@mes  .informations  en  fonction  du  temps  sur  une  maquette 
en  mouvement  tel  que  l'un  au  moins  des  elements  de  E  soit  different  de  zero  ;  en  outre  si  et  w 
n'ont  plus  ndcessairement  la  m$me  direction.  Ces  conditions  s'expriment  par  : 

E*0  IflAVl^o 


Une  vue  d'ensemble  des  systemes  d'acquisition  utilises  actuellement  dans  les  installations  au  sol  est 
prdsentee  dans  lq  suite  de  ce  paragraphe  sous  leurs  deux  aspects  :  detection  et  conditionnement- 
enregistrement. 

11.4.1“  La  detection  des  forces  est  effectuee  de  faqon  classique  depuis  plus  de  vingt  ans  par  des 
barreaux  d'acier  dquipds  de  jauges  de  contrainte,  gdndralement  en  silicium  pour  les  mesures  instation- 
naires . 

-  Le  ddveloppement ,  plus  recent,  des  techniques  d'asservissement  a  largement  contribud  a  la  creation  et 
k  la  fabrication  industrielle  de  nombreux  detecteurs  de  grandeurs  cindmatiques  (position,  vitesSe, 
acceleration) . 

On  peut  citer  k  titre  d' exemple  : 

-  les  resolveurs.  detecteurs  de  position,  constituds  par  de  petits  transformateurs  a  secondaire  mobile, 
alimentds  par  un  courant  alternatif  de  400  Hz.  Les  clinometres  fabriques  suivant  ce  principe  ont  une 
plage  de  mesure  de  0  a  360°,  dans  une  bande  de  frequence  variant  du  continu  k  plusieurs  dizaines 

de  Hz. 

-  Les  gvrometres  miniaturises  qui  peuvent  Stre  loges  dans  les  maquettes  de  vol  libre. 

-  Les  acceldrometres  utilises  en  grand  nombre  dans  les  etudes  de  resistance  des  structures  aux  vibrations. 
Leurs  elements  sensibles  sont  gdndralement  des  cristaux  de  silicium  pidzo  rdsistants  ou  des  cdranui- 
aues  piezo  eiectriaues.  La  borne  inferieure  de  leur  bande  passante  part  du  continu  pour  les  premiers 

et  est  de  l'ordre  du  Hz  pour  les  seconds,  alors  que  la  borne  superieure  peut  depasser  plusieurs 
dizaines  de  kHz  avec  l'un  ou  1' autre  element.  Les  dimensions  de  ces  elements  miniaturises  permettent 
de  les  installer  sans  difficulte  k  l'intdrieur  des  ailes  et  des  empennages  des  maquettes. 

Les  signaux  eiectriques  dmis  par  les  detecteurs  sont  presque  toujours . transmis  aux  instruments  de  mesure 
par  des  conducteurs  de  dimensions  aussi  reduites  que  possible,  ndanmoins  des  rdcepteurs-emetteurs  de 
tdldmesure  et  de  tdiecommande  miniaturises  sont  actuellement  embarquds  sur  des  maquettes  essaydes  en 
vol  libre.  Les  sections  des  fuselages  de  ces  maquettes  sont  dgales  ou  supdrieures  k  un  ddcimktre  carrd. 

Un  dtage  de  prdamplification  est  gdndralement  interpose  entre  la  sortie  du  ddtecteur  et  l'entrde  de  l'd- 
metteur. 

Des  chatnes  de  teldtransmission  k  30  voies  sont  installdes  dans  les  maquettes  essaydes  k  l'Institut  de 
Mdcanique  des  Fluides  de  LILLE  [  19  ]. 

II. 4. 2—  La  premikre  fonction  du  condi tionnement  est  d'amplifier  les  tensions  de  sortie  (quelques  pV 
k  100  mV)  des  signaux  dmi3  par  les  ddtecteurs  pour  les  rendre  compatibles  avec  les  entries  des  autres 
organes  du  systkme  d'acquisition,  gdndralement  de  ±  ImV  k  ±  10  V. 

La  sortie  des  amplificateurs  comporte  un  premier  dtage  de  filtrage  qui  dlimine,  en  particulier,  les  bruits 
d'origine  dlectrique  ou  dlectronique. 
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$2nsuite,  suivant  le  mode  de  traitement  adopts,  lea  signaux  ainsi  obtenus  sont  : 

-  soit  enregistrEs  sur  bande  magnetique  analogique  pour  $tre  traitEs  en  temp3  diffErE  ; 

-  aoit  branches  aur  des  calculateurs  analogiques  specialises  ; 

-  3oit  encore  munerises  direct ement  puis  transmis  aux  registres  d'un  ordinateur  ou  enregistrea  sur  bande 
magnEtique  numerique  ;  la  cadence  des  convertisseurs  analogiques-numEriques  modernes  eat  de  40  kHz. 

Lea  ordres  de  commutation  et  d' Echantillonnage  sont  donnes  par  des  circuits  logiques  intEgres  dans  les 
calculateurs  analogiques  specialises  et/ou  par  des  ordinateurs  de  gestion  utilises  pour  1* ensemble  des 
essais  de  la  soufflerie. 

H.  5-  TRAITEKENT  DES  INFORMATIONS  ET  PROCEDURES  EXPERIMENT  ALES 

L'objet  du  traitement  des  informations  est  la  recherche  des  solutions  des  equations  (l)  ou  (5),  ecrites 
au  paragraphe  II.  2  rappeldes  ci-dessous  : 

i  l  )  R+G*  +  s„  , 

(5  )  [^^l-SG*-ASl]k*,-  , 

dans  lesquelles  R  (done  Cr  )  ou  ,  Cjj  g  sont  les  inconnues. 

Les  informations  d' essais  sont  affectees  de  bruits  k  caractkre  plus  ou  moins  aleatoire  et  d' intensity  rela- 
tivement  faible  lorsqu'ils  sont  dus  k  la  turbulence  de  l'ecoulement  de  la  soufflerie  ;  par  contre  leur 
intensity  devient  importante  s'ils  sont  provoquds  par  des  singularity  de  l'ecoulement  sur  la  maquette 
(ddcollements  et  ondes  de  choc  instationnaires). 

L'une  des  fonctions  essentielles  du  traitement  des  informations  est,  par  consequent,  de  depouiller  les 
signaux  utiles  des  bruits  aldatoires  et  si  possible  de  les  evaluer.A  cet  effet,  ces  dernieres  annees,  ont 
dte  ddveloppes  parallelement  des  mdthodes  et  moyens  de  f iltrage  analogiques  et  numeriques. 

Grftce  aux  pro  grey  constants  de  l'informatique,  il  apparalt  maintenant  que  les  moyens  de  traitement  entiere- 
ment  numerique ,  beaucoup  plus  precis,  sont  appeles  a  remplacer  les  moyens  analogiques  ou  hybrides  actuelle- 
ment  en  service,  dds  que  des  temps  de  rdponse  assez  brefs  pour  satisfaire  les  besoins  de  traitement  en  temps 
rdel  auront  dte  obtenus.  Cette  generalisation  ne  sera  effective  qu'aprfes  la  mise  au  point  de  mini-ordina- 
teurs  conqus  spdcialement  pour  une  fonction  particulikre,  dejk  dprouvee  a  1' occasion  des  realisations  de 
traitement  analogique  telles  que  par  exemple  :  le  dispositif  (dampometer)  de  mesure  des  coefficients  d'amor- 
tissement  reduit  en  oscillations  libres  [A  7]  ou  la  commande  automatique  d'une  gouverne  anti-flottement  dont 
le  schdma  est  indique  fig.  1 ,  page  13  de  la  reference  [48]. 

Dans  la  pratique,  les  mesures  sur  des  maquettes  en  mouvement  avec  et  sans  ecoulement  aerodynamique ,  sont 
traitees  le  cas  dchdant  lindairement  d'aprks  liquation  (5)  suivant  des  mdthodes  exposbes  ci-dessous  en 
II.  5. I  .  Dans  le  cas  contraire  et  faute  de  mieux  la  recherche  des  solutions  de  (l)  n'est  pour  le  moment 
entreprise,  k  de  rares  exceptions  prbs,  evoquees  en  II. 5.2,  qu'k  partir  des  mesures  stationnaires  classi- 
ques.  Cette  procedure  est  appliquee  dgalement  aux  mesures  obtenues  avec  des  balances  en  rotation  uniforme 
parallble  ou  non  k  l'ecoulement. 

Or,  les  installations  au  sol  disposent  a  l'heure  de  moyens  de  traitement  assez  puissants  pour  aborder  la 
recherche  des  solutions  de  (l). 

II.  5.1-  Methodes  lindaires 


Elle3  ne  s'appliquent  qu'a  la  determination  des  matrices  Cp^get  Cr.e  dont  les  ElEments  sont  des  constantes. 
II.  S.  1.1  -  Oscillations  libres 

Le  principe  du  traitement  est  rappelE  dans  le  cas  simple  des  oscillations  de  tangage  h  un  seul  degrE  de 
liberte  A©  ,  avec  rappel  dlastique. 

Dans  ce  cas,  compte  tenu  de  : 


A  0  =  Ao 


Les  transposdes  du  vecteur  d'etat  et  de  sa  ddrivde 
(5)  s'ecrivent  : 


AE  =  [Aa  Aq*  3 


•t 

par 

et 


A0  =  Aq  , 

rapport  au  temps  qui  figurent  au  premier  membre  de 

il*-  [Ad*  Aq*]  . 


de  m@me,  d'aprks  l'expression  de  M}~,dcrite  en  1.4.5  ,  les  elements  des  deux  matrices  Cp.|:  et  C||  j 
ddtermindes  par  les  mesures  faites  au  cours  de  ceo  oscillations  peuvent  s'dcrire  sous  la  forme  d'un  vecteur 
de  transposde  : 


[Cma  Crnq+Cmj] 


Au  second  membre  de  (5)  le  tenseur 
\/t  ;  en  outre  : 


Par  ailleurs 
avec  : 

Sx.p  rigidite  de  1' articulation 
constants,  de  l'en3emble  mecanique 
et  A©  sa  Vitesse. 


nt  et  la  matrice  k"'  se  rdduisent  respectivement  aux  scalaires  ly  et 

AGx=  0 


Elastique,  co$  et  Ct  pulsation  propre  et  amortisseaent  rEduit,  supposes 
constituE  par  la  maquette  et  1'articulation,  A Q  Elongation  du  mouvement 
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Dans  ces  coalitions  1 ’Equation  (r0  prend  la  forme  ie  liquation  lindairs  ’i>  second  ordre  suivante: 


(  6  ) 


[Ae  Ael] 


Cma 

lCmq+cmo 


=[-A&-2;.0)!A6-«u;,A»]  iy/|v’s 


et  sous  forme  sans  dimension,  exprimee  par  le  symbole  (  )* 


(  7  ) 


[Att  At)  ]  fCmo 


Cm  a 

Cmq+ Cn 


=  [-A0-2;,osA0  -<o8JA0]*pryJ 


La  solution  de  (7)  est  de  la  forme  : 


Ae=  A0g«st 


Lorsque  1* elongation  AO  a  la  forme  d'une  oscillation,  S  valeur  propre  complexe  de  liquation  (7) 
(7)  a  pour  expression  : 

S=  u£ (-;♦  i/T^5)  ...e  U»t *  Vijt/  V 


<oJ  et  C  sont  respectivement  la  pulsation  propre  et  l'amortissement  r^duits  de  l'oscillation  ; 
1'  elongation  A0  peut  @tre  eliminee  de  (7)  qui  devient  : 


(  8  ) 


[1  s] 


'ma 


=  [“S2-2;#t.j®*s“wS#2]  |iry2  av«c  iy*2=  ly/r 


j 

i 


La  procedure  experimental  consiste  : 

-  k  determiner  en  absence  d'ecoulement  :  S*  g  par  des  mesures  statiques  ou  en  oscillations  avec  une 
fausse  maquette  d'inertie  connue,  puis  oiget  s*  en  oscillations  avec  la  maquette  a  essayer  i 

-  les  reponses  en  oscillations  dans  l'ecoulement  aonnent  u>0  et  ^  •  d'ou  l'on  deduit  S  (grandeur 
complexe) . 

On  dispose  alors  des  parametres  necessaires  au  calcul  de  Cmtt  et  Cmq +  Cmj  d ' apres  (8). 

Pratiquement,  les  frequences  et  les  amortissement  reduits  sont  mesures  par  des  appareils  specialises 
(dampometer)  a  sorties  numeriques 

Le  premier  d'entre  eux,  a  notre  connaissance ,  a  ete  conqu  et  realise  avant  1954  [51]  par  OLSON  et 
ORLIK-RUCKEMANN. 

Les  moyens  d' essais  en  oscillations  libres  ont  ete  principelement  etudies  pour  la  mesure  des  faibles 
amortissements  de  tangage  des  missiles  e lances  en  supersonique  et  hypersonique.  Les  valeurs  de  C 
obtenues  aux  essais  sont  de  l1 ordre  de  : 

S<0,05. 

L'inter$t  de  cette  methode  reside  d'une  part  dans  la  simplicity  de  la  mesure,  limitee  au  releve  de  1' elon¬ 
gation  en  fonction  du  temps  et  d' autre  part  dans  la  rapidite  d' execution  des  essais.  En  effet,  le  montage 
presente  figure  All  permet  d'obtenir  plusieurs  series  d 'oscillations  au  cours  de  la  duree  de  l'eccu- 
lement  sain  d'une  soufflerie  intermittente. 

II.  5. 1.2-  Oscillations  entretenues  et  forcees 

Les  m^thodes  de  traitement  pratiquement  identiques  pour  les  deux  modes  d ' oscillations  ont  pour  obje^  les 
mesures  de  la  frequence,  des  amplitudes  et  phases  du  mouvement  et  des  torseurs  d' effort. 

Lorsque  le  mouvement  est  a  un  seul  degre  de  liberte,  sa  phase  est  prise  comme  origine  ;  mais  gen^ralement , 
le  mouvement  est  a  plusieurs  degree  de  liberte  et  1* origine  des  phases  est  alors, par  convention,  celle  du 
degre  de  liberty  d' amplitude  maximale. 

Les  procydures  experimentales  et  les  methodes  de  traitement  sont  classiques  ;  elles  sont  rappelees  dans 
les  communi cat ions  1  a  9  du  symposium  FLP-AGARD  d'Athbnes  (1978).  II  a  done  paru  suffisant  pour  fixer  les 
idees  de  presenter  seulement  ici  figure  All-6  deux  diagrammes  qui  resument  clairement 

les  methodes prat iquyes  par  la  D.F.V.L.R.  et  I’industrie  aerospatiale  de  la  R.F.A.  [42] 

-  La  figure  AIK6  montre  un  exemple  d'utilisation  des  multiplieurs  analogjques  (RCV)  a  trois  encrees 
lesquelles  sont  raccordees  simultanement  au  signal  a  mesurer  et  a  deux  signaux  sinusofdaux  purs  decales 
en  phase  de  90°  ;  les  resultats  des  multiplications  du  signal  a  mesurer  par  ces  deux  signaux,  disponibles 
sous  forme  analogique  aux  deux  sorties  de  l'appareil  sont  appliques  a  un  voltmetre  integrateur  a  3ortie 
numyrique  (iRVDM).  La  resolution  de  (5)  est  ensuite  faite  par  ordinateur. 

Cette  methode  est  utilise?  couramment  par  le  NAE  (voir  communication  9  au  symposium  city  ci-dessus),  et 
dans  les  souffleries  de  recherche  do  1'ONERA  a  Chalais-Meudon. 


I 


i9 


t*  figure  II. 1  tlrd*  da  la  rdf  Irene*  jtoJ  aoatr*  1m  val*ur*  d«*  co*fflol*nt«  d*  ac«*at 
d*  tajutuc*  *t  d*  lac*t  par  rapport  4  la  vitas**  angulair*  d*  lae*t,  ottam**  av*e  la  >alane*  du  KAI 
rsprdasntds  fig.AlI-< -2.  las  assures  oat  its  sffsc tud*s  dans  la  soufflsri*  tzaaasoniqu*  d*  la  (ASA  A 
Aass  (  1,83  ■  x  1,63  ■  )  A  Hach  jK  m  0,7  dan*  un  iatarvall*  4' incidence  d*  0*  A  40*. 

Oss  discontinuitds  iaportantaa  but  css  deux  co*fflci«ats  apparaissant  sntr*  16*  *t  20*| 
us  changeaent  d*  signs  est  obssrvd  sur  1*  coefficient  croisd. 


cmr  'cm4  cos  <*o  i.r 


it  a  n  >* 

*LFH*  (»CB) 

Dynamic  Pitching  Moment  Derivative  due  to  Yawing. 


i  it  a  x  hi 

HLM*  (Id) 

Damping-in-Yaw  Derivative.  M  -  0.70f  0  -  -5°. 


Mg.  II. 1 
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-  La  figure  A 11-6  donne  le  schema  d'un  ays t feme  de  traitement  numdrlaue  par  analyse  de  POURIE"  liaitde  fe 
la  frequence  fondamentale .  Lee  signauz  analogiques  dais  par  les  ddtecteurs  sont  conditionals  (Relais-unit), 
puis  aiguillds  par  un  systfeme  logique,  qui  les  branche  sequential  lament  sur  un  convertiaaeur  analogique- 
numdrique  et  les  transmet  fe  l'ordinateur.  Celul-ci  effectue  l'analyse  de  FOURIER  et  la  recherche  des 
solutions  de  l'dquation  (5)  correspondent  au  mouvement  imposd  fe  la  maquette . 

Cette  mdthode  plus  modems  que  la  prdcddente  eat  appliqude  fe  l'ONERA  dans  les  souffleries  industrielles  de 
Modane  et  du  Fauga.  Les  rdsultats  prdsentds  dans  la  communication  11  du  symposium  d'Athfenes  ont  dtd  obtenus 
de  cette  faqon. 

II.  5.  2  -  Mdthodes  non  lindaires 

11.5.2.1-  Les  questions  posdes  par  la  recherche  des  qualitds  de  vol  satisfaisantes  dans  les  domaines 
particuliers  parcourus  de  plus  en  plus  frdquemaent  par  les  adronefs  actuals ,  tels  que  rafales  atmoaphd- 
riques,  vol  de  transition  des  ADAC/ADAV,  manoeuvres  fe  grande  incidence  des  avions  de  combat  en  transso- 
nique,  etc...,  ne  peuvent  Stre  dtudides  par  les  mdthodes  lindaires  exposdes  au  paragraphs  prdcddent. 

Les  variations  de  R  ,  gdndralement  instationnaires  sont  souvent  dues  fe  des  ddcollements  et  recollements 
prdsentant  des  phdnomfenes  d'hystdrdsis.  Elies  ddpendent  non  seulement  des  valeurs  instantandes  de  E 
et  c  ,  mais  aussi  de  leur  dvolution  dans  le  temps.  En  consdquence,  les  mesures  de  R  ,  obtenues  en 
fhisant  evoluer  une  maquette  dans  ces  domaines  avec  des  valeurs  de  E  et  E  particuliferes  imposdes  a 
priori,  ne  donnent  que  des  rdponses  insuffisantes. 

Ces  considerations  expliquent  le  ddveloppement  actuel  des  installations  d'essais  au  sol  de  maquettes  en 
vol  libre,  mentionndes  au  paragraphs  11.3.4  ,qui  apportaient  ddjfe  depui3  de  nombreuses  anndes  une  impor- 
tante  contribution  dans  l'dtude  des  qualitds  de  vol  fe  basse  Vitesse  (vol  en  approche  et  vrille).  A  l'heure 
aotuelle  ces  installations  prdsentent  un  intdrdt  accru,  car  elles  permettent  d'dprouver,  dans  de  bonnes 
conditions  de  cotit  et  de  ddlai,  des  mdthodes  de  traitement  non  lindaire,  couramment  utilisdes  dans  les 
techniques  d'asservissement  [52],  mais  qui,  jusqu’A  ces  derni&res  anndes,  n'avaient  4t4  appliqu 4ea 
qu'exceptionnellement  en  a^rodynamique  exp^rimentale  ;  de  nombreuses  Etudes  sont  poursuivies  maintenant 
dans  cette  voie,  exposdes  en  grande  partie  dans  les  communications  n°s  19  &  32  de  la  reunion  AGARD 
d*  A  thanes  en  1978. 

En  effet  l’historique  de  la  mdcanique  des  fluides  montre  que  si  des  dtudes  sur  les  4coulements  instation¬ 
naires  non  d£coll£s  ont  4t4  entreprises  et  d^velopp^es  depuis  plus  de  cinquante  tins  [53]  et  [54]  ,  celles 
des  dcoulements  instationnaires  ddcolles  vient  fe  peine  d'etre  abordde.  Les  mdthodes  de  traitement  non 
lindaires  ddveloppdes  depuis  une  quinzaine  d' anndes  permettent  dfes  fe  prdsent  l'dtude  de  ce  problfeme  sur 
le  plan  expdrimental . 

n  5.2. 2  -  Eremple  de  traitement  d'un  phenomfene  d'hystdrdsis  adrodynamique . 

L'analyse  des  oscillations  de  tangage  autoentretenues  d’une  maquette  de  corps  de  rentrde  de  rdvolution  fe 
jupe  tronconique  en  dcoulement  supersonique  dlevd  ( =  4,5)  presentd  par  X.  VAUCHERET  en  1965  [55]  et 
rappelde  pages  10-8  de  la  rdference  [20]  est  fe  notre  connaissance  l'une  des  premiferes  applications  d'un 
traitement  non  lindaire  d' informations  fe  un  phdnomfene  de  ddcollement-recollement  d'un  dcoulement  adrody¬ 
namique  . 

La  maquette  dtait  disposde  sur  le  montage  reprdsentd  figure  AH  et  calee  fe  une  incidence  °  de 
7  degrds  (  &-  u  ) , 


En  transition  naturelle.  1' dcoulement  fe  cette  incidence  dtait  decolld  sur  toute  la  partie  cylindrique  de 
l'extrados  et  ne  prdsentait  aucun  ddcollement  sur  l'intrados. 

A  1' incidence  nulle  (  0  =  0°),  le  ddcollement  dtait  limitd  fe  la  moitid  arrifere  de  la  partie  cylindrique. 


La  maquette  libdrde  atteignait  aprfes  un  certain  nombre  d' oscillations  un  cycle  limite  de  1,5  degrd 
d'amplitude  environ.  Ce  cycle  limite  stable  se  retrouvait  quelle  que  soit  l'incidence  du  lficher,  y  compris 
1' incidence  nulle. 


amplitude 

(dagrda) 


Fig. II. 2 


V  Evolution  da  1'saplitud*  an  fonctlon  du  noabra  d’ oscillations  aat  reprdnentde  figure  11,2 
ainsl  qu*  calls  da  la  ataa  aaquetta  an  transition  ddclenchde  i  on  observe  qua  catta  Evolution 
aat  lindaire  sans  aucnna  taadanca  A  us  cycle  linlta. 


Fig.  11.3 


la  court*  schdaatlqua  da  Cm  en  fonctlon  da  a.reprdeentde  fig. II. 3  ,  trsduirs.it  cette  hypothAae. 

Le  point  figuiatlf,  coneiddrd  k  partlr  du  point  A  correspondent  A  1' Elongation  naxiaala  («=  0L=1,5*), 
parcourt  d' abort  le  aegnent  de  droite  AC,  dont  la  pent*  negative  (  Cn^cO)  indique  un  aoaent  da  rappel 
adrodynanique  stable  j  sals  la  droits  AC  ne  pasae  paa  par  1' origins ,  car  la  dlseyndtrie  de  l'dcouleaent, 
due  au  ddcollenent  iapoee  k  la  aaquatte  un  Cmo  ndgatif . 

Au  point  C  l'dcouleaent  eat  en  partis  recolld,  la  pent*  Cmo  de  CD  eet  plus  grande  en  valeur  abeolue 
qua  celle  de  AC  i  male  CD  ne  paaae  paa  par  l'origine  car  2orsque  la  naquette  re passe  A  l'lncidence  nulls, 
l'dcouleaent  n'a  pas  encore  totalanent  recolld. 

L’dcouleaent  ddcolle  de  nouveau  en  D,  la  droite  support  du  segnent  DE  a  la  seas  pants  qus  SC  par 
raison  da  ayadtrle. 

Da  alas  la  trajst  da  retour  SFBA  aat  syndtriqua  par  rapport  A  l'origine  da  l'allar  ACDB. 

La  surface  da  la  boucls  lialtde  par  la  contour  BCDF,  dgala  A  t 

~£Cmdo 

aat  positive  at  proportlonnalla  A  l'dnergls  correspondant  A  l'entxetlen  du  cycle  linlta.  la  choir  da  la 
fone  nulti lindaire  du  aoaent  de  rappel  adrodynanique  a  pernis  d'intdgrer  par  noroeaux  l'dquation  du  nou- 
veaent  dcrlta  sous  foraa  lindaire . 

11.5. 2.3-  Rennique  aur  l'intarprdtation  de  non  lindaritds  donnant  lieu  A  das  cycles  liaitaa. 

Ia  courba  representative  de  la  norne'du  torseur  das  forces  adrodynaaiques  [Cp|en  fonctlon  de  la  none* 
du  vectaur  d'dtat  IE*  Id'un  cycle  linlta  as  prdaente  gdndralenent  sous  la  fone  d’une  oenrba  fanda. 

Ia  caa  da  la  courba  fanda  sinplement  connaxa,  sane  point  double  ,  eat  discutd  id. 

Ia  variation  da  E<  dtant  pdriodiqua,  il  eat  possible  da  ddvelopper  an  add*  da  Fourier  |E*|  at  ICrI* 
si  on  as  linlta  A  la  frequence  fondaaentale ,  la  courba  raprdaentative  ICyq I ,  IE3  a  la  fora*  d'uns  ellipse 
seeblabl*  A  calls  qua  l'on  paut  obaarvsr  an  rdglae  lindaire, 

Le  aouvanant  dtant  auto-entretenu  ,  certains  auteurs  interprttant  laa  informations  d'oasaia,  recueilllee 
gdndralenant  aprda  flltraga,  an  tradulsant  l'dnergls  ndcessaira  A  l'entretlen  du  aouvanant  ooaaa  rdaultant 
d ' aaortisaaaanta  ndgatlfa  ou  mils,  c'ast  A  dir*  tala  qua  t 

Cmq+Cmd>0  ou  cnr  “  cn$  eo*«  *  0  . 

Urn  tails  interpretation  qul  parent  da  fairs  cofncidar  artificlallaaant  dee  rdeultata  expdriaentaux 
at  dsa  rdaultata  tbdoriquaa  da  ealculs  lindariads,  eat  aanifestauent  ineorraota,  oar  alia  naaqua  la  nature 
da  l'dcoulsaant  adrodynanique  rdel  dont  laa  affats  ralAvant  d'un  processus  non  lindaire  inatationnalre . 


0.5. 2.4  -  Remarque*  aur  l'dtat  actual  das  ndthodaa  da  trai tenant  daa  caa  da  vol  A  grand*  incidence  an 
dcoulenent  ddcolld. 

A  notre  cenaalasaaca,  il  n'exiate  pas  A  l'haura  actuelle  da  ndthoda  gdndrala  pernettaat  de  ddgnger 
d'una  adria  d'erpdrieaca*  ainplaa  ,  an  aouffleri*  ou  dans  un*  installation  au  sol  da  vol  libra,  daa 

infonations  aufflaantsa  pour  prddira  cvec  certitude  une  Evolution  quelconque  d’un  adronef  eh  rdgioe  ddcolld. 

Cy'l-  •• 


42 


In  effet  l«a  adthodee  d'dtude  de*  cm  de  vol  k  grande  Incidence  an  dcoulesent  ddcolld  m  prdsentent 
actuelleaent  seua  deux  Mpecta  i 

1)  lee  essaie  dene  lee  lnatallatlons  de  vol  llbre  qul  perasttent  de  elauler  avec  une  aaquette, 11 bre 

de  toute  contraiate  iaposde  par  une  suspension,  le  aouveaent  de  1 1 adronef  en  prdsence  d'un  toreeur  R 
fonctlon  lnstatlonnalre  de  E  et  de  eon  dvolution  daaa  le  teape.  eont  effectude  k  dee  noatarea  de  Mach  at 
de  Reynolds  tron  falblee.  * . . . 

2)  lee  eaMl'e  dare  lea  souffleriee  aux  noebree  de  Mach  de  l'adronef  rdel  et  4  dee  noabrea  de  Reynolds 
plus  1* portents  qua  ceux  de  1)  eont  conduits  d’aprte  lea  adthodee  llndairee  expoadea  en  II. 5.1.  Celles-cL 
ndceaaltent  la  prdaence  d'une  suspension  et  dans  certains  cm  d' organa*  soteurs,  qul  ispoeeirt  dee  contnalntee 
k  la  aaquette  ( Sx  #  O)  ,  s'oppoeant  k  la  slaulatlon  de  toute  dvolution  E(t)  de  l'adronef  rdel  ne  vdrlf lent 
pas  l'hypothbee  de  llndarltd. 

las  indications  qui  suivent  but  l'dtude  dee  perturbations  des  quell tde  de  vol  sigaaldM  en  I. 6.2.2 
illustrent  les  re Barques  prdcddentes  i 

-  1 ‘apparition  du  trssbleaent  (buffet  onset  )  est  no  tde  pendant  1m  assures  jtetlonnalres  claselques, 
fig.  I  .  6  ,  et  1m  spectres  de  puissance  des  pressions  sur  l'aile  eont  relevdM  au  cours  de  eon  ddvsloppe- 
aent  . 

-  le  doaaine  d' incidence  et  de  ddrapage  correspondent  au  balanceaent  d'ails,  k  1* inversion  de  la 
conaanie  de  gauchisseaent  et  au  ddpart  en  vrille  avec  parts  de  controle^eat  caractdrisd  par  les  rale ure 
des  deux  coefficients  ddduites  llndalreaent  des  assures  stationnalrea,  solt  : 

Cnpo=  Cnpco*°  (  Cnj  dynaaique  ) 

LCOP  =  Cnf  _  (  lateral  control  departure  parses  ter). 

<CngL.  »»«leetl«*»  da  lacat  S  dareulia  da  la  gauaarna  da  gauehiaaamanl  ) 

Lorsque  Cn^D  et  LCDP  e'annulent  et  deviennent  ndgatife,  11  y  a  une  forte  probability  pour  que 
l'adronef  se  trouve  en  eltuation  de  ddpart  en  vrille  avec  perte  de  controls. 

Dam  le  cae  particulier  du  HARRIGR  (3  5]citd  en  1. 62.2, l'dtude  de  eoufflerie  n'a  permis  de  ddfinir  qu'un 
domaine  limite  d'apparition  des  perturbations,  presente  fig.  It. 4  ,  ayant  sensiblement  le  ra§me  aspect  que 

celui  observe  en  vol. 


D'une  aanifcre  gdndrale,  11  apparai t  que  lee  donndes  recueilliee  par  les  adthodee  utillsdee  actuelle¬ 
aent  aa  foumiasent  que  das  inforaatiom  insuff lsantes  pour  aettre  en  oeuvre  de  faqon  satlafalsante  1m 
elaulateure  d'dtude  des  qualltde  de  vol  et  d'entralneaent  dee  pllotes. 


II  Mt  done  indispenMble  d'eseayer  de  coabler  cette  lacuna  de  l'expdriaentation  au  sol.  Pour  attelndre 
cat  objectif,  il  faudrait  pouvoir  rdalieer  un  aontage,  en  eoufflerie  transeonique  par  exeaple,  eonpu  pour 
alauler  avec  una  aaquette  le  aouveMnt  de  l'adronef  rdel. 

Un  tel  aontage  coaporterait  plus i sure  degrds  de  libertd  aaaervia  k  une  lol  de  aouveaent  caleulde  en 
team  rdel  d'aprte  lee  Msures  lnstantandee  dee  torseura  R  ou_  Ex  de  la  aaquette,  de  fa?on  i 


1)  k  annular  S*  dans  Is  cm  oil  les  tsmeurs  de  Basse  rdduite  de  l'adronef  et  de  Is  aaquette  ont  la 
sdaa  valsur, 

2)  4  lapossr ,  la  cm  dchdant,  les  valsurs  ds  S*  correspond  ant  4  la  difference  entre  ces  deux  tense  ure. 

2 

la  tense ur  le  imh  rdduite  a  pour  expression  i  _ 


avec; 


ns 


-Ui  o 

o  ;'i/mr2 
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CONCLUSION 


Un  tour  d'horison  but  l'Etat  actual  des  Doyens  au  sol,  destines  &  la  determination  des  parang tree  aEro- 
dynamiques  de  stability  des  aEronefs  modemes,  a  EtE  effectuE  aprbs  un  bref  rappel  des  bases  thEoriques 
de  la  dynamique  du  vol,  ezposdes  de  faqon  concise  avec  des  synboles  matriciels.  Les  ElEments  des 
matrices  sont  exprimes  conformEment  aux  recommendations  des  normea  ISO  aur  la  mEcanique  du  vol. 

L'extension  de  telles  normes  Internationales  aux  matrices  elles— mftmes  faciliterait  considErablement  la 
redaction  et  la  lecture  des  publications  de  plus  en  plus  nombreuses  sur  ce  sujet. 

II  ressort  de  cet  examen,  effectuE  principalement  d'aprEs  les  communications  prEsentEes  aux  trois  reu¬ 
nions  de  l'AGARD  organises  par  le  FMP  a  Valloire  (Prance),  A  Rhode-Saint  GenAse  (Belgique)  en  1975  et 
par  le  FDP  A  AthAnes  (GrAce)  en  1978,  que  les  centres  d'essais  aerodynamiques  sont  actuellement  EquipEs 
pour  determiner  dans  des  conditions  satisfaisantes  les  paramAtres  de  stabilite  tant  qu'ils  restent 
justiciables  des  mEthodes  linEaires. 

Or  ces  mEthodes  ne  sont  gEnEralement  pas  applicables  dans  les  cas  de  vol  comportant  des  dEcollements 
instationnaires  brutaux  de  I'ecoulement  sur  l'aile.  Dans  ces  cas  particulars,  les  phEnomAnes  aErody- 
namiques  &  un  instant  donnd,  ddpendent  non  seulement  des  valeurs  instantandes  du  vecteur  d'dtat,  mais 
aussi  de  1' evolution  antdrieure  de  ce  vecteur. 

Seules  les  installations  de  vol  libre  permettent  d'obtenir  avec  une  maquette  la  similitude  temporelle 
de  cette  dvolution.  Mais  le  domaine  de  Mach  et  de  Reynolds  couvert  dans  les  installations  existantes 
reste  trAs  insuffisant. 

II  est  done  ndeessaire  de  tenter  de  combler  cette  lacune.  Les  progrAs  rdeents,  d'une  part  des  techniques 
d'asservissement  des  systemes  mecaniques  et  d'autre  part  des  moyens  de  calcul  permettent  d'envisager 
la  conception  de  nouvelles  methodes  d'essais  dans  des  souffleries  modemes  &  grand  Reynolds  en  cours 
d' dlaboration. 

Le  montage  utilisd  pourrait  disposer  de  plusieurs  degrds  de  libertd  asservis  A  une  loi  de  mouvement 
calculde  d'apres  les  mesures  instantandes  faites  sur  la  maquette.  II  serait  alors  possible  de  rEaliser 
avec  une  bonne  approximation  le  comportement  de  1'aEronef  au  cours  des  di verses  Evolutions  prEvues 
dans  le  domaine  non  linEaire. 
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ANNEXE  All 


-  RAPPEL  DE  NOTIONS  THEORIQUES  SUR  U  DYNAMIQUE  DU  SOLIDE  INDEFORMABLE  EN  KOUVEMBNT  DANS  L' ATMOSPHERE  - 
t  -  Repdres  - 

Le  mouvement  par  rapport  d  la  terre  i'un  solide  indgformable ,  dans  l’atmosphdre  supposde  immobile  eat 
calculi  en  utilisant  deux  repdres  orthonormds  de  sens  direct  0  xJLZo  et  GxyZ  definis  de  la  fa? on 
suivante  : 


Ox.y.2.  dans  le  referential  lid  d  la  terre  - 
Le  point  0  eat  fixe  sur  la  terre  supposde  plane  et  oonfondue  aveo  le  plan  . 

La  verticals  de  0  est  le  support  du  veoteur  unitaire  Zg  orient d  vers  le  baa. 

Les  vecteurs  unitaires  X0et  y0  forment  avec  Zo  un  triddre  de  sens  direct. 

GxyZ  dans  le  rdfdrentiel  du  solide  - 

Son  origins  est  usuellement  le  centre  de  gravitd  (ou  centre  d'inertie)  G  du  solide.  Aucune  condition 
n'est  imposde  a  priori  dans  le  choix  des  deux  directions  qui  fixent  la  position  du  repdre  dans  le  solide, 
mais  s'il  possdde  un  axe  de  symetrie  cet  axe  est  pris  pour  support  du  vecteur  unitaire.  x  qui  est  orientd 
de  l'arridre  vers  l'avant  ;  le  point  avant  est  choisi  conventionnellement.  Si  le  solide  possdde  un  plan  de 
symdtrie,  les  vecteuxs unitaires  X  et  Z  sont  disposds  dans  ce  plan  de  fa?on  que  X  soit  orientd  de 
l'arridre  vers  l'avant  et  Z  oriente  vers  le  bas.  Les  points  avant  et  bas  sont  choisis  conventionnellement. 

2  -  Variables  d'dtat  - 

Les  variables  d'dtat  inddpendantes  qui  caractdrisent  le  mouvement  du  solide  dans  l'atmosphere  a  une  altitude 
et  un  instant  donnd,  sont  representds  par  le  vecteur  d  six  composantes  appeld  vecteur  d'dtat  (ou  torseur 
distributeur  des  vitesses1  )  .  Ce  vecteur  considdrd  dans  le  repdre  GxyZ  est  reprdsentd  par: 

Ex  s  V  v  11  qui  s'derit  sous  forme  transposde  Ejt  —  [  u  v  w  p  q  rj» 

avec  : 

V  vitesse  de  G  par  rapport  d  0  , 

ft  vitesse  angulaire  du  solide  par  rapport  d  0  x  y  Z  . 

3  ~  Equations  du  mouvement  - 

Les  six  dquations  qui  expriment  l'dquillbre  des  forces  d'inertie  et  des  forces  extdrieures,  ainsi 
que  l'dquillbre  des  moments  correspondents  par  rapport  au  centre  d'inertie  du  solide,  sont  reprdsentdes 
de  faqon  concise  dans  le  repdre  lid  au  solide  par  une  seule  dquation  matricielle: 

(D  «7»Eex  , 

avec  au  premier  membre  : 


O  =  Dm  E» 


D  est  une  ma trice  carrde  d'ordre  6  de  la  forme  : 


D  = 


d  ;  0 
:  «* 


dans  laquelle  : 

d  xatrice  carrde  d'ordre  3  ,  d'opdrateurs  ayant  pour  expression 


d=  Jfl  +  fiA  = 


1 

0 

°1 

0  -r  q 

0 

1 

0 

+ 

r  0*p 

0 

0 

1 J 

-9  p  o 

•  ma trice  carrde  d'ordre  3  dont  les  dldments  sont  mils, 

m  ma trice  carrde  d'ordre  6  (  tenseur  de  masse  )  de  la  forme  : 


•  0 


0  :  In 
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m  est  la  msae  du  solide  et  In  le  tenseur  d*inertie,  <SJest*4«dlre  : 


1 

0 

0 

lx 

-Ixy 

-Ixz 

0 

1 

0 

et 

In  ■■ 

-Ixy 

ly 

-iyx 

0 

0 

1 

>% 

X 

1 

-lyz 

l2 

Ex  veoteur  d'etat  ddfini  en  2  ci-dessus  , 
F.x  torseur  dea  forces  ext4rieures, 


Ce  node  d'expression  a  fetfe  prfeser.tfe  il  y  a  plus  de  vingt  ana  par  l'Ingdnieur  Gfendral  E.  BILLION  1  et 
repris  3ous  une  forme  plus  gdnferale  par  la  Professeur  P.  GERMAIN1 

Le  dtfvfeloppement  de  1' Equation  (l)  conduit  aux  six  Equations  suivantes  : 

Equations  des  quantltds  de  mouvement  : 


n>(.lH+qw  -  r  v  )  -  Xe 

<2)  m(il  +  ru  _  pw)  -  Ye 

m  (i"  +  pv  _  qu)  =  Ze 
at 


+  ( •*  -  ly )  qr  +  Iyr  ( r2  -  q2 )  +  C  Ixy  r  -  1**  Q  i  P  =  Le 


+  (1 »  —  I*  )  rp  +•  iXz  ( p2-  r2  )  +  (lyz  P-  IXy  f)q  =  We 


+  (ly  _  lx)pq  +  l*y  (q2-  p2  )  +  (lxzq  -  lyz  p)  r  =  N0 


i]  BILLION  E. 

1  —  Probifemes  modemes  de  dynamique  du  vol  (1956) 

2  -  Dynamique  du  vol  (1958-1959) 

CENTRE  D' ETUDES  SUPERIEURES  DE  MECANIQUE,  PARIS 

j]  GERMAIN  P. 

Mdcanique  des  milieux  continue 

Tome  1  Thdorie  Gdndrale,  chapitre  I,  paragraphs  1.2,  page  12 
MASSON  et  Cie,  PARIS  (1973) 


Equations  du  moment  clndtlque  : 


(3) 


lx 

d£ 

-  Ixy 

-  lx* 

d  r 

dt 

v  dt 

dt 

ly 

dq 

-  lyx  iL 

-  Ixy 

22 

dt 

d  t 

dt 

lt 

dr 

-  *xx  lE 

-  1  yx 

dq 

dt 

"  dl 

dt 

Al-3 


ANNEXE  AI.2 

CHANGEMENT  DE  REFERENTIEL 


Les  repferes  utilises  couramment  en  dynamique  du  vol  sont  orthonormds  et  ddfinis  dans  l’espaoe  ouclidien 
&  trois  dimensions.  Le  centre  d'inertie  G  de  l'adronef  est  pris  comme  origine  des  reperes  considdrds 
dans  la  prdsente  annexe. 

Les  changements  de  rdfdrentiela  qui  font  l'objet  des  trois  premiers  paragraphes  A  I. 2.1.,  A  1.2.2. , 

A  1.2.3.  sont  effectuds  dans  des  repferes  de  mSme  origine  sur  l'adronef,  definie  par  un  tenseur  de  masse 
donne.  Ces  repdres  se  ddduisent  done  l'un  de  1' autre  par  une  rotation  autour  de  G  .  alors  que  les 
reperes  considdrds  dans  le  dernier  paragraphe  A  1.2.4.  se  ddduisent  l'un  de  l'autre  par  une  translation, 
leurs  vecteurs  unitaires  restent  identiques. 

Dans  oe  qui  suit  une  matrice  de  changement  de  repere  par  rotation  est  reprdsentde  dans  l'espaoe  a  trois 
dimensions  par  la  lettre  minuscule  pu  ,  Norite  en  caractere  gras.  L'indioe  infdrieur  caractdrise  l'opd- 
ration  rdalisde  par  la  matrice  de  transfert.  La  matrice  rdalisant  la  m£me  operation  entre  torseurs  a  6 
dimensions  est  symboliade  par  la  lettre  majuscule  l*(  ,  affectee  du  mdme  indice  ;  elle  a  pour 

expression  : 


*?>  = 


?»i  0 

•IV 


Al.2.1-  Passage  du  repere  normal  terrestre  6  X0y020  au  repere  li4  a  l'adronef  Gxv  2 

Ce  passage  peut  Stre  effectue  par  les  3  rotations  successives  ,  ©  ,  <D  ,  repr^sentees  dans  le 

sens  positif  figure  page 

-  La  rotation  autour  de  l'axe  vertical  Z0  amkne  le  repere  G  Xey0Zo  en  Gxhyh20  ;  leurs 

vecteurs  unitaires  sont  relies  par  : 


X 

*o 

cosy 

sinV 

0 

Vh 

Vo 

avec 

Pr  = 

-sin,^, 

cosT 

0 

*o 

.  — 

zo 

0 

0 

1 

© 

autour  de 

l'axe  yh 

amhne  le  repere 

6  *.,*,*. 

en  G  x 

i  leurs 

re lids 

par  : 

X 

*h 

cose 

0 

-s'ln© 

vh 

& 

ii 

Vh 

avec 

Pa  = 

0 

1 

0 

Zq 

sin0 

0 

cos© 

Z  est  dans  le  plan  de  symetrie  x  z  . 

-  La  rotation  <D  autour  de  l'axe  x  araene  le  repere  G  X  yh  Z,  en  G  xyz  ;  leurs  vecteurs  uni¬ 
taires  sont  relies  par  : 


X 

X 

’  1 

0  o 

V 

z 

-  - 

* 

If 

vh 

*s 

avec  P(j,  = 

0 

0 

cos$  sin4> 

-sin4>  cos4> 

Les  vecteurs  unitaires  des  repferes  GxyZ  et  Gx0JbZ0 
passage  suivante  : 


sont  par  consequent  relies  par  la  matrice  de 


X 

*o 

Y 

ll 

Vo 

z 

zo 

*0  “  Py 


AM 


1 


dont  1' expression  ddveloppde  s'dcrit  : 


po 


cosOcosT 

-cositsin’V  +  sin<6sin0cosY 
ein'lisin  'f  +  cos't>sin©cosH' 


cos8sinY 

oostf'cos'V  +  s  in<J>s  inesin'Y 
-siwbcosH?  +  nos4>sin©sin'i' 


-sin© 

sin4>cos© 

cosicos© 


et  la  relation  inverse  : 


Po,  =  W*>=P° 


A I.  2.2-  Passage  du  repere  lie  a  l'aeronef  G  xy  z  au  repere  a&rodynamique  G  Xaya  Z» 

Ce  passage  est  effectue  par  les  deux  rotations  successives  01  et  P  represent^es  dans  le  sens  positif 
sur  la  figure  i  page  6 

-  La  rotation  Ql  autour  de  Oy  am&neGxyz  en  6  x,y  Z,  ;  leurs  vecteurs  unitaires  sont  lies  par  : 


■*s' 

X 

"cos  a 

0 

-sina  ’ 

Y 

II 

y 

avec 

pa  = 

0 

1 

0 

2a 

Z 

_sira 

0 

coca  . 

x,  et  Z.  sont  dans  le  plan  X  Z  • 

-  Puis  la  rotation  ()  autour  de  Gz.  am&ne  Gx,y  Z,  en  G  xtytZ#  ;  leurs  vecteurs  unitaires  sont 
lies  par  : 


Xa 

Xs 

‘cos  p  -sinp  0 

Va 

V 

avec  = 

sin p  cosP  0 

0  0  1. 

.  2  a 

za 

Les  vecteurs  unitaires  des  reperes  G  x  y  Z  et  G  xa  y, Z,  sont  par  consequent  relics  par  la  matrice  de 
passage  suivante  : 


x." 

x" 

r  cosacoep 

S3  n  p 

si  nacos  p- 

V* 

V 

ii 

V 

avec 

H'e~  Pp Pa=  -cosasinp 

cos  p 

-slnasi n fi 

z 

L-slna 

0 

coco. 

et  la  relation  inverse  : 


X 

X. 

'cosocos  P 

-cosasinP 

-sina ‘ 

Y 

z 

.  - 

11 

f 

Y. 

z. 

avec  n  =  p  p.  =  p~'  s 

-a  -P  *36 

sin  P 

.sinacoeP 

cos  P 

-sinosin  P 

0 

coca. 

Al.2.3-  Relations  entre  les  elements  d'une  matrice  exprimee  dans  deux  reperes  deduits  l’un  de  1* autre 
par  une  rotation 

Soient  et  les  deux  reperes  de  vecteurs  unitaires  respectifs  X  ^  y1  Z ,  et  Vt  *1  $ 

P  la  matrice  de  passage  de  a  . 

Ces  vecteurs  sont  relies  par  : 


V2 

*  P 

Vi 

*2 

*1, 

A I -5 


Soient  v,  et  Wj  Xes  expressions  respectives  d'un  veoteur  quelconque  dans  et  ,  de  la  forme 

■*,'  *2 

v^m,  y,  ®t  v?  rn  m,  y, 

.*'J 

m,  et  mutant  des  matrices  carries  d'ordre  3  j  d'aprta  (l): 

v2=  pv,  , 

ou  bien  , 

M  f*'l 

(  2  )  m,  y,  •  pm,  y,  , 


r  *i 

Vi 

=  p-’  Vi 

*1 

*2 

et  conpte  +enu  de  : 


on  retrouve  la  relation  bien  connue? 

(  3  )  m2  =  p  m,p-'  , 

Si  les  deux  repferes  sont  orthonormds,  la  relation  (3)  deviant  s 

(  4  )  m2  =  p  m,  pi  . 

AI. 2.4-  Expressions  du  vecteur  d'etat  et  des  derivees  aerodvnaaigues  de  stabilite  en  fonction  du 
centrage  de  l'adronef 

Les  accroissements  du  vecteur  d'etat  et  les  matrices  de  coefficients  a^rodynamiques  de  stability  a  deux 
centrages  de  l'a^ronef  definis  par  les  points  P  et  G  »  sont  lies  par  les  relations  etablies  ci- 
dessous. 

Les  deux  reperes  de  l'aeronef  G°xyz  et  G  xy  z  ont  m8mes  vecteurs  unitaires  et  sont  ddduits  l'un  de 
1' autre  par  la  translation  G^G. 

Les  elements  de  la  vitesse  angulaire  sont  les  oSmes  dans  les  deux  reperes,  m is  ceux  de  la  Vitesse  sont 
diff brents  si  la  vitesse  angulaire  est  diff Create  de  z6ro,  soit  n  *•. 

-  A  un  petit  ecart  de  l’aeronef  de  sa  position  de  regime  correspond  un  accroissement  du  vecteur  d'etat 
qui  a  pour  expressions  : 


AeS  a  AV°UAQ 

AE,  «  AV°-GS3'A  AQ  S  Av°  j(1-<Fg)AQ 


dans  G°  x  y  z  • 
dans  G  X  y  Z  . 


II  en  results  la  relation  suivante  : 

(l)  AE«»riAEx  n*'AEx  a  AE?  . 


n  « 


n*.  b*aib  , 

0  s  I 


II  est  obtenu  en  changeant  le  sjgne  du  vecteur  b  ddfini  ci-dessous. 

X  ,  y  ,  Z  etant  les  ooordonnees  de  G  dans  G*X/Z  ,  b  a  pour  expression  transpose 


*  "B  =  ~  [ X  y  z]  d'oti 


0  Z  -y 

bA  =  -z  0  x  . 

y  -x  0 


-  Au  m§me  ecart  de  l'adronef  correspond  un  accroissement  du  torseur  des  coefficients  a^rodyr imiques  sans 
dimension,  d' expressions  transposes  : 


ACr  *  [Acs  ACy  AC}  ACJ  AC*,  ACR] 


dans  G**y  Z  • 
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ACr=[ACx  ACy  ACZ  \Cf  ACm  ACn]  dans  Gxyz 

avec  : 


[AC, 

'ac;' 

'AC;  ' 

AC| 

-1 

0 

0 

Ac;' 

\Cy 

s 

ACy 

ACm 

= 

A  Cm 

•  *  * 
+■  b/\ 

0 

1 

0 

ACy 

i  r 

L-'^z 

AC| 

_ACn 

ACn 

0 

0 

-1 

AC| 

Les  relations  precedentes  s ' expriment  de  fagon  condensee  par  : 

ACr *  8 \Cr 

avec  : 


B*= 

l  i  0 

: 

— *1 - 

i 

K*=  tf  A 

'  -1 

0 

0  o  ' 
1  0 

K*|  1 

0 

0  -1 

Les  transport  des  moments  de  G°  en  Q  impose  de  changer  le  signe  des  composantes  de  forces  suivant 
X  et  2  ,  car  ces  composantes  exprimees  conf ormement  aux  conventions  des  normes  sont  negatives 
lorsqu’elles  sont  orientees  dans  le  sens  positif  de  ces  axes. 

Les  changements  de  signe  sont  effectues  au  moyen  de  la  matrice  diagonale  ecrite  a  droite  de  b/\  . 


cr.e.AE* 

cr,e„  AE, 


B 


c8.es  Aes* 
,c8.t5  ae$* 


1' indice  superieur  ( )*  represente  une  expression  rendue  sans  dimension. 

Compte  tenu  de  l'equation  (l)  ci-dessus,  les  matrices  des  coefficients  de  stability  sont  reliees  par  : 


B* 


C8,E5 

c8,fe» 


rf 


La  relation  (2)  ci-dessus  s'applique  a  tout  repere  ddduit  de  G'xyZ  par  une  translation,  quelles  que 
soient  les  valeurs  des  angles  U  et  P 

II  parait  cependant  utile  de  l'illustrer  par  des  exemples  simples  qui  se  presentent  dans  la  pratique.  Le 
plus  courant  est  celui  d'un  changement  de  centrage  suivant  l'axe  x  . 

a)  relations  entre  les  coefficients  de  stability  longitudinale  en  fonction  du  centrage  suivant  l’axe  de 
"stability"  x,  . 

Les  matrices  bA  et  K*  ont  pour  expressions  : 


b*A-l 


r~s 

\y  x. 

*'  -X.  0 


K*  = 


'  / 
V* 


A1 


d’ou  les  deux  matricv  \  expliciteeo  dans  (2) 


J  i  0 

w  am  mm-mmmm  u  ••  m  • 

I - V  I 

!•/«:!  i 

l**'*»?  o ' 


1  I  X* 

. f . 

0  :  1 


sorte  que  la  matrice  des  coefficients  de  stability  a  pour  expression  en  fonction  de  xj 


CRL,EU  1 


CRL  est  reprdsentdr par  une  relation  de  m$rac  forme. 

b)  Relations  entre  les  coefficients  de  stability  transversale  en  fonction  du  centrage  suivant  l'axe  x 
du  rep&re  G  x  y  z  . 

11  soffit  da^s^  ce  cas  dc  remplacer  xa  par  x  pour  obtenir  les  expressions  des  matrices  b  ,  K  , 
B  et  n*-  et  le  calcul  donne  la  relation  : 


2CX 

0 

C»« 

0 

C*n+  x.Cxa 

0 

0 

0 

0 

0 

2CZ 

0 

Czo 

0 

CJq  +  x?Cza 

0 

0 

0 

0 

0 

0 

0 

Cmo-x^CZa 

0 

Cmq  +  **(  Cmo  ~  C  zq )  “ 

.0 

0 

0 

0 

0 

CRT,ET  ss 


0  0 

0  % 
0  0 


0  0 

0  C$r-«*CyJ 
0  0 

0  C|r  -  X*C,p 

0  0 


Cnd-x’Cyf  0  CRp-x^Cyp  o  Cnr-X*(C^- C)r)  ♦  X**C,r 


ANNEXE  AI.3 

EXPPESSIONS  DE  fl 

•  •  • 

Al.3.1  -  Expression  dc  *2  on  fonction  de  O  ,  0  V  . 

•  «  .  .  % 

Le3  vi tenses  angulaires  ^  ,  0  ,  peuvent  6tre  considert^es  comme  des  vecteurs  O-,  0X, 

*PX  t  sur  In  base  Gxy  z  . 

-  4^  a  pour  support  l'axe  x  ,  d’ou  son  expression  tranopos^e  : 

[<t>  o  o] 

0  a  pour  support  l'axe  yh  f  d'ou  son  expression  : 


et  s a  trano ponce 


0X  b  0  [  0  coc<J>  _nin$J 

-  a  pour  support  l'axe  t9  t  d'ou  son  expression  :  _ 


H*  ■  p  o 

**  O  , 
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<  t  sa  transposee  : 

*[  -sin®  sinOcoc®  eos^cos®] 

Le  vecteur  n  est  egal  a  la  somme  des  trois  vecteurs  ecrits  ci-dessus  : 

£1  •  <i>x  +  0*+'H* 

qui  s'ecrit  sous  forme  developpee  : 


p 

[<j> 

’l 

0 

-sin® 

n= 

q 

s 

p<i|e 

avtc 

PnS 

0 

cos«i> 

6in«I>cosG 

r 

[y 

0 

-sin«l> 

cosmos® 

et  la  relation  inverse  : 


p 

t 

sin<l>tg®  cos«l>tg® 

® 

= 

pn 

q 

•vac 

P«1  = 

0 

cosfl>  -sin4> 

y 

r 

0 

sin<J>/cos®  cos<l>/cos® 

A1 . 3 . 2  -  Expression  de  LI  en  fonction  de  P,  >  o  .  0  en  regime  de  vol  rectiligne  . 

Un  regime  de  vol  rectiligne  eat  caracterise  par  une  direction  fixe  par  rapport  a  0  Xojfe  lo  de  l'axe  X, 
du  repere  aerodynamique  G  x, 

Cette  condition  unique  a 'exprime  par  : 

Xa constant 

Xa  est  le  vecteur  unitaire  de  l'axe  longitudinal  G  X, . 

Seule  la  direction  de  la  Vitesse  de  G  reste  fixe,  son  module  peut  varier  de  fapon  quelconque  en 
fonction  du  temps  ;  il  en  est  de  m@me  des  composantes  de  fi  .On  se  propose  de  trouver  l'expression  de 
en  fonction  des  vitesses  angulaires  a  ,  fi  et  pt 

-  Soit  6*  le  vecteur  vitesse  angulaire  4  dans  le  repere  G  xy  2  i  d  est  porte  par  l'axe  transversal 

Gy  ,  sa  transposee  a  pour  expression  : 

6x=  [0  i  0] 

•  • 

P*  porte  par  x#  ,  f  est  portee  par  Za  ,  raais  orientee  en  sens  inverse  ;  en  effet  a  P>0 
correspond  une  valeur  r,  <  0  . 

II  resulte  des  deux  alineas  ci-dessus  que  : 

Q  m  +  paa 

Corapte  tenu  de  l'expression  de  la  raatrice  p  donn^e  en  annexe  A  1.2.2.  H  a  pour  expression  : 


P. 

0 

-p 


Y 

'  COS  OCOB  fi 

0 

slna' 

A 

•v.c  %  = 

slop 

1 

0 

d. 

slna  cosp 

0 

-COSO. 

et  la  relation  inverse  : 


Y 

a 

a 

a 

ii 

avtc 

p  = 

'  cosa/coB0 
-cosa  tgP 

0 

1 

sina/ cosf- 
-sma  tgP 

Y 

R 

R 

sinacosP 

0 

-COSO 
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ANNEXE  AI.4 

REUTION  ENTRE  L 1 ACCROISSEMENT  DU  VECTEUR  D'ETAT  AEx  DIKENSIONNEL  ET  SON  EXPRESSION 
SANS  DIMENSION  AE*. 

Soient  V  et  Va  ,  les  expressions  du  veoteur  vitesse  de  6  ,  respectivement  sur  la  base  (  x  ,  y  ,  X  ) 
et  sur  la  base  de  "stability"  (  X,  ,  y,  ,  X,  )  qui  ont  pour  traneposdes  : 

V  =  [u  v  w  ]  =  [Vcosacosp  ysinp  Velnacoap]  Va  *  fv  0  0] 

Leurs  accroissements  relatifs  ; 


4Y  -  r  iu  Aw  ~i 

V  «-  V  V  V  J 

sont  relies  par  : 


#t  AVj  =  [y  AP  Aa] 


<  av: 


avec 


rcosocoB p 

bJ  nP 


p;  - 


LsinacosP 


coca  BinP 
cob  p 
-slnasin  p 


-einacos  p 

0 

cosa  cob  P 


et  la  relation  inverse  : 


Yv*  = 


avec 


B  = 


COS  OCOB  P 
-COBOSiDp 

-Blna/cosp 


slop  einacos P" 
coep  -slnasinP 
0  cosa/cosp. 


Les  deux  formes  du  veoteur  d'etat  AEX  (dimensionnel)  et  AE*  (sans  dimension)  sont  ddfinies  d'aprfes 


AE«-  AVUAQ 


AE*  =  AVaUAft 


La  relation  entre  A  E*et  AE  X  ot  leurs  derlvdes  par  rapport  au  temps  peut  s ' ecrire  sous  la  forme  : 


AE* 

AE* 


=  Px 


V  AE»I 

.V5 


et 


AEX| 

AE 


=  Px 


V  AE* 


i, 

I  I 


La  ma trice  ct  non  inverse  ont  pour  expressions 


Px 


f  1  f . 

r 

r- . •.  ' 

*  j  j  •  j 

•t  R ?m 

S’ 

!  0  i 
•  « 

•  « 

r . \*\  : 

i  ®  \  * 

Mt  * 

i  0 

1. 

Z  V 

iZ . 

Ll .  «J 

J 

V  est  le  soalaire  module  du  vecteur  V  vitesse  de  6 
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ANNEXE  AI5 


REPRESENTATION  DU  TORSEUR  IES  FORCES  D ' INERTIE  D’ENTRAINEMENT  ET  DE  L’ ACCELERATION  CORRESPONDANTE 

PAR  UNE  FORME  LINEAIRE  EN 

Le  torseur  A  JJj  qui  figure  dans  I' Equation  (2)  represents  le  produit  de  la  aatrioe  antisymdtrique  AD 
par  le  torseur  mEg  ,  soit  : 


Ac^)«  AD  m  Ex  i 


dans  ce  produit 


'aha 

0 

0 

AftA 

AD 


La  propridtd  du  produit  vectoriel  permet  d1 eorire  : 


•t 


m  E, 


0 


m 

----+ - H 

0  I  In 


A  Jfi  =  - 


et  apr&s  1‘ introduction  des  grandeurs  sans  dimension  : 


mVA  j 

0 

A«' 

0 

mV  A 

0 

0  ! 

InftA 

AH 

0 

InA 

An 

Hj 

*xP*"*xy1*'***  r* 

V*  =  [  cosocobP  ein  (J  GinacoeP]  ! 

Infl  „ 
m/V 

u* 

l“!j 

1 

"S3* 

-Ixy|?+  Iy  q*-Iyz  r* 
-Ixyp*-Iy2q*+  Iir* 

on  obtient  1* expression  qui  figure  au  second  membre  de  (3)  : 


m^ka, 


0 


0  !h*a 


ae1 


dans  la quelle  mvy  l  est  un  scalaire,  k  est  explicitde  au  paragraphe  1.3. 1.1.  et  ! 


e,  = 


-1  n; . ! 

W\  0  i 

0 

-sinacos  P 

slnP 

■  0 

-H*  Hj' 

i  \  1\  : 

11 

*< 

> 

sinacosp 

0 

-COCO  COE p 

.  H*A  = 

0  -hS 

i  0  \l\i 
(. . .A  1 

-sin  P 

0O80CO6 p 

0  j 

-Hy 

0 

Le  rflle  de  0j  est  de  faire  ooncorder  les  signes  des  termes  d' effort  de  -A^  avec  ceux  des  termes 
correspondents  de  AR  imposes  par  les  normes  dans  k  (  1.3. 1.1.).  La  presence  de  8i  permet  d^dviter 
d' dventuelles  erreurs  de  signe  si  l'on  s'dtait  contents  de  changer  les  signes  correspondents  de  V *A 


A  ^ 


ft  peut  fitre  repr£sent£  dans  1* equation  (3)  par  une  forme  lin^aire  sans  dimension  analogue  &  celle 


adoptee  pour  Ar  en  remar quant  que 

de  sorte  que  l'on  peut  poser  : 

AJft-  afnjsrAe* 


|sv*k  -  *  v* 


•VtC 


0  I  V*A 

wl&K 


Al-I ) 


peut  6tre  d^oompos4e  d'aprfes  : 


-S*L  =  n 


|  0 

slnocosp 

o' 

0 

1° 

0 

0 

! 0 

-COSOCOBp 

0 

e 

j 

0 

, 

£?= 


0 


L.  • 

sin  p  ■ 


Slop 


^rT=u 


0  0 

0 

sinacosP  0  -cob 

0  0 

r - t 

0 

!  g  t_' 

My 

0  0  0 

0 

Hz  0  -Ht 

. 

0  0  0 

r  \ 

O 

0 

o  -hJ  o 

0 

0  0  0 

o 

*XX 

O 

Ces  expressions  se  simplifient  d'une  manifere  evidente  dans  les  cas  oil  a  ,  P  et  ft  3ont  nuls . 

L1  acceleration  d'entralnement  est  calculde  en  multipliant  A  ^  pour  la  matrice  inverse  de  la  masse 
apparente,  soit  : 

(m  +  maj“'A  =  (  m  +  m,]"'  Jpy  AE* 

On  retrouve  dans  cette  expression  la  forme  lindaire 


ANNEXE  All. 1.1 


ROTATION  UNIFORME 


BASSE  VITESSE  -  SOCIETE  MACCHI 


1  nit  support  de  la  maquette, 

2  rail  de  guidage, 

3  bras  toumants  at  chariots  mobiles  but  le 
rail. 

4  ferruxe  d'attache  de  la  maquette  et  de 
rbglage  de  l'assiette  longitudinale  et  de 
1' angle  de  roulis. 


5  dynamombtre  b  cinq  composante3  libes  aux 
axes  tournants  de  la  balance. 

6  pro'.ongateur  interchangeable, 

7  ferrure  d'extrbmitb  du  mSt  c8tb  rail, 
rbglage  de  l'azimut, 

8  dynamombtre  mesurant  1’ effort  parallble  a 
l'axe  du  mfit. 


Installbe  dans  una  sorfflerie  &  veins  circulaire  ouverte,  dans  laquelle  la  Vitesse  de  l’bcoule- 
nent  set  de  *V0  k  50  ms-1  ,  La  vitesse  de  rotation  maximale  de  la  balance  est  de  2,5  tours  par  seconde, 
soit  150  rpm.  L'envergure  de  la  maquette  peut  attelndre  1,30  n. 


Rewaroue.-Les  axes  *0y_®0»  reprbeentbs  sur  la  figure  oi-dessus,  sont  libs  h  la  balance.  L'axe  x  ,  parallble 
b  la  direction  de  l'ecoulenent  est  fixe  et  horizontal,  done  confcndu  avec  l'axe  longitudinal  du°trlbdre 
normal  terrestre.  Four  que  les  autres  notations,  6  ,  V  ,  z  ,  y  solent  conformes  aux  rscommandations  ISO, 
TT~7aut  suppoeer  que  la  figure  est  reprbsentde  dans  la  position  oii  le  mSt  est  vertical,  et  le  calage  de 
l'angle  ®  effectub  b  l'azimut Y  -  0,  e'est-i-dire  x^  et  xh  confondus. 


Mf Ere  nee.-  I*  St*11  Behaviour  and  Spin  estimation  Method  by  Use  of  Rotating  Balance  Measurements 

Stall/Spin  Problems  of  Military  Aircraft.  AGARD  CP  199  nov.  1975. 


ANNEXE  A1I.1.2 


ROTATION  UNIFORME 

TRANSSONIOUC  SUFER80NIQUE  -  NASA  Aims 


1  -  axe  de  rotation  de  l'ensemble  dard-maquette, 

2  -  maquette , 

3  -  darxt  dynamome  triques ,  droit  et  eoude, 

4  -  dispositif  de  reglage  de  1' angle  9  , 

5  -  dispositif  de  reglage  de  l' angle  9,  , 

6  -  contrepolds  d'equilibrage  statique  de  la 

partie  toumante. 


7  -  axe,  parallele  a  l'axe  de  rotation, 

autour  duquel  la  position  du  contre- 
poids  est  ajust^e, 

8  -  moteur  hydraulique, 

9  -  conduites  d'arrivee  et  de  depart  du  fluids 

moteur  sous  pression 

10  -  support  fixe  de  la  balance 

11  -  mSt  porte-dard  de  la  soufflerie. 


souffleries  :  3,65  m  de  diametre  pressurisee  et  3,35  x  3,35  m 


transsonique 


0,*  $  1,* 


rotation  anirrtee  par  moteur  hydraulique,  deux  sens  de  rotation 


Q  s  *2  rad/s 


envergure  maximale  de  la  maquette  b  =  0,7m 

La  mlse  en  incidence  a  et  en  derapage  p  est  effectu^e  en  fonetion  des  angles  »  et  <!>,  dont  les 
valeurs  sont  fixdes  par  une  commande  Si  distance  sans  necessiter  de  do  montage  de  la  maquette* 


Angle  of  atteck  and  sideslip  envelope 
for  various  stlng-nodel  angles. 


Un  dard  droit  et  deux  dards  coud^s  permettent  d'opdrer  de  fag on  continue  dans  les  3  plages  d' incidence 

raprdaentdea  sur  la  figure  ci-dessus 


Rdfdrence  :  MALCOLM  O.M.  et  DAVIS  Sanford  S.  (Ames  Research  Center  NASA)  New  NASA  Ames  Wind  Tunnel 
Techniques  for  studying  Airplane  Spin  and  Two-dimensional  unsteady  Aerodynamics. 
DYNAMIC  STABILITY  PARAMETERS  N°  3  FDP  A SARD  235  1978. 


ANNEXE  AD. 2 


Coefficient  d 'nmortlsteaent  de  etructurc 
Evolution  en  cnplltude 


ORIGINE  t  VAUCHERET  X.  &  BROUSSAUD  P. 

Dispositif  d' oscillations  libres  pour  souffleries  &  rafales 
La  Recherche  Aerospatiale  n°  107,  p.  49-52,  1965 
et  ONERA  T.P.  n°  285  ( 1 965) 


ANNEXE  AH.4.1 


OSCILLATIONS  ENTRETENUES 


TRANSSONIQUE  OU  SUPERSONIQUE  -  R-A.E,  Bedford 


l - 1  m 


O 

1 


Pig.  All.4.1  •  Structure  de  la  balance 

1  maquette  • 

2  ferrure  d' attache 

o 

Les  f leches  indiquent  les  emplacements  des  jauges  de  mesure  dea  deformations. 


Fig.  AH.4.1  b  Principe  du  dispositif  d'entretien  des  oscillations  transversales  de  la  maquette  A®,  Ay, 

Ay,  . 

1)  bielle  centrde  sur  l'axe  longitudinal  G  *  ,  portant  k  ses  extrdmitds  les  bras  b,  et  bj  perpendi- 
culaires  au  plan  des  lames  Gx  I  4  la  position  moyenne  de  1' oscillation. 

2)  bielle  parallfele  au  plan  des  lames  disposee  4  la  distance  bj  de  ce  plan  et  inclinee  de  45°  sur  Gx  . 

3)  tringlerie  de  transmission  de  I'effort  imposd  Al;  par  l'oscillateur  electromagndtique  (4). 

La  detection  de  I'effort  transmis  par  (4)  peut  Stre  effectude  par  des  jauges  disposdes  sur  b,  ou  bt 
AF.  force  ddlivrde  par  l'oscillateur  dlectromagndtique. 

Conditions  de  fonctionnement 

masse  de  la  maquette  55  kg  intervalle  de  frequence  4^  f  ^10  Hz 

force  normale  maximale  4  000  N  taux  de  fatigue  des  lames  27,6  hb  (20  tons/inclr) 

amplitude  des  oscillations  Wz  =  i  2  deg. 

A*  =  *  2  deg. 

Ay  -  ±  4  mm 

Souffleries  :  RAE  subsonique  (3,9  x  2,7  m1 2 3),  supersonique  (2,4  x  2,4  m^)  et  A.R.A  transsonique  (2,7  x  2,4  m2) 

Rdf.  :  J.S.  THOMPSON  et  R.A.  FAIL  R.A.E.  Oscillatory  Derivative  Measurements  on  Sting-Mounted 
Wind  Tunnel  Models  A.O.A.R.D.  CP  n°  17,  1966 

FAIL  R.  R.A.E.  Experimental  Determination  of  Stability  and  Control  Aerodynamic 
Determination  V.K.I.  Aircraft  Stability  and  Control  Lecture  Series,  1975 


£TtA 


sections  du  d&rd 


3  lame  de  ressort  avant  ddformde  en  torsion  et  flexion 

4  lame  de  ressort  arri&re  ddformde  en  flexion 

5  bSti  support  dard 


ANNEXE  AO.4.2 


OSCILLATIONS  ENTRETENUES 
TRANSSONIQUE  ET  SUKRSONIQUE  -  NA.E. 


en  kl^vation 
coupe  du  dard 
par  le  plan  xz 


vue  en  plan 
coupe  du  dard 
par  le  plan  xy 


-  Balance  de  stability  du  N.A.G.  sur  maquette  montde  en  dard  et  animke  d’oscillations  entretenues 
par  un  moteur  dlectrodynamique. 

a)  vue  perspective  delatke  maquette  disposke  en  tangage 

b)  schdma  de  principe 


1  maquette  7 

2  ferrure  d'attache  be  la  maquette  8 

3  ressort  cruciforme  de  roulis, 

4  etrier  ' 

5  ressort  cruciforme  de  lacet 

6  ressort  A  lames  croisdes  de  tangage  (ou  de  ^ 

lacet) 


extrkmitk  avant  du  dard 

ressort  k  lames  parallfeles  donnant  un  accroisse- 
ment  de  rigidity  en  tangage  (ou  en  lacet) 

circuit  £lectrique  a.  enroulement  unique  parcouru 
par  un  courant  alternatif  k  frequence  rkglable 

aimant  permanent  en  deux  demi-coquilles 

dard 


OriginatORLIIt-RWlLUIAIW  K.J.  and  HARTf  I.S,  Exparlsenta  on  Cross-Coupling  and  Translational 

Aceslszstion  Dsrivativss  A CARD  CP  235  H*8  1978 


ANNEXE  AIL  5 


A2-7 


OSCILLATIONS  FORCCES 

TRANSSONIQUE  OU  SUFERSONIQUE  -  CA.LAF-AN. 


2 

Balance  de  stability  de  la  aoufflerie  transsonique  supersonique  (2,4  x  2,4  m  )  du  CAL3PAN. 

1 
2 

3 

4 


8  volants  fixes  sur  lea  axes  des  excentriques  Ei 

et  Ex 


maquette 

bloc  dynamombtre  et  accblbrombtre 

dard  oscillant  autour  de  l'axe  Pyo  normal  k 

1 '  eeoule  merit 

bielle  tranamettant  le  mouvement  de  l'excentrique 
E,  a  1' extremity  arribre  du  dard 


bielle  intermediaire,  oscillant  autour  de  «iy. 
fixe,  tranamettant  &  Py.  le  mouvement  de  l'excen¬ 
trique  Ex 

bfiti  fixe  lie  a  la  structure  de  la  aoufflerie, 
support  de  l'axe  fixe  (*y0  normal  k  l’dcoulement 
bielle  tranamettant  le  mouvement  de  1 1 excentrique  E* 


Le  rayon  et  la  phase  relative  de  l'excentrique  E ,  par  rapport  a  E2  sont  rbglables  a  distance.  Le  mouve¬ 
ment  de  la  maquette,  b  deux  degres  de  liberte,  consiste  en  une  rotation  autour  de  l’axe  Py„  normal  b 
l'bcoulement  et  a  une  rotation  de  cet  axe  autour  de  l'axe  t»y0  fixe  et  parallble  b  Py0  . 

Lea  conditions  de  fonctionnement  sont  : 

frequence  :  J  ^  f  ^  12  Hz  ;  acceleration  lineaire  maximale  (  -  20  g)  :  li’l  196  m/sec2  ; 
accblbration  angulaire  maximale  10)^200  r/ aec^  ;  charge  normale  maximale  IZ|  Eg  5  000  N  . 

Ces  conditions  permettent  d'obtenir  simultanement  a  une  frequence  de  7  Hz,  des  amplitudes  de  translation 
A  I  =  ±  100  mm  et  de  rotation  Af)=  ±  6°. 

Ce  dispnsitif  offre  la  possibility  d'impoaer  a  une  maquette  montee  sur  dard,  une  oscillation  d'amplitude 
importante  autour  d'un  axe  virtuel  situe  b  l'interieur  de  celle-ci,  ce  qui  perraet  d'obtenir  les  mesures  du 
torseur  des  forces  aerodynamiques  dans  de  bonnes  conditions  experimentales,  en  particulier  en  rendant 
minimales  les  valeurs  du  torseur  des  forces  d'inertie. 

Origine  :  STATLER  Irving  C.  TOFTS  Ovan  B.,  HIRTREITER  Walter  J. 

A  New  Capability  for  Measuring  Dynamic  Air  Loads  in  a  Wind  Tunnel 
Journal  of  Aircraft  Vol  3  N°  5  pages  418  b  425,  sept.  1966 
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comptes-rendus  des  reunions  FDP  et  FMP  AGARD,  comptes-rendus  des  reunions  FDP  et  FMP  AGARD, 

tenues  pendant  ces  quatre  demieres  annees.  Un  expose  tenues  pendant  ces  quatre  dernieres  annees.  Un  expose 

critique  des  methodes  experimentales  d’analyse  des  critique  des  methodes  experimentales  d’analyse  des 
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phenomenes  alrodynamiques  instationnaires  met  en  Evidence  les  lacunes  des  methodes  phenomenes  aerodynamiques  instationnaires  met  en  evidence  les  lacunes  des  methodes 
et  moyens  actuels  notamment  en  ce  qui  conceme  les  informations  necessaires  a  la  et  moyens  actuels  notamment  en  ce  qui  conceme  les  informations  necessaires  a  la 
mise  en  oeuvre  correcte  des  simulateurs  de  vol  dans  les  cas  de  comportements  non  line-  mise  en  oeuvre  correcte  des  simulateurs  de  vol  dans  les  cas  de  comportements  non  line- 
aires  (decollements,  ondes  de  choc).  Des  suggestions  sur  des  orientations  possibles  a  aires  (decollements,  ondes  de  choc).  Des  suggestions  sur  des  orientations  possibles  a 


